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Analysis of fouling thermal resistance of 
feed-water heaters in steam power plants 


Dariusz Butrymowicz, 
Jerzy Gtuch, 

Tomasz Hajduk, 
Marian Trela, 

Andrzej Gardzilewicz 


Abstract 


In this paper a problem is discussed of thermal degradation of shell-and-tube recuperative feed-water heaters due to heat 
transfer surface fouling. Application (to an example steam cycle of power plant) of DIAGAR software system intended for the 
analysing of impact of degradation of a particular heat exchanger on whole recuperative system performance, is presented. 
On the basis of the systematically performed simulative calculations it was concluded that the degradation of feed-water 
heaters of the lowest extraction pressure results in an additional load imposed on the heaters of higher extraction pressures, 
moreover it was estimated that the degradation of the heaters due to presence of sediments of the thermal resistance reaching 
5*10-* m?*K/W, has resulted in the drop of turbine set efficiency by about 0.3%. Such drop depends on a given 
configuration of degree of degradation of particular feed-water heaters applied in a considered turbine set. 


Keywords: steam power unit, fouling, feedwater heaters 


Introduction 


Power efficiency improvement of steam power 
generation systems is today one of the most crucial 
research problems and concerns both stationary and 
marine power plants. The problems associated with their 
operation including thermal degradation of heat power 
systems resulting from the fouling of their surface, 
become more and more important. 


In this paper is discussed the problem of impact of 
thermal degradation of heat exchangers working as a low- 
pressure heat recuperative system on power plant efficiency. 

An insufficient attention has been paid so far to the 
problem because it has been assumed that the problem 
of the fouling sedimented on heat transfer surface and 
thermal degradation associated with it, would not be very 
important due to the fact that the working medium of 
very high purity flows on both sides of such exchangers. 
As results from the authors' experience which has been 
gathered so far, the above presented opinion does not 
comply with real state ofthe matters, therefore the authors 
have initiated themselves a systematic research focused 
on the assessment of degree of thermal degradation of 
heat exchangers resulting from the fouling of heat transfer 
surface [1, 2]. In this paper the problem ofthe assessment 
of impact of surface fouling of low-pressure feed-water 
heaters on power plant efficiency, is discussed. The 
undertaken analysis deals with a supercritical power plant. 
Such choice is justified by the fact that in the power plant 
in question number of heaters forming heat recuperative 
system is the largest. 


The below presented analysis is based on calculations 
which have to be very exact in view of a high precision 
necessary to determine differences in efficiency of 
considered cycle configuration variants as well as when 
various component devices are applied. For the planned 
scope ofthe investigations to compare impact of thermal 
degradation of different feed-water heaters on energy 
efficiency of a considered cycle, is necessary. For such 
analysis two quantities are in particular important, namely: 


- specific heat consumption defined as the ratio of 
heat flux delivered to boiler's steam cycle and 
turbine's internal mechanical power; 


- turbine's mechanical power. 


The modeled object was thermally calculated and 
its characteristics were determined. To the calculations 
DIAGAR software was applied [3,4,5,6]. Some 
procedures of the program had to be so modernized as 
to make it possible to calculate the cycles in the range 
of supercritical steam parameters. Simultaneously, 
calculation procedures of heat exchangers were so 
modified as to make it possible to take into account set 
values of thermal resistance of the fouling sedimented 
on heat transfer surfaces. Values of thermal resistance 
were assumed on the basis of the experimental tests on 
thermal resistance of the fouling accumulated on pipes 
of one of Polish electric power plants, performed on the 
authors' original laboratory stand [7,8]. The work was 
aimed at determination of expected loss of efficiency 
resulting from fouling sedimented in such heaters. The 
loss was assessed in relation to the power plant of clean 
shell-and-tube feed-water heaters. 
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Fig. 1. Schematic diagram of the analyzed cycle of thermal power plant 
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Analyzed system of power plant 


In this paper the influence of thermal degradation 
of low-pressure feed-water heaters on an example 
thermal cycle of steam power plant, was analyzed. In 
order to solve the problem to apply a complete model 
of thermal cycle is necessary; this is however a very 
complex problem. Therefore it was decided to use 
a thermal cycle model of steam power plant, which was 
at disposal of these authors. Results of the analysis can be 
considered representative in a wide range of operational 
parameters of steam power plants, both for marine and 
land applications. 


In Fig. 1 schematic diagram of the considered example 
thermal cycle of 600 MW power, is presented. The 
remaining heat-and-flow data decisive of operational 
characteristics of the cycle are given below: 


- Live steam pressure: 28.5 MPa 
- Live steam temperature: 600 °C 

- Interstage superheating temperature: 620 °C 

- Condensation pressure: 5 kPa 

- Degassing pressure: 1.15 MPa 


The cycle is intended for co-operation with: 
- electric drive of feed-water pump; 

as well as 
- turbine drive of feed-water pump. 


Moreover from the heat-and-flow data the operational 
characteristics of the thermal cycle’s component devices 
such as: pressure losses, condensate overcooling, ram 
temperature rise etc, were calculated. 


Thermal degradation of heat exchangers 


Thermal degradation of heat exchangers is defined 
as the lowering of exchanger’s heat transfer efficiency, 
which results from operational factors such as: heat 
transfer surface fouling, inert gas presence, exclusion of 
a part of heat transfer surface from operation, or other 
operational causes. In this work special attention has 
been paid to the degradation caused by heat transfer 
surface fouling. 


Fouling thermal resistance depends on definite 
operational conditions of heat exchanger. With a view of 
different chemical composition of water, e.g. in different 
electric power stations, sediments of a different thickness, 
structure and chemical compostion can be met. Moreover, 
even in a given electric power station thickness and 
structure of sediments can be different as results from 
observations of state of heat transfer surfaces. Hence any 
generalization of only once measured values of fouling 
thermal resistance for all heat exchangers operating in 
a given electric power station, is very uncertain. 


Influence of thermal degradation resulting from the 
fouling of heat transfer surfaces as well as inert gas 
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presence necessitates one-dimensional model of heat 
exchangers to be applied. In such model, elaborated by 
Butrymowicz and Trela [1,2], a change in heat transfer 
conditions along water flow through exchanger is taken 
into account. With a view of necessity of taking into 
account changes of the conditions in particular runs, such 
model can be considered quasi-two-dimensional. In the 
model, impact of inert gas content in condensating steam 
as well as bundle effect are to be taken into account 
locally. With a view of inert gas presence to take into 
account local change of the saturation temperature T, 
as well as the phase separation surface temperature T,, 
is necessary. 


In the analyzed model the following balance 
differential equations have been used for calculations: 


' dT, , 
m,c Ga, DT) (r, -7,)-0 0) 


wt pw 


in which T,, stands for local temperature of water, and 
the heat transfer coefficient k is given by the following 
equation: 
$= l (2) 


1 : D 1 D 
ER AT. R R 2 R. 2 
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where the following notation is assumed: a, and 
a „— heat transfer coefficients on the side of condensating 
steam and cooling water, respectively; T,, -pipe outer 
surface temperature; D -bundle coefficient; T,. -steam 
saturation temperature; Rẹ, Rẹ R, -fouling thermal 
resistance of: outer surface, inner surface and pipe 
material, respectively; D,, D;, D,, -pipe outer, inner and 
logarithmic mean diameter, respectively; R, -additional 
thermal resistance due to inert gas presence. Thermal 
resistance of inert gas layer can be determined by using 
Berman relation [10] well experimentally confirmed and 
often applied to numerical calculations [9]. 


In the case of the analyzed model the following set of 
equations is considered: 


* the differential equation (1); 
e and the following algebraic equations: 


- the equation which combines the heat transfer 
coefficient on the condensating steam side, a, the 
temperatures T,,, Tp, and the steam mass flow 
rate m, [1,2], 


to? 


- the equation which expresses the thermal 
resistance, R,, resulting from inert gas presence [2], 


- the energy balance equation for an elementary 
segment of heat exchanger pipe: 


dm, lh, +C,, (T, zT )F mc, dT, (3) 


w pw 


- the equation which combines local steam 
saturation pressure and air content: 


4 
p,- of + 0.622 a.) e 
Tf, 


v 


where: m, — air flow mass rate; p, — steam partial 
pressure; h,, — specific evaporation enthalpy, c, — specific 
heat capacity of superheated steam. On the basis of Eq. (4) 
local value ofthe steam saturation temperature T, can be 
determined provided distribution ofthe air flow mass rate 
m, is known. Dependent variables are the following: the 
water temperature T,,, the steam mass flow rate m, , the 
phase separation surface temperature T,, the pipe outer 
surface temperature Tto as wellas the steam saturation 
temperature T,.. All the quantities are local. 


In the case of multi-run heat exchangers, heat 
transfer conditions on the condensating steam side 
undergo discrete changes, which is covered by the bundle 
coefficient B. Hence the solving of the set of equations 
in question is the most convenient if the differential 
equation (1) is written in the form ofa difference equation. 
In the case of occurrence of the convective cooling of 
superheated steam (without condensation) the considered 
set of equations requires to be extended, that is shown 
in [11]. 


The one-dimensional model of heat transfer within 
recuperative heat exchanger is deemed most suitable for 
diagnostic analysis purposes because it makes it possible 
to take into account local effects important from the 
heat transfer point of view, at a relatively low input of 
numerical calculations. 


a) 70 


0 10 20 30 40 5 0 80 90 10! 
L/L,,X100 [%] 
Fig. 2. Temperature distribution in XN1 heat exchanger along water path, 


Lw — the total water path length. In calculation it was assumed: 
a) Ry -L5 10°4 m?K/W; = 1*10? kg/s; b) Ry = 0; mz 0. 


In Fig. 2 are given example results of calculations of 
a recuperative heat exchanger absorbing steam from the 
last extraction of LP turbine of 200 MW steam turbine 
set, i.e. the feed-water heater VIII. The heat exchangers 
of the lowest extraction pressure are vertically installed 
in condenser casings, and their heat transfer surface 


is made by a bundle of horizontal pipes. The basic 
geometrical quantities of the example feed-water heater 
VIII for which heat transfer calculations were performed 
are the following: — number of pipes: 670; — pipe outer 
diameter: 19 mm; - pipe inner diameter: 17 mm; 
— pipe material: brass; — active length: 3.642 m; — active 
breadth: 0.358 m; — number of vertical pipe 
rows: 52; — number of water runs: 4; — cross-section 
area of steam inflow to heat exchanger: 2.548 mx0.200 m; 
— pipe triangular spacing: 44 mm. For numerical 
calculations the following operational parameters of 
the heat exchanger were assumed: — water mass flow 
rate: 63.74 kg/s; -water inlet temperature: 28.51?C; 
— steam inlet temperature: 93.06 °C; — steam inlet 
pressure: 0.271 bar. In the calculations it was assumed 
— by analogy to power condensers — that the steam mass 
flow discharged from deaerators is associated with the 
minimum local condensation temperature higher by 
2 +5 K than the water inlet temperature. 


From the point of view of diagnostic assessment 
of heat exchanger operation its thermal degradation 
is revealed directly by the increasing of the so called 
ram temperature rise, i.e. difference between saturation 
temperature and heated water temperature at outlet from 
the exchanger. 


As results from the presented diagrams, the distribution 
of the water temperature T, in the exchanger working in 
air-presence conditions differs greatly from logarithmic 
one characteristic for operation of a clean exchanger free 
from presence of air. The air content in condensating 
steam, amounting to m,= 1*10 kg/s which is 
equivalent to the ratio (m, /m,) = 0.31% assumed 
for the calculations, results in a significant increase 
of the ram temperature rise 67, =T —T,,,. For the 
non-degraded exchanger the increase of the ram 
temperature rise equal to 3,5 K is obtained, and in the 
case of a significant degradation due to the heat transfer 
surface fouling the increase by about 18 K is obtained. 
On the basis of other analyses [2] it can be stated that 
the degradation resulting from dirtiness of heat transfer 
surfaces only causes the increase up to 9.0 K. 


Tab. 1. Design data of low-pressure recuperative heat exchangers 


Values of geometrical characteristics 
of non-degraded recuperative heat 
exchangers 
Name of device Heat transfer | Number of pipes 
surface area in a heater 
[m2] [pieces] 
Recuperative heat exchanger VIII 435 704 — horizontal 
Recuperative heat exchanger VII 900 3183 — vertical 
Recuperative heat exchanger VI 1053 2874 — vertical 
Recuperative heat exchanger V 1060 3183 — vertical 


By taking into account the available research 
results on fouling thermal resistance [1,2,7,8,11], 
for purposes of this work two characteristic values 
of fouling thermal resistance can be proposed: 
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Fig. 3. Numerical scheme of 600 MW steam turbine unit with electrical 
drive of feed-water pump . 
Notation: 10 -boiler; 15 — pipeline with valves; 20 — control stage; 30 
— group of HP stages I; 40 -extraction I; 50 -group of HP stages II; 60 
-extraction II; 70 — interstage superheater; 80 -group of IP/LP stages I; 
90 -extraction III; 100 -group of IP/LP stages II; 110 -extraction IV; 120 
-group of IP/LP stages III; 130 -extraction V; 140 - group of IP/LP stages 
IV; 150 -extraction VI; 160 -group of IP/LP stages V; 170 -extraction 
VII; 180 -group of IP/LP stages VI; 190 -extraction VIII; 200 -group 
of IP/LP stages VII; 210 — outlet diffuser TG; 250 — main condenser; 
260 — condenser cooling; 270 — main condensate pump; 280 — degasifier 
steam condenser; 299 — pipeline of extraction VIII; 300 — recuperative 
heat exchanger VIII; 309 -pipeline of extraction VII; 310 -recuperative 
heat exchanger VII; 319 - condensate pump; 320 — water distribution; 
329 -pipeline of extraction VI; 330 -recuperative heat exchanger VI; 
339 -pipeline of extraction V; 340 -recuperative heat exchanger V; 

346 — pipline; 349 -pipeline of extraction IV; 350 — odgazowywacz; 360 
-pompa zasilajqca; 370 - water distribution; 371 — injection of water to 
interstage superheater; 380 -recuperative heat exchanger III; 399 -pipeline 
of extraction II; 400 -recuperative heat exchanger II; 409 -pipeline of 
extraction I; 410 -recuperative heat exchanger I; 429 -pipeline of extraction 
IIT; 500 — auxiliary turbine; 505 — auxiliary turbine diffuser; 510 — auxiliary 
condenser; 520 - auxiliary condenser cooling; 530 — condensate auxiliary 
pump; 540 -pipeline 
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e R, =2.5 10 ^m? -K/W ] -assumed to be typical 
for heat exchangers; 


eR; =5.0-10“[m* - K/W] -assumed to be 
relevant to heat exchangers of an elevated degree 
of degradation. 


In this work the two above mentioned values of 
thermal resistance were used in assuming various 
degradation degree configurations of low-pressure 
recuperative feed-water heaters. 


Calculations of the power plant thermal cycle 


The presented heat-and-flow calculations of the cycle 
were performed with the use of DIAGAR software 
[3,4,5,6]. Each of the considered cycle configuration 
necessitates the same numerical scheme to be applied. 
The numerical scheme of the cycle for the electric drive 
of feed-water pump is presented in Fig. 3. 


The estimated values of geometrical characteristics 
of feed-water heaters were based on the characteristics 
of the heaters installed in 200 MW power units, namely: 
the same inner and outer diameter and arrangement of 
pipes. The characteristic data of the heaters are given in 
Tab. 1. 


The collected results of the calculations are graphically 
presented in Fig. 4. There are also given the applied 
configurations of thermal resistance levels of the 
fouling in LP recuperative system. The results indicate 
the losses revealed by introducing to the calculations 
the heat transfer surface fouling of feed-water heaters, 


Relative change of quantity 
o 


Thermal degradation configuration No. 


oT relative hei umption 


Fig. 4. Changes of operational parameters of 600 MW power unit at 
various combinations of degree of degradation of low-pressure heat 
exchangers. Notation: 1 — clean exchangers; 2-moderately degraded heat 
exchanger VIII; 3 -moderately degraded heat exchanger VII; 4 -moderately 
degraded heat exchanger VI; 5 -moderately degraded heat exchanger V; 
6 — maximum degraded heat exchanger VIII; 7 -maximum degraded heat 
exchanger VII; 8 -maximum degraded heat exchanger VI; 9 -maximum 
degraded heat exchanger V; 10 -moderately degraded heat exchangers VIII 
— V; 11 -maximum degraded heat exchangers: VIII and VII; 12 -maximum 
degraded heat exchangers: VIII and VI; 13 -maximum degraded heat 
exchangers: VIII and V; 14 -maximum degraded heat exchangers: VIII, 
VII and VI; 15 -maximum degraded heat exchangers: VIII, VII and V; 16 
-maximum degraded heat exchangers: VIII — V. 


Tab. 2. Ram temperature rise of low-pressure recuperative heat exchangers 


Values of ram temperature rise [K] 
Name of device Re =0 ped n ics r i 
Recuperative heat exchanger VIII 1.62 3.55 5.72 
Recuperative heat exchanger VII | 0.63 3.82 7.88 
Recuperative heat exchanger VI | 0.73 4.08 7.92 


Changes of steam mass flow rate of heat 
exchanger VIII of SOOMW power unit. 
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Fig. 5. Changes of steam mass flow rate of heat exchanger VIII at various 


combinations of degrees of degradation of low-pressure heat exchangers; 
code of degradation -see Fig. 4 


Changes of steam mass flow rate of heat exchanger VII of SOOMW power unit 
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Relative change of steam mass flow rate in XN2 


Code of combinations of fouling of heat exchangers 
== - change of flow rate in XN2 
Fig. 6. Changes of steam mass flow rate of heat exchanger VII at various 


combinations of degrees of degradation of low-pressure heat exchangers; 
code of degradation -see Fig. 4 


Changes of steam mass flow rate of heat exchanger VI of 600MW power unit 
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Fig. 7. Changes of steam mass flow rate of heat exchanger VI at various 
combinations of degrees of degradation of low-pressure heat exchangers; 
code of degradation -see Fig. 4 


Changes of steam mass flow rate of heat exchanger V of 600MW power unit 
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Fig. 8. Changes of steam mass flow rate of heat exchanger VI at various 
combinations of degrees of degradation of low-pressure heat exchangers; 
code of degradation -see Fig. 4 


expressed in the form of relative changes in turbine-set 
mechanical power as well as relative changes in specific 
heat consumption. 


As the presented results yield, in the case of the 
fouling of a single feed-water heater almost the 
same level of specific heat consumption increase and 
turbine-set power decrease is obtained for all the 
feed-water heaters and the typical value of R,, (Case 
2, 3, 4 and 5), and for the elevated value of R,, certain 
diversity of the effects is observed since for the feed- 
water heater VIII (of the lowest extraction pressure) 
much lower changes are achieved. However, from 
the operational point of view the cases are more 
interesting of different configurations of the fouling 
of feed-water heaters, having a significant value of the 
thermal resistance Ry, - Case 10 through 13 -where for 
both the analyzed quantities the same level of changes 
was obtained. It allows to conclude that in the case of 
a significant degradation of the feed-water heater of the 
lowest extraction pressure, degradation of any of the 
remaining heaters of a higher extraction pressure does not 
play any important role. The greatest changes are obtained 
in the case 16 for degradation of three feed-water heaters. 
The changes in ram temperature rise in particular 
feed-water heaters are presented in Tab. 2. Along with the 
increasing of the fouling thermal resistance also the ram 
temperature rise increases, and in the conditions of the 
considered geometry especially high values of the ram 
temperature rise are obtained for feed-water heaters of 
higher extraction pressures. 


In Fig. 5 through 8 the degradation effects in the 
form of changes in steam mass flow from recuperative 
extractions to feed-water heaters are presented for 
various configurations of degradation in the heaters. 
From the obtained results it can be generally concluded 
that the fouling of the feed-water heaters of lower 
extraction pressure values results in the increasing of 
thermal load imposed on the heaters of higher extraction 
pressure values because mas flow of steam sucked-in 
by the heaters then increases. Thermal degradation of 
a feed-water heater certainly makes that a smaller steam 
mass flows to it. 


Heat transfer surface degradation of feed-water heaters 
impacts cycle's operational characteristics first of all by 
increasing the specific heat consumption. The calculation 
results given in Fig. 4 show that the highest impact on 
the worsening of cycle's efficiency (i.e. increasing the 
specific heat consumption) is introduced by the feed-water 
heater VIII, the first in the row of the recuperative feed- 
water system. It specially concerns moderate degradation 
levels, that leads already to necessity of cleaning them. 
This indicates that special attention should be paid just 
to cleanless of surface of the feed-water heater. Such 
conclusion is also confirmed by change in the steam mass 
flow passing through particular feed-water heaters, which 
constitutes one ofthe most important parameters deciding 
on performace of the exchanger itself. 
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Summary 


In this paper theoretical analysis of impact of thermal 
degradation of recuperative low-pressure feed-water 
heaters on change of specific fuel consumption as well 
as output power of turbine set, has been presented. 
Degradation of feed-water heaters of the lowest 
extraction pressure is the most important. The revealed 
level of the changes should be considered significant. 


In this stage of calculations of steam power plant 
cycles by means of DIAGAR software, carried out for 
the design load conditions and with taking into account 
cleanless of heat transfer surfaces of recuperative low- 
pressure feed-water heaters the following conclusions 
can be drawn: 


1. Increasing the ram temperature rise is a direct 
effect of presence of sediments on surface of recuperative 
feed-water heaters (characterized by certain value of the 
thermal resistance Rf). In the range up to Rf = 0.0005 
m2K/W, the resistance results in a multiple increase of 
ram temperature rise. 


2. The maximum degradation of all recuperative low- 
pressure feed-water heaters (at Rf = 0.0005 m2K/W) may 
resultinthe specificheatconsumption increase by about 0.3%. 

It is the value which can be deemed significant for 
operation of high-output power generation units. 


3. Thermal degradation of feed-water heaters of 
the lowest extraction pressure is accompanied by an 
additional thermal load imposed on the feed-water heaters 
of higher extraction pressures. Hence in the feed-water 
heaters of higher extraction pressures greater relative 
changes in sucked-in steam mass flow rate than in the 
heaters of lower extraction pressures, are observed. The 
relative changes in mass flow rate of steam sucked-in by 
the heater of the lowest extraction pressure do not exceed 
16%, whereas the changes for the heater of the highest 
extraction pressure are contained in the range up to about 
30%. The above described effects directly result in a drop 
of turbine-set efficiency, hence in a relevant increase of 
specific fuel consumption. 


Therefore it is necessary to fil up an important 
information gap in the area of real thermal resistance 
values resulting from the fouling of recuperative 
low-pressure feed-water heaters used in steam power 
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plants as the available data in question have not 
a comprehensive character. 


This scientific research project (the ordered project 
No. PBZ-MEIN-4/2/2006) has been financed from 
science budget resources for the years 2007-2009. 
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Specificity of automatic control of micro- 
turbines (steam or gas -driven and expanders) 
in dispersed generation system of heat and 
electric power 


Zygfryd Domachowski 


Abstract 
This paper presents specific problems of automatic control of steam micro-turbines and expanders intended for the dispersed, 
combined generating of heat and electric power. The investigations concern ensurance of certainty of energy supply and its 
required quality. 


Keywords: micro-turbine, heat-and-electric power plant, structures, automatic control. 


Introduction 


Dispersed electric power generation has many 
important advantages: it makes it possible to use local 
primary - energy sources (of water, wind, biomass or 
geothermal ones), it lowers electric power transmission 
losses, increases electric supply safety (in the case of state 
of emergency of large electric power system), improves 
stability conditions of electric power system (except 
of windmill electric generator farms). Application of 
biomass to electric power generation has an additional 
advantage: it facilitates the developing of coordinated 
economy, i.e. combined generation of heat and electric 
power. 


On the other hand the dispersed electric power 
generation introduces certain problems dealing with 
cooperation of dispersed electric power sources with 
an external (large) electric power system. Apart from 
technical requirements concerning the switching-on of 
such dispersed sources to an external electric power 
system, the problems result from necessity to ensure 
suitable quality of electric energy (its frequency and 
voltage). The bigger they are the greater power developed 
by such source and the greater variability of its power. For 
this reason the greatest operational (control) difficulties 
are produced by windmill electric generator farms — their 
nominal power reaches from a few dozen to several 
hundred MW, and power fluctuation may be significant. 
In view of required quality of frequency control of 
electric power system the windmill electric power 
generating plants requires to be aided by conventional 
ones driven by water, steam or gas turbines [4], or by 
diesel engine. 


Mini heat-and-electric power plants working on 
biomass, develop nominal electric power reaching from 
a few dozen, or even several, to a few hundred kW. 
Therefore from the point of view of frequency control 
their role in electric power system is incomparable with 
that of windmill electric generator farms. Nonetheless 
they should be also supported by conventional electric 
power sources due to the following reasons. Firstly, 
as a rule they do not take part in frequency control of 


electric power system since their controlled quantity 
is heat power, and electric power is an additional 
product whose quantity is associated with amount of 
produced heat. Secondly, total electric power of mini 
heat-and-electric power plants installed over a given 
territory, can be significant. Hence further considerations 
are focused only on selected variants of automatic control 
systems of mini heat-and-electric power plants and on 
specific operational aspects associated with them. 


Specificity of a seperate electric power system fed 
by a microturbine 


Mikroturbines can be used, and are used to delivering 
electricity and heat to: 


— dwelling houses, especially individual houses, 
— schools, 

— hospitals, 

— hotels, 

— office buildings, 

— small firms, 

— consumers of heat and electric power on ships. 


Insuch arole they are topics of various publications and 
patents, however it mainly concerns gas microturbines. 
From the point of view of electric power generation role 
of both types of microturbines (gas and steam driven) 
is analogous, therefore role of their automatic control 
systems is also the same. It makes that such considerations 
associated with a concept of automatic control of gas 
microturbine [2, 3, 6+13] can be also related to a steam 
microturbine. 


Steam microturbine first of all intended for heat 
delivering, delivers as much electric power as it results 
from amount of delivered heat. Therefore internal electric 
power system (internal electric network) fed by it, should 
be connected to an external electric power system in order 
- in the case of lack of balance between amount of electric 
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power produced by the microturbine and that absorbed 
by the internal electric power system -to supplement 
possible shortage from the external system or to deliver 
possible surplus to it, see Fig. 1. 


external electric 
power system 


electric power,P internal electric 


primary power system 


energy 


——— 


heat-and-electric 


power plant heat power system 


——— 


l 
l 
l 
| 
| 
heat, Q l 
| 
I 
| 
l fed object 


Fig.1. Simplified schematic diagram of feeding an object with heat and 
electric power 

And, in land applications possible periodical breaks of 
connection between internal and external electric power 
system due to various emergency situations should be 
taken into account. On ships the internal electric power 
system is the only one. Its connection to an external 
electric power system may be only considered in the 
case of ship’s stay in port. In such situation if amount 
of generated electric power (resulting from delivered 
heat amount) is smaller than the demanded (let’s say 
critical) for feeding e.g. a refrigerator, necessary lighting, 
pumps, computers, then steam pressure control at outlet 
from turbine should be replaced by electric power 
control (current frequency control). For this reason the 
U.S. Patent [9] shows various structures of connection 
between a micro-gas turbine intended for the feeding of 
computer network with electric power and internal and 
external electric power systems. 


In accordance with the concept described in the 
patent the internal electric power system is split into two 
parts: that unconditionally demanding to be fed (critical 
power demand) and the remaining part (non-critical 
power demand). The part of the internal electric power 
system, which is subjected to the critical feeding, is not 
connected to any external electric power system, but the 
remaining part of the internal electric power system is 
connected to it. Hence, in accordance with the concept 
presented in the patent in question in normal conditions 
the microturbine feeds separately (by separate converters) 
both the parts of the internal electric power system, see 
Fig. 2. In normal conditions to the internal system's part 
connected to an external electric power system only 
a surplus of electric power generated by the microturbine 
over that absorbed as critical demand, is transferred. 
Moreover, if the amount of electric power transferred 
by microturbine to that part ofthe internal electric power 
system is greater than the demanded, then the resulting 
surplus is directed to the external system; if it is smaller 
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Fig.2. General concept of connection between micro-turbine and internal 
and external electric power systems 
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the demanded shortage is absorbed from the external 
system. 


When internal electric power system is switched-off 
from external one then in the case of application of 
steam backpressure microturbine it is necessary to pass 
from steam pressure control at outlet from the turbine to 
electric power control (frequency). Schematic diagram 
of feeding the internal electric power system is presented 
in Fig. 3 (such concept is not covered by the patent in 
question. It concerns micro-gasturbines only). 
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Fig. 3. Concept of feeding the internal electric power system after its 
switching-off from the external system 


When microturbine is excluded from operation (for 
instance due to its failure, unserviceability, overhaul or 
repair) then the feeding of the internal electric power 
system is ensured by its connection to the external system, 
see Fig. 4. To the part of the internal power system, which 
constitutes non-critical load, electric power is transferred 
from the external system directly, and to the part which 
constitutes critical load — by means of a ,,non-critical” 
two-way converter (i.e that which allows for sending 
electric power in two directions). However, the direct 
feeding of both the parts of internal power system, 
i.e. without making use of any converters, is deemed more 
appropriate, see Fig. 5. 
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Fig.4. Concept of feeding the internal electric power system after switching- 
off the micro-turbine 
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Fig.5. Concept of feeding the internal electric power system directly from 
the external system after switching-off the micro-turbine 


Large objects can be fed, also for reliability 
reasons, by several microturbines operating in parallel. 
An example of such co-operation realized by four 
microturbines is presented in Fig. 6. 


Electric current parameters in seperate electric 
power system 


It is assumed that current parameters in the internal 
electric power system, there by parameters of electric 
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Fig. 6. Example of feeding the internal electric power system, by means of 
four micro-turbines 


power delivered by microturbine are the same as standard 
ones ina low-voltage electric power distribution system, 
1e. 50 Hz frequency (in Europe), and 230V voltage. 
However values of rotational speed of microturbine rotors 
are contained in the range of 45000+120000 rpm, that 
corresponds to electric current frequency in the range of 
750-2000 Hz. Therefore it is necessary to apply an 
electric power converter at output from microturbine- 
driven generator so as to adjust parameters of the current 
delivered by the microturbine (frequency, voltage) to 
appropriate quantities relevant to low-voltage electric 
power system. 


External 
electric 
power 


A concept of such converter is presented in Fig. 7. 
Alternating current of frequency in the range of 750-2000 
Hz, produced by microturbine -driven generator, is rectified 
(by rectifier as shown in Fig. 7). Direct current can 
feed consumers which require such feeding, as well as 
electric accumulators. The feeding of internal electric 
power system requires an inverter which converts direct 
current into alternating one of the parameters necessary 
for feeding the consumers in the internal electric power 
system (50 Hz frequency, 400V phase-to-phase voltage), 
to be applied. 


Purposefulness of direct current (DC) application 
to seperate electric power systems fed with electric 
power obtained from renewable sources, has been also 
considered [1]. Such concept is presented in Fig. 8. In 
seperate electric power system DC voltage amounts to 
about 800V. Such system is fed from DC voltage sources 


internal electric power external electric power 
generator rectifier inverter system system 
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Fig.7. Structure of connection of micro-turbine to internal and external 
electric power systems by means of converter and inverter 


micro-turbine 


(e.g. dispersed renewable sources) as well as -through 
rectifier -from external electric power system. An inverter 
is necessary to make delivering electric power from such 
seperate system to an external AC power system possible. 
For feeding AC consumers in the seperate internal DC 
power system individual inverters at their inputs are 
necessary. 


Ship electric power system is an example of seperate 
system. Application of direct current to such system is 
deemed worth paying attention. 


15 kV line b) 


a) 15 kV line 


Transformer 150.4 kV Rectifier 15/0.4 kV. 


0.4 KV line. 


Fig.8. A concept of replacement of alternating current (a) by direct current 
(b) in seperate internal electric power system 


Automatic control of heat generated by 
heat-and-electric power plant connected to 
external electric power system 


Automatic control of heat-and-electric power plant 
is aimed at ensuring generation of demanded amount of 
heat of determined parameters. Amount of simultaneously 
generated electric power is associated with generated 
heat amount, but not subjected to control. Nonetheless 
in such conditions it should be intended to generate as 
large amount of electric power as possible. Hence, if 
heat delivered to an object is used both to space heating 
and service-water warming-up then two steam flows 
of different temperature values should be controlled. 
Simplified schematic diagram of automatic control 
system of heat-and-electric power plant, relevant to 
such situation, is presented in Fig. 9. Then to apply 
a multi-stage -extraction -backpressure turbine would 
be appropriate. 


primary — GÓC 
energy > heat-and-electric Qi 
power plant Qe 


heat flow rate 
controller 


Fig. 9. Simplified schematic diagram of automatic control system of heat 
flow rate for space heating (Q1) and service-water warming-up (Q2) 


In the case of application of backpressure turbine, only 
heat flow rate at turbine outlet is subjected to control (by 
controlling outlet steam pressure). It must correspond 
to the higher value out of the two temperature values: 
of service water and of space heating medium. In the 
case of floor heating it is reduced to choice of heating 
water temperature, from which limitation for amount of 
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generated electric power results. For this reason if a heat 
accumulator is applied to warming-up service water then 
steam pressure at outlet from turbine could be controlled 
with regard to space heating medium temperature. It 
would be then necessary to increase periodically the 
pressure in order to supercharge the heat accummulator. 


When delivered heat is used only for service-water 
warming-up (in the periods when accommodations are 
not heated) then steam pressure at outlet from turbine 
has to correspond to a required value of service-water 
temperature. 


If also other renewable sources (solar collectors, 
geothermal water) were used for delivering heat to an 
object then heat-and-electric power plant will be put 
in operation to supercharge heat accumulators, only if 
necessary. 


Automatic control of heat generated by 
heat-and-electric power plant not connected to 
external electric power system 


In the case of lack of connection between internal 
electric power system and external one as well as in the 
case of a break of the existing connection (e.g. in states 
of emergency of external electric power system) the 
heat-and-electric power plant is expected to deliver both 
required amount of heat and required amount of electric 
power [5]. In such situation two variants corresponding 
to two types of turbine assumed to be used in a given 
heat-and-electric power plant, should be considered. 


Delivery control of heat and electric power if 
backpressure turbine is used 


In backpressure turbine only one quantity can be 
subjected to automatic control, i.e. only one signal can 
be controlled. It can be either flow rate of delivered 
heat — the controlled signal is then steam pressure at 
outlet from turbine, or electric power — the controlled 
signal is then rotational speed of turbine-set rotor. In the 
considered situation (lack of connection between internal 
and external electric power systems) turbine-set electric 
output power must be controlled in order to ensure the 
required feeding of electric power consumers in internal 
electric power system. Amount of heat delivered by the 
turbine depends on its output power only. 


If heat flow rate at outlet from the turbine is lower 
than the demanded (both for space heating and service- 
water warming-up) then the lacking heat amount is 
directly delivered from boiler by starting-up a reduction — 
and-cooling station (see the simplified schematic diagram 
of control of delivery system of heat and electric power, 
shown in Fig. 10). 


During heating period, steam pressure at outlet from 
the turbine should correspond to a demanded temperature 
of working medium (in the case of floor heating); in this 
case the controller R, controls the reduction-and-cooling 
station SRS,. Heat used to warming-up service water 
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boiler 


> 


to boiler 


Fig.10. Simplified schematic diagram of control of delivery system of 
heat and electric power in heat —and —electric power plant fitted with 
backpressure turbine. 

Notation: T — turbine, G — electric generator, OC — heat consumers, 
SRS — reduction-and-cooling stations, f — frequency, p, — steam pressure 
behind the turbine, t Q^ water temperature in heat accumulator, 

R — controller. 
would be then delivered by means of the reduction-and- 
cooling station SRS, which feeds heat accumulator, being 
switched-on periodically. Outside the heating period steam 
pressure at oulet from the turbine should correspond to 

a demanded temperature of service water. 


When heat amount resulting from quantityof 
generated electric power is greater than the demanded 
then steam surplus at outlet from the turbine is sent 
directly to boiler, bypassing the heat consumers, see Fig. 
10; in this case the controller Ry controls the reduction- 
and-cooling station SRS,. 


Delivery control of heat and electric power if 
extraction-and-backpressure turbine is used 


In the extraction—and—backpressure turbine two 
quantities can be controlled for the reason of possible 
impact on two kinds of control valves: before the turbine 
and behind the extraction. Therefore, apart from electric 
power, steam pressure in the extraction or behind the 
turbine can be also controlled. Fig. 11 presents a concept 
of control of steam pressure in the extraction when the 
steam absorbed from the extraction is used to warming-up 
service water, and the steam taken from turbine's outlet 
is used to space heating. Flow rate of the steam used for 
warming-up service water, is subjected to control.To 
this end, the valves behind the extraction are used. Flow 
rate ofthe steam at outlet from the turbine, used to space 
heating, is not controlled. 


When steam flow rate at outlet from the turbine 
is not sufficient to maintain temperature demanded 
for accommodations then the lacking amount of 
heat is absorbed from the extraction. To this end, the 
reduction-and-cooling station SRS,, controlled by the 
behind the turbine pressure controller R,,, see Fig.11. 


When steam flow rate at outlet from the turbine is 
greater than the demanded for space heating then the steam 
surplus is sent directly to boiler -the controller R, in Fig. 11. 


to boiler < Y < < 


Fig.11. Simplified schematic diagram of control of delivery system of heat 
and electric power in heat —and —electric power plant fitted with extraction- 
and-backpressure turbine of controlled flow rate of the steam used for 
warming-up service water. 

Notation: T — turbine, G — electric generator, OC 1^ heat consumer for 
space heating, OC, — heat consumer for service water warming-up, 
SRS — reduction-and-cooling stations, f — frequency, pp — steam pressure 
behind the turbine, p, -Steam pressure in the extraction, 

R — pressure controller. 

In another variant the control valves behind the 
extraction can be applied to control steam flow rate at 
outlet from the turbine, used for space heating. Then, 
steam pressure at outlet from the turbine serves as the 
controlled signal. Such concept of delivery control of 


heat and electric power is presented in Fig. 12. 


Steam pressure behind the extraction is subjected to 
control. Its change resulting from a change of flow rate 
ofthe steam used for space heating, triggers off action of 
the controller, Rus onto control valves behind the turbine 
extraction. The steam used for warming-up service water 
is periodically taken from the turbine extraction. The 
service -water temperature controller R, is applied to this 
end, see Fig.12. 
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Fig.12. Simplified schematic diagram of delivery control of heat and 
electric power in heat-and —electric power plant fitted with extraction-and- 
backpressure turbine of controlled flow rate of the steam used for space 
heating. 

Notation: t ue service water temperature in accumulator, 

R, - temperature controller; as to the remaining symbols — see Fig. 11. 


Outside heating period, steam from turbine 
outlet should be directed to service-water heater 
when turbine extraction is shut-off; it corresponds 
to the backpressure turbine control whose concept 
has been presented in the preceding section. 


Final remarks 


Application of steam microturbines to the dispersed, 
combined generating system of heat and electric power, 
bringsmanyspecificproblemsintheareaofautomatic control. 
From this point of view various operational situations both 
dealing with microturbine set connection to an external 
electric power system and with its operation in a separate 
electric power system, have been reviewed. In both 
the above mentioned variants structure of microturbine 
automatic control system must ensure feeding the so 
called critical consumers with electric power. À concept 
of microturbine's automatic control system which fulfils 
the requirement has been presented in this paper. Structure 
of automatic control system for microturbines, both of 
backpressure and extraction-backpressure type, intended 
for the feeding heat consumers oftwo types and different 
steam pressure values, has been proposed and discussed. 
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Polish aspects in the modern design solutions 
of steam turbine flow systems 
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Abstract 
This paper presents new solutions of turbine stages which have been applied to high-power steam turbine flow systems. The 
structures, elaborated in Poland, were designed with the use of the authors’ original computational software and verified by 
measurements performed on real turbine sets. 


Keywords: steam turbine, power cycles, flow systems, efficiency, errosion. 


Introduction 


Most steam turbines operating in Poland were 
produced by Zamech Works, Elblag. For power industry 
needs domestic solutions were elaborated and also 
licence designs, mainly from Russia, were adapted. 
For past years it was managed to elaborate a relatively 
modern production process based on own computer 
design systems. The programs have been based first 
of all on co-operation with Polish scientific research 
institutes. It would appear that after taking over the 
Zamech Turbines Division by the ABB, ALSTOM, 
the worldwide recognized concerns, realization of 
Polish designs will be eliminated. In line with the new 
philosophy the works in Elblag took over immediately the 
most modern design and technological solutions of the 
concern and started to use foreign numerical calculation 
programs. However the fact of keeping Polish engineers 
still at work in the enterprise, has resulted in that Polish 
turbine solutions have been still coming into being. 
In this paper are described two design solutions elaborated 
in co-operation with the Fluid Flow Machinery Institute, 
Polish Academy of Sciences, Gdansk. They deal with 
new turbine stages: high-loaded ones and those before 
regenerative extraction. 


Examples of realization of new stages of 
steam turbine 


High-loaded stages 


Impulse stages intended for operating under high 
loads, which ALSTOM, Elblag, has applied to 7CK-65 
turbine set in the Zielona Gora Electric Power Station, 
can be deemed the best example of the modern solutions 
presently realized in Poland. Their superiority over other 
solutions of the concern has consisted in that the blading 
system of the turbine of 65 MW electric power and 135 MW 
heat power, could be accommodated in one common 
cylinder, that is schematically shown in Fig. 1 [1]. 
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Fig.1 Machinery devices of 7CK-65 steam turbine unit installed in Zielona 
Gora Electric Power Station [1]. 


It was achieved at maintaining its high efficiency, 
improved availability, and obviously, material savings. 
The minimized number of 23 stages resulted from the 
application of new profiles permitting to operate under 
high loads, to majority of rings of guide vane and rotor 
blades. Geometry of the profiles had to be changed the 
higher loads to be realized. They had to guarantee a little 
changed kinematics of steam flow as compared with 
conventional solutions. This results from a change of inlet 
and outlet angles as shown in Fig. 2, acc. Yamamoto [2]. 


As can be observed in Fig. 2 the rotor blade profiles 
of highly loaded stages make flow passages longer, 


(CEG 


turning angle -113.5deg 


(a) conventional turbine cascade ] (b) ultra-highly loaded cascade 


Fig. 2 Comparison of rotor blade ring cascades: a) conventional and b) 
highly loaded, acc. [2]. 


a) b) 


Fig.3 Geometry of guide vane and rotor blade rings of the turbine stages: 
a) conventional and b) high-loaded 
that is resulted from a distinctly larger change of flow 
turning angle in blade rings. The so large changes of flow 
direction reaching even 160° are realized first of all due 
to change in suction (concave) side contours of profiles, 
as described in detail by Amos [3]. 


The larger flow turning angles result obviously in 
a greater loss of energy. In the tests performed on the 
new designed stages the turning angles were limited to 
the range of 130+140°. In this case the guide vane ring 
profiles operate at outlet angles in the range of 10+12° 
and the rotor ring profiles — in the range of 16+18°. 
Forms of the new profiles resulted from optimization of 
conventional impulse profiles elaborated in ALSTOM, 
Elblag, in co-operation with the von Karman's Institute, 
Brussels. Their geometry is schematically shown in 
Fig. 3. 


The performed numerical analyses confirmed that 
the new profiles have been applicable to the assigned 
aims [4]. The relevant flow calculations were carried 
out by using the CFD FlowER program and 3D model 
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Fig.4. Comparison of energy losses in function of the angle al for 
cylindrical blade ring cascades of; a) PKH3-40 profiles and b) PLK-50 
profiles, acc. [4]. 
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Fig.5. Comparison of energy losses for P2-35 and P2V-35 profiles in 
function of the inlet angle beta 1, acc. [4] 


which takes into account viscosity of medium [5]. The 
calculations dealt with single ring of guide vane and 
rotor blades for geometry and operational conditions 
appropriate for operation in heat power turbine. It was 
observed that the losses within the cascades of guide 
vane blades of the new, PKH3 profiles, at the small 
values of the outlet angle a, ~ 11°+13° and the velocities 
corresponding with higher loads, were lower by more than 
1%, as compared with those recorded in the PLK typical 
impulse blade rings, Fig. 4. 


In the PKH3-profile blade rings the lower level 
of flow losses resulted not only from the performed 
efficiency-orientated optimization of form of the 
profiles but also change of their geometry. In the new, 
thicker” guide vane blade profiles it was possible to 
limit breadth of the ring at maintaining similar values of 
strength parameters. Such operation leads additionally 
to limitation of secondary and edge losses at the same 
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Fig.6 Change of efficiency in function of loads for the group of three 
cylindrical stages with PKH-3/P2V and PLK-P2 profiles at different values 
of the inlet angle, a1, acc. [4]. 


POLISH MARITIME RESEARCH, Special issue 2009/51 15 


length of blades. It should be stressed that values of the 
optimum-efficient relative pitch of the new-profile rings 
exceeded 0.8, whereas values of the conventional relative 
pitch of guide vane rings were kept on the level of 0.7. It 
corresponds with a lower number of blades. 


A little lower losses than those for P2 conventional 
profiles were also obtained for new rings of rotor 
blades of P2V profiles. The lower losses concerned 
flows at smaller inlet and outlet angles, characteristic 
for higher values of stage load. Comparison of the new 
and old profiles showed that the lowering of losses 
was especially distinct for values of the inlet angle 
B,< 27°, that corresponded with values of the outlet 
angle B,< 18°. The observed differences exceeded 2% 
in favour of the new profiles, however the level of 
minimum losses for greater angles corresponding with 
lower loads on conventional rings, was lower by 0.5%, 
(see Fig. 5). 


The positive results of analysis of steam flow through 
rings of PKH3/P2V profiles encouraged to perform 
successive numerical tests. They dealt with investigations 
of turbine stages. In this case was searched for an optimum 
efficiency-orientated geometry of three cylindrical stages 
fitted with the new and conventional profiles, suitable 
for operational conditions of the first group of stages 
of the designed heat power turbine. It was revealed that 
the stages with PKH3 P2V profiles were more efficient 
than those with PLK P2 profiles at the same assumed 
diameter D,, = 800 mm, in the range of greater loads at 
the applied angle a, — 11.5? (see Fig.6). In the diagram 
the loads are determined by the ratio u/co. The curve of 
efficiency change in the new stages appeared more ,,flat” 
as compared with the characteristics of typical impulse 
stages. The superior features are not observed at greater 
values of inlet angles, i.e. on the level of 13.59. 


The relative pitch values of the calculated stages 
at which the highest efficiency was obtained, was 
determined on the level of 0.8 for the new, PKH3 guide 
vanes, and for the old, PLK ones on 0.7, at the similar 
pitch values for both rotors, on the level of 0.7. It resulted 
from the optimization process of mutual orientation of 
blades. 
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Fig.7 CAD system for high-loaded turbine stages, acc. [6]. 
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Also in this case a favourable impact on efficiency was 
given by the smaller breadth assumed for the new guide 
vane profiles, that additionally made length of turbine 
flow passage, shorter. As above mentioned, the turbine 
stages with the new blade rings which permit higher 
loads to be applied, were practically installed in 7CK-65 
turbine construction. Relevant design calculations were 
carried out with the use of the CAD system elaborated in 
Poland, whose schematic diagram is shown in Fig. 7 [6]. 


Data concerning thermodynamical parameters of the 
stages were determined from thermal cycle calculations 
of the turbine unit, with taking into account parameters 
of boiler, water cooling system, as well as the electric 
and heat power assumed for the unit. From preliminary 
calculations of the turbine (OD model) it has resulted 
that cylindrical stages could be used for three groups of 
stages. The following basic gabarites were determined 
for them: diameters, blade heights, as well as number of 
stages. The data were used as the basis for detail design 
calculations (synthesis) and check calculations (analysis) 
which were rperformed by means of simplified programs 
for 1D and 2D model. In the programs solutions of 
non-viscous-flow equations were corrected by means of 
coefficients obtained from CFD numerical calculations 
for blade rings with new profiles, in accordance with the 
concept presented in [7]. 


Efficiency optimization procedures built-in the 
programs, made it possible to determine passage 
geometry, as well as number of blades and their most 
favourable orientation in particular stage rings. In the 
calculations the producer’s design standards concerning 
diameters, overlaps as well as solutions of sealing were 
taken into account [8]. 


The increasing of stage loads made it possible to 
reduce number of cylindrical stages in the turbine from 20 
to 16. It was achieved at a high flow efficiency confirmed 
by numerical calculations of 3D model [9]. As obtained 
from the tests, the efficiency was -by a little more than 
1% -different from the maximum value available in the 
case of conventional impulse stages. 
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Fig.8 Axial cross-section of 1st group blade system of 7CK-65 heat 
power turbine, acc. [8] 


In Fig. 8 the computer drawing of the flow system of 
Ist group cylindrical stages is presented. 


The drawing was automatically generated on the 
basis of heat-and-flow calculations just after design 


verification of particular rings by means of strength and i 
dynamics calculations [8]. As a result of the lowering 

of number of cylindrical stages it was possible to resign 

easily from an additional cylinder originally planned 

in the turbine. A photo of the turbine set and its axial 
cross-section covering also two steam -flow control stages 

intended for heating the extractions is shown in Fig. 9. 


The calculations for the designed flow system of the 
turbine were verified experimentally in the heat-and- 
power station. To this end the routine guarantee tests 
of the entire turbine set were supplemented by thermal 


measurements in compliance with the concept given 
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Fig.9 Photo and axial cross-section of 7CK-65 steam turbine set, [4]. 
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Fig.10 Schematic diagram of the measuring system intended for the testing 
of the cylindrical stages of 7CK-65 turbine, acc.[11]. 
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Fig.11 Run of steam expansion line for 7 CK-65 turbine, acc. [10]. 


in [10]. Schematic diagram of the measuring system 
intended for the testing of cylindrical stages is shown in 
Fig. 10, acc. [11]. 


The measured pressure and temperature values made it 
possible to determine steam enthalpy and expansion line 
in this part of the turbine (Fig. 11). The tests covered series 
of four measurements in various operational conditions. 
The performed analysis showed that the measurement 
errors in estimating efficiency of two first groups of the 
turbine’s stages did not exceed 1%. The third group of the 
turbine’s stages was excluded from the analysis because 
during the tests its operational parameters appeared 
greatly variable. It was caused by the fact that to this 
group the steam leaking from seals, largely changing the 
blade system parameters, was delivered. 


Tab. 1 presents comparison of calculated and measured 
values for 1st groups of cylindrical stages in nominal 
conditions. 


Tab. 1 Comparison of calculated and measured values for Ist groups of 
cylindrical stages in nominal conditions 


As measured As calculated 


Inlet pressure [MPa] 6.287 6.287 


Inlet temperature [oC] 498 498 
Outlet pressure [MPa] 1.553 1.553 
Outlet temperature [oC] 310.1 307.5 


Efficiency of the group of stages [-] 0.879 0.892 


Mass flow rate [kg/s] 51.05 51.35 


Only 1.5% differences in efficiency were obtained at 
a very high conformity of the calculated mass flow rate 
with that measured with the use of an area reducer. They 
were only a little smaller from measurement errors. As the 
analysis showed, the higher value of calculated efficiency 
resulted from the fact that non-stationary processes were 
not taken into account [9] during determination of the 
correction coefficients used in the design programs. 
The 7CK-65 one-cylinder turbine of 16 stages fitted 
with the new profiles has faultlessly operated for over 
seven years. 
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New solutions of stages before extraction 


In construction of turbine stages to introduce clearance 
(gaps) between moving and motionless elements in order 
to ensure their co-operation is necessary. Steam fluxes 
flowing through such gaps do not transfer work to rotor, 
which results in loss of energy. The phenomenon is 
especially intensive in low-pressure parts of turbine 
where large gaps over unshrouded rotor blades are 
present and fast flowing steam makes medium outflow 
to regenerative extraction difficult. In the blade system 
whirl and stagnation zones are then formed, that was 
confirmed by probing measurements carried out on real 
turbines (see Fig. 12), acc. [12]. In order to eliminate the 
intensive fluxes it was proposed to introduce, to the stages 
placed before regenerative extraction, a ring which will 
direct the leaking steam to the heat exchanger directly, 
Fig.13. The relevant solution was claimed to Polish patent 
office [13] and sucessfully implemented in LP parts of 
200 MW turbines with Bauman's stages [14]. 
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Fig. 12 Low-pressure part of 200 MW turbine with measurement 
instrumentation. Separation zones in leakage area are shown in grey. wg [12]. 


Advantages of application of the ring, with a view 
of turbine efficiency and operation, consisted in: 


* elimination of whirl zones in the stage behind 
extraction both in meridional and circumferential 
plane; 


* elimination of the mixing of two fluxes: of main 
flow and leaking flow; 


*the making use of higher energy of leaking steam 
in the first, usually underheated, regenerative 
extraction; 


*the lowering of moistness losses and errosion 
damage by effective separation of water drops to 
the extraction; its low cost as compared with other 
solutions. 
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Fig.13 New solution of turbine stage before extraction, acc. 
Patent no. P 160-805 


Both the performed theoretical investigations 
and measurements have demonstrated superiority of 
application of the new solution. Owing to it, shrouding 
seals -much more expensive and difficult in the case of 
realization of long blades -could be eliminated. 


Fig.14 Photos of 200MW turbine Bauman’ stage in versions before (left) 
and after (right) its modernization, acc. [14]. 


Photos of the turbine blade system in the versions 
before and after introduction of the ring are presented in 
Fig. 14. Salt sediments in the stagnation zones, eliminated 
after introduction of the ring, can be there observed. 


In the years 1994-2006 the solution was implemented 
in fourty 200 MW turbines in Poland and Ukraine. The 
average power increase of 600 kW/unit was obtained. The 
rings operated in the turbines faultlessly, that was confirmed 
by the observations carried out during operation [15]. 


The successful implementations have contributed to 
application of the considered solution also to the most 
modern design of LP parts fitted with ND 41AALSTOM 
outlet, used for modernization of 200 MW turbines. In 
this case the decision on the application has been taken 
not so much due to profits in efficiency as to decrease 
faster erosion of inlet edge of blades of the last rotor 
stage, spread over 2/3 of their length (see Fig.15). This 
region was reached by water drops seperated in the next- 
to-last stage and intensively splashed by leakage flux. 
The recorded erosion damage are specially dangerous in 
unhardened areas shown in Fig. 15, acc. [16]. 


Fig.15 Photos of inlet edge erosion of the blade of the last rotor stage after 
40 000 h service, acc. [16]. 


Fig. 16 presents the concept of introduction of the 
rings to the modernized LP part of 200 MW turbines, 


acc. [17]. 


extraction 
chamber 


Fig.16 Concept of introduction of the ring to the LP part of 200 MW turbine 
modernized in ALSTOM Works, acc. [17] 


After performing some additional analyses the patent 
was implemented to the power unit no.5 of Kozienice 
Electric Power Station. The analyses concerned detail 
heat-and-flow calculations in which values of its 
operational parameters before and after modernization 
were compared to each other. The CFD numerical 
calculations were carried out by means of the FlowER 
software [5]. For both the cases identical steam flow 
boundary conditions which resulted from the probing 
measurements carried out on the real power unit [18], 
were assumed. They covered distributions of static 
pressure, total pressure, temperature and steam flow 
angles before and behind stages, for a set value of steam 
mass flow rate. Description of the applied calculation 
methods can be found in [19]. 


The calculated total pressure distributions and 
streamlines prior assembling the rings in the diffuser 
between 3rd and 4th stage are presented in Fig. 17. The 
images distinctly show how the leakage blocks the steam 
flowing toward extraction. 


In the performed analysis concerning the modernization 
the ring’s form most favourable for efficiency was 
determined by using multi-variant numerical calculations. 


The following quantites were optimized: diameter of 
the pipe the ring was made of, its length, inlet edge 
orientation angle and its location with respect to rotor 
blade. The location was crucial for final choice as it had 
to guarantee that any intensive leakage is eliminated 
from the main steam flux at only minor disturbance of 
it. Some characteristic images of the distributions of 
flow parameters for the calculated variants, are given in 
Fig. 18. 


Total pressure distribution and streamlines in the 
diffuser and extraction chamber for a selected variant are 
presented in Fig. 19. The leakage flows to the extraction 
chamber and after being braked it rather smoothly splits 
into the flow to the heat exchanger and the continued flow 
through the blade system. For the ring’s position selected 
for nominal conditions, 3 kPa rise of pressure in the 
extraction chamber with respect to the initial version, was 
recorded. Also, the inflow to the guide vanes ring of the 
last stage has been favourably ordered, that is illustrated 
by the uniformly arranged stremlines. As calculated, the 
phenomenon resulted in the stage efficiency increase by 
over 1%. The shift of the ring inlet from the edge of the 
rotor blade should ensure effective sucking-off the water 
accummulated in this zone due to action of centrifugal 
forces, and to direct it into regenerative extraction and 
reduce erosion this way. 


The performed calculations made it possible to assess 
forces acting on the ring, which were the basis for strength 
and dynamic verification. The recorded loads appeared 


m 


Fig. 17. Total pressure distrbution (lefi) and streamlines(right) in the 
interstage diffuser, close to regenerative extraction, prior assembling the 
ring, acc. [18]. 
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Fig. 18. Distributions of steam flow total pressure close to extraction in 
turbine for selected calculation variants. 


Fig.19. Distrbutions of total pressure (left) and streamlines(right) in the 
interstage diffuser, close to regenerative extraction, afier assembling the 
ring, acc. [18]. 
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rather small and danger of resonances effectively 
eliminated by proper choice of number of ribs fastening 
the rings to clamp [20], which was taken into account in 
the final design. 


The significant changes in extraction chamber pressure 
and the last stage efficiency were used to assess energy 
profits to be obtained from application of the patent in 
question. They were assessed by calculations of thermal 
balance of the entire turbine set, performed with the use of 
the thoroughly tested DIAGAR software [21]. Main data 
for the above mentioned calculations were taken from 
measurements of thermodynamic parameters prepared 
on the basis of the guarantee tests [22]. 


Schematic diagram of thermal cycle of the power unit, 
i.e. 200 MW power unit no.5 at Kozienice Electric Power 
Station, is presented together with depicted measurement 
points in Fig. 20. 


STEAM GENERATOR 


Fig. 20. Schematic diagram of heat-and-flow cycle of the power unit, i.e. 
200 MW power unit No.5 at Kozienice Electric Power Station, with depicted 
measurement points, acc. [22]. 


Changes of output power and specific heat consumption 
in the cycle after implementation of the patent to the 
modernized turbine fitted with ND41 outlet in function 
of output power and condenser pressure, are shown in 
Fig. 21. 


Changes of output power of the power unit no. 5 
after implementation of the patent. 
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Fig.21 Changes of output power and specific heat consumption for various 
loads of the turbine set after implementation of the patent in question, acc. 


[19]. 


As can be observed in the figure, within the turbine load 
range of 120-220 MW and the condenser pressure range 
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of 3+6 kPa the decrease of specific heat consumption 
amounts to 11-16 kJ/kWh, which results in the output 
power increase of 200—400 kW. The equivalent increase 
of the turbine-set efficiency is in the range of 0.1+0.2% 
which can not be measured during routine guarantee tests. 


The efficiency rise effects can be only determined 
indirectly during measurements after modernization. 
It concerns the easily measurable rise of feed-water 
temperature in the first heat exchanger, by 39C. 
Measurements of distribution of thermodynamical 
parameters can be used in interstage spaces, which can be 
realized by the probing method only. Such measurements 
are planned to be performed soon on the turbine set no. 5 
at Kozienice Electric Power Station (see Fig. 22). Within 
the scope ofthe measurements control of erosion damage 
growth in the last stage rotor blades as well as inlet area 
of the rings, is provided for. 


» 4 4 
Fig.22. The ring placed inside 225 MW turbine 
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Selected problems of determining an efficient 
operation standard in contemporary 
heat-and-flow diagnostics 


Jerzy Gtuch 


Abstract 
In this paper are specified causes of ambiguities and difficulties in determining an efficient operation standard for land and marine 
turbine power plants. Two of them, namely: lowering values of working media inlet parameters and changes in configuration 
of connections between component devices of thermal cycles, have been analyzed. As a remedy for resulting difficulties 
a thermal cycle computational model of modular structure has been proposed. As showed, the tuning of the computational 
model with respect to correct measurements by elaborating individual corrections for calculation coefficients of computational 
methods used for particular devices ensures good accuracy of determination of unserviceability symptoms in these conditions. 
Introduction to the modular computational model — of control modules making use of valve opening indicators in DCS system 
ensures good accuracy also in the case of changes in structure of complex power object. 


Keywords: steam turbines, operation of turbines, power units, thermal measurements, thermal diagnostics, diagnostic relations 


Introduction 


Heat-and-flow diagnostics belongs to such forms 
of process diagnostics, in which the diagnostic method 
based on application ofa standard, turns out to be useful, 
Fig. 1. As it belongs to process diagnostics its main task is 
to assess efficient run of energy transformation in turbine 
thermal cycles [2, 5, 8, 10, 12, 13, 15, 21]. Such standard 
is that of efficient work of operationally undegraded 
object [5, 10, 15]. 
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Fig. 1. Schematic diagram of diagnosing process based on a standard 
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Heat-and-flow diagnostics is intended for the 
determining of efficient operation of power plants 
fitted with thermal turbines. Turbine engines are applied 
to land and marine power plants. In the case of steam 
turbine power plant, its equipment with many measuring 
instruments is characteristic both in land applications, 
[4,7,9, 10, 11, 14, 16, 20], Fig. 2, and marine ones [3, 6], 
Fig. 3. The equipment is available in the frame of 
Distributed Control System (DCS) [1, 14, 16, 17, 18, 19]. 
It makes good prospects for diagnosing objects of both 
kinds by using an efficient operation standard. 


Inefficient work of steam power plant results from 
operational degradation of geometry of component 
devices of thermal cycle. Control of efficient run of energy 
transformations in power plant constists in examining 
signatures consisted of inefficient work symptoms. 
Every symptom determines deviation of value of 
heat-and-flow parameter corresponding to it (mass flow, 
pressure, temperature), or characteristics (e.g. efficiency, 
steam flow capacity of particular component elements) 
from a reference value characteristic for undegraded 
object. 
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The above presented description seems to be simple 
but only in the case of object working in stable conditions 
and with its unchangeable structure. Yet steam turbine 
power plant (both land and marine one) is characteristic of: 


* varying conditions of the loading of the power 
plant, and of its environment, 


* dependence of power plant operational parameters 
on changes in its loading and environment, 


* dependence of power unit's operational parameters 
on operational degradation of geometry of 
component devices. 


It leads to difficulties in making diagnosis. 
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Fig. 2. An example of thermal cycle configuration and arrangement of 
measurement instruments in a high -power unit of an electric power station 
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Fig. 3. An example of arrangement of control -measurement instruments 
(represented by colour circles) in steam power plant of one of 
contemporary steam turbine ships [6] 


To obtain an unambiguous diagnosis, symptoms and 
signatures should be so formed as to get them depending 
on all characteristic states of operation. As acquisition 
of measurement results is performed with taking into 
account all the states it means that the same requirement 
is imposed on the standard. Therefore the standard-based 
diagnostics requires to have: 


e an efficient operation standard which reflects 
dependence of changes of parameters on current 
load state, 


e an efficient operation standard which reflects 
structure of a given object. 


Such aims are fulfilled by a standard of functional type, 
having modular structure, which makes it possible: 


*to identify measurement points, 
*to imitate (reproduce) object's structure and, 


eto calculate distribution of measured parameters in 
function of power unit's load and environmental 
conditions. 


Hence every change in mode of operation leads 
to determination of a new reference state, i.e. a new 
standard. 


The first ofthe considered situations which complicate 
diagnosis making, arises when periodical limitations 
of power unit operational parameters occur or when 
configuration of connections between the unit's devices 
should be periodically changed. The second situation 
takes place when a fault measurement gauge should be 
replaced by a new one. 


The third case concerns a machine in which only some 
of its degraded component devices were repaired. 


A seperate problem is to elaborate a standard, 
i.e. a reference state for on-line diagnostics. This 
paper is aimed at analyzing the first of the considered 
disturbances in operation of steam turbine plants. 
In the case of land power plants limitations of 
parameters may occur if there is a need to operate 
a non-fully serviceable object provided it is approved 
by its manufacturer. In practice are more and more 
often met circumstances to change configuration of 
feed component devices, which occur in the case of 
significant power limitation of the power unit. In the 
case of ship power plants the necessity to continue ship 
voyage under limited steam boiler parameters sometimes 
appears. In some cases it is necessary to keep ship 
propulsion system working at delivery of power from 
only one turbine cylinder. For these considerations 
200 MW power unit cycle whose schematic diagram is 
presented in Fig.2, was selected. For this cycle influence 
of limitations of parameters and selected changes 
in configuration were analyzed, and selected results 
of simulation calculations have been presented here. 


Investigated object 
In Fig. 2 schematic thermal diagram of the 200 MW 


power unit is presented. The arrangement of measurement 
instruments intended for measuring the heat-and-flow 
parameters such as: pressure, temeperature, mass flow 
rates and — additionally— also electric ouput power, is 
given. Majority of the parameters is dependent on load 
state and environmental conditions of the object. And, 
a part of them just indicates the load and environmental 
conditions, and to consider them as independent 
parameters, is possible. For the 200 MW condensation 
steam power unit the following parameters belong to the 
dependent ones: 


* electric generator output power or live steam mass 
flow rate, alternatively, 


* live steam pressure, 

* live steam temperature, 

* interstage superheat temperature, 
* condensation pressure, 


* pressure in the degasifier (since the constant 
-pressure degasifier is applied), 


erate of water jet to live steam superheater, 


erate of water jet to reheat steam superheater. 
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Fig. 4. Illustration of numerical calculation scheme of the power unit [5, 10] 


The measurement points ofthese quanities are indicated 
by irregular asterisks in the schematic diagram, Fig.2. 
Results ofthe measurements can be used to calculate the 
cycle efficiency characteristics. Measurement results of 
the dependent parameters together with the so determnied 
characteristics can be taken as current values in 
heat-and-flow diagnostics. They next serve for determining 
inefficient operation symptoms. 


The efficient operation standard can be determined 
computationally. For this purpose the steam power 
plant calculation model [5, 8, 10, 15] as well as the 
DIAGAR numerical calculation program [5] intended 
for balancing thermal cycles and performing their 
diagnostics, can be used. Fig.4 shows graphically the 
calculation model of the power plant presented in 
Fig. 2. Every square symbolizes a calculation procedure 
for particular component devices of the cycle, which 
makes use of geometrical and flow data of a given device, 
whereas the connecting lines reproduce structure of the 
power plant. With the use of the schematic diagram 
measurement points of heat-and-flow parameters can be 
identified. After introduction of independent parameters 
to the calculations, by hand or directly from DCS 
system, dependent values, including those associated 
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with measurement points, are calculated. Reference 
(standard) values are determined by taking into account 
undegraded geometry of component devices. 


The calculation program fulfils role of the efficient 
operation standard. It is a kind of functional standard. It 
makes it possible to calculate reference values in function 
of changes of independent parameters. Structure of the 
model is flexible and allows to introduce modifications. 
Partial calculation models of the devices are susceptible 
to tuning [8, 10]; they can be tuned on the basis of correct 
experimental results and are characterized by accurate, 
i.e. convergent with the experiment, calculation results 
of heat-and-flow quantities [10, 15]. The tuning of model 
consists in determination of corrections to experimental 
coefficients appearing in calculation methods. The 
corrections take into account all phenomena which 
cause differences between measurement results and 
relevant calculation results. Such corrections for two flow 
coefficients: u, and u,, and two velocity coefficients: @ 
and y, were prepared for turbine stages, as well as the heat 
transfer coefficient k — for heat exchangers [10, 15]. The 
calculation model tuned with their use satisfies the above 
mentioned requirements for efficient operation standards. 


Application of efficient operation standard in the 
case of external limitations imposed on power 
unit operational parameters 


In service conditions such situations -usually resulting 
from minor failures -happen that power generation 
objects operate at limited inlet parameters of working 
media because of cost of long -lasting standby periods. 
This has been observed in land electric power stations. 
In the cases known to this author live and reheat steam 
temperatures were decreased (sometimes even by 30°C) 
because of overruning repairs of boiler superheaters. In 
marine conditions similar phenemena can also occur if 
they happen during voyage which should be continued 
without generating excessive risk to safety of the ship and 
its crew. In each of the above mentioned events to get 
approval to do it from the side of manufacturers of crucial 
devices (first of all — boiler and turbine) is necessary. 

The above mentioned case of limitation of live and 
reheat steam temperatures is exemplified by a fragment 
of line of expansion within turbine, Fig. 5. The nominal 
temperature ty, was lowered to that t’)). This resulted in 
change of run ofthe expansion line from the nominal state 
(continuous purple line) to the transient state (dashed 
line). Only minor influence on operation of the valves 
and control stage can be observed. In both the states the 
efficiency characteristics of the non-regulated part of 
the turbine are only slightly different. However pressure 
distribution within the turbine is changed distinctly. 
Such state of the power unit can be reproduced in the 
calculation model of the DIAGAR program. Hence, an 
efficient operation standard can be elaborated on the basis 
of the calculation model. 


The tuning of the calculation model at limited 
inlet parameters is excessively time-consuming and 
unprofitable. Periods of such power plant operation 
are limited and the time-consuming process of the 


24 POLISH MARITIME RESEARCH, Special issue 2009/S1 


tuning cannot be then started. Therefore tuning-up and 
corrections for the nominal state of operation at different 
loads should be used. 


Fig. 5. Simplified schematic diagram of steam flow changes at inlets to HP 
and IP cylinders in the case of limited superheating temperatures 


The corrections of the above mentioned coefficients 
of turbine stages and the heat transfer coefficients k for 
heat exchangers can be performed either individually 
or globally. The global approach consists in introducing 
the correction in function of the independent parameters 
of the cycle, as expressed by the below given relations 


(1)+(5), [5, 10]: 


The individual corrections (6)+(10) are functions 
of the local parameters, i.e.: generalized losses — for 
the velocity coefficients, as well as blade geometrical 
parameters and Reynolds number — for the flow 
coefficients, and functions of the heating steam and feed 
water inlet parameters — for the heat transfer coefficients 
k for heat exchangers [10, 15]: 
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By making use of the global corrections and owing to 
compliance with experiment, to perform more accurate 
calculations in conditions of power plant normal load 
as compared with application of local corrections, is 
possible. However the global corrections have also 
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disadvantages resulting from their elaboration on the 
basis of the narrow changeability range of independent 
parameters, Fig. 6. 
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Fig. 6. Example of experimental ranges of live steam changes in conditions 
of normal turbine load to obtain global correction relations,[15] 


The narrow range results from requirements of 
producers of equipment, especially boilers and turbines, 
who recommend to maintain values of the parameters 
relatively constant. It concerns live steam pressure as well 
as temperatures of primary and secondary superheating. 
And, the corrections result from regression relations. If 
the range used for determining the regression relations is 
exceeded the obtained calculation results will be loaded 
by large errors reaching even a few dozen percent [5, 10]. 


The local corrections make it possible to perform 
more accurate calculations in conditions of large changes 
of independent parameters, which distinctly exceed the 
ranges used in tuning the model. It results from only 
minor exceedance of ranges of independent variables 
of correction functions, Fig. 7, even at large changes of 
independent parameters of power plant [5, 10]. In this 
case linear correction functions are especially favourable 
[15]. The expected value of calculation errors does not 
exceed 1% [15]. 
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Fig. 7. Example areas of application of individual correction relations for 
HP and IP cylinders of steam turbine: (a) — for velocity coefficients, Q; 


(b) — for flow coefficients, u1, [10, 15] 


Hence, the application of local correction makes it 
possible to calculate reference values and values of 
symptoms with a satsifactory accuracy in conditions 
of significant limitations of power plant independent 
parameters. As a matter of fact symptom determination 
error is greater than that in conditions of normal load but 
it makes it possible to use the so determined symptoms 
for diagnosing efficient operation of steam turbine cycles 
even in the so changed operational conditions. 


Application of efficient operation standard in 
the case of cycle configuration changes 

In operational practice cycle configuration changes 
can be met more often than the lowering of independent 
parameters. For land power plants for instance it concerns 
mode of feeding the degasifier, Fig. 8, or unsealing 
the valves in bypasses of heat exchangers, Fig. 9. 
Configuration changes often result from short current 
overhauls of component devices of power unit or also 
from producing low power. For marine power plants 
a change in utilization of steam for overall ship purposes 
orchange in electric output power ofturbogenerators may 
lead to changes in configuration of the cycle. 


HA ge cee Bee 
Fig. 8. A fragment of schematic diagram of power unit numerical 
calculations, which illustrates reproducing the cycle configuration change 
which consists in choice of various extractions (No. II marked 35/349, 
or No. III marked 75/349) for feeding the degasifier (marked 375/440), [5, 10] 


In the case of possible power plant operation at 
configuration changes of connections between component 
devices, calculation procedure of reference state 
parameters should take into account their occurence. 
The above mentioned numerical calculation program 
DIAGAR [5] fulfils such condition. Due to its modular 
structure which makes it possible to reproduce power 
unit's structure it can contain modules intended for the 


Fig. 9. A fragment of schematic diagram of power unit numerical 
calculations, which illustrates reproducing the cycle configuration change 
which consists in activation of the bypass of HP regeneration exchangers 

— yellow line, [5, 10] 
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controlling of working media flow to various devices of 
the power unit. The control modules of the program can 
be commenced manually by DIAGAR system operator, 
or automatically by reading valve opening indicators of 
DCS system. Values of dependent parameters necessary 
to determine inefficient operation symptoms are then 
calculated by such system adequately to an operational 
state of steam turbine power plant. 


Significance of reproduction of power object structure 
is illustrated by the numbers which determine selected 
parameters of powerplantandits efficiency indices, presented 
in Tab. 1. They are obtained by means of simulation 
calculations of the chosen power unit of 200 MW 
output. Characteristics of component devices comply 
with reference conditions and are not subjected to 
changes; the changes concern only inlet parameters or 
configuration. The first column deals with the operation 
in conditions of nominal configuration of the power unit, 
which consists a.o. in feeding the degasifier from the 
extraction III and closing all bypasses of heat exchangers. 
Three successive columns deal with changes in the cycle. 
The second one reproduces the operation of the cycle at 
lowered temperatures of live and reheat steam. The third 
one reproduces the operation of the cycle at feeding the 
degasifier from the extraction II. And the fourth one 
reproduces the operation at partial opening the bypass 
of HP regeneration exchangers. 


Tab. 1. Changes of parameters and characteristics of 200MW power unit in 
the case of introduction of configuration changes with respect to nominal one 


Changed inlet parameters or configuration 


Feeding the 
degasifier 


Nominal 
configuration 


Lowered 
temperatures: 

t and t 

00 M 


Descrition Symbol Units HP bypass 


from the active 
extraction II 


Power N MW 196.700 196.700 196.700 196.700 


Live steam 


M p, MPa 12.500 
pressure 00 


12.500 12.500 12.500 


Live steam t *c 539.000 500.000 539.000 539.000 
temperature 00 


Reheat steam o, 


539.000 
temperature M 


500.000 539.000 539.000 


Condensation 


p kPa 9.770 9.770 9.770 9.770 
pressure c 


specie heut q wd 9221.000 9309.000 9238.000 9266.000 
consumption j kWh 


HPcyl. 


RT n - 0.695 0.695 0.695 0.695 
efficiency 


IP cyl. 


efficiency n - 0.889 0.889 0.889 0.889 


LP cyl. 


ot n 0.801 0.801 0.801 0.801 
efficiency LP 


Feed water 
temp. afler the L ety c 
last HP heater "T 


219.000 219.000 219.000 210.000 


Terminal temp. 
differ. of the last r 
HP heater 


Cc 0.107 0.107 0.107 0.107 


Values of the parameters in Tab. 1 indicate that the 
only quantity susceptible to changes is the specific heat 
consumption qj. In a content-related way it is justified 
because only the parameters assumed nominal and only 
the nominal configuration ensure to reach maximum 
value of object’s efficiency (i.e. the lowest specific heat 
consumption). The remaining quantities and characteristics 
determined by calculation are not subjected to changes. 
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It shows that the use of the DIAGAR system calculation 
model as an efficient operation standard is correct even 
in changed external conditions and for operation of the 
cycle with changed configuration of connections between 
component devices. 


Summary 


The above presented analysis and selected results of 
simulation calculations indicate that there are applicable 
solutions which make it possible to solve certain problems 
concerning the power plant reference state of power 
objects, including marine ones. The problem is important 
for heat-and-flow diagnostics carried out with the use of 
object’s efficient operation standard. 


In this paper the problem of influence of limitation of 
values of working media inlet parameters and changes in 
configuration of connections between component devices 
of thermal cycles — with respect to conditions deemed 
nominal — on determining a reference state in the case 
of correct geometry of thermal cycle component devices, 
was analyzed. 


It was stated that the use of the modular calculation 
model as an efficient operation standard makes it possible 
to determine the reference values in the described 
changeable conditions. Such model tuned to a diagnosed 
power plant on the basis of correct measurements at 
nominal configuration will comply with accuracy criteria 
required for that standard also in the case of changed 
cycle configuration. 


In the case of limitation of working media inlet 
parameters, corrections to experimental coefficients of 
calculation methods should be then functions of local 
parameters in the surrounding of the device modeled by 
them. Then even large changes in power plant independent 
operational parameters make only slight impact on values 
of variables of the corrections. In simulation calculations 
no significant errors in determining characteristics of 
component devices, are observed, Tab. 1. 


The modular calculation model for reproducing 
a modified structure of power plant should contain special 
modules to control calculations, depending on current 
cycle configuration. Control parameters can be inserted 
by hand. However the use of information received from 
contemporary DCS systems makes it possible to automate 
the process by gaining data from valve opening indicators. 
In the case of correct geometry of component devices, 
modification of object’s structure does not change their 
correct characteristics, that is illustrated by the values 
given in Tab. 1. 


Comparison of the values for the nominal configuration, 
given in Tab. 1, with those for the changed configuration 
and independent parameters shows that the using of 
the modular calculation model as an efficient operation 
standard is correct even in the case of changed external 
conditions and the operation with changed configuration 
of connections between component devices. 
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Nomenclature 

b, BC - profile chord, 

i — enthalpy, 

k — heat transfer coefficient, 

l — blade length, 

p — pressure, 

Poo  — live steam pressure, 

Pe — condensation pressure, 

Pm — reheat steam pressure, 

qi — specific heat consumption, 

t — working medium temperature or 
cascade pitch, 

too — live steam temperature, 

tagp — feed water temperature after the last HP 
heater, 

tm — reheat steam temperature, 

X, — generalised primary profile losses, 

X, — generalised secondary profile losses, 

X, — steam dryness at the end of isentropic 
expansion, 

^ — deviation between actual and reference 
values, 

Tap — —terminal temperature difference of the 
last HP heater, 

9 — nozzle blade velocity coefficient, 

y — rotor blade velocity coefficient, 

u — blade flow coefficient, 

n — efficiency. 

Indices: 

1 — stands for nozzle blade cascade, 

2 — for rotor blade cascade, 

HP — for HP turbine cylinder, 

IP — for IP turbine cylinder, 

LP — for LP turbine cylinder, 

obl  — for calculated values, 

p — for steam, 

w — for water. 
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Design analysis of Tesla micro-turbine operating 
on a low-boiling medium 


Piotr Lampart, Krzysztof Kosowski, 
Marian Piwowarski, Lukasz Jędrzejewski 


Abstract 


This paper presents results of the design analysis of a Tesla bladeless turbine intended for a co-generating micro-power plant 
of heat capacity 20 kW, which operates in an organic Rankine cycle on a low-boiling medium. Numerical calculations of flow 
in several Tesla turbine models were performed for a range of design parameters. Results of investigations exhibit interesting 
features in the distribution of flow parameters within the turbine interdisk space. The calculated flow efficiency of the investigated 
Tesla turbine models show that the best obtained solutions can be competitive as compared with classical small bladed turbines. 


Keywords: Tesla turbine, axial flow bladed turbine, flow efficiency, CFD calculations, RANS model 


Introduction 


It is a rather common opinion that bladed turbines are 
unbeatable as compared with Tesla turbines in the range 
of large power outputs. However, at very small outputs 
some superiority of bladeless turbines can be expected, 
Fig. 1. In a majority of published cases, rotational 
speeds of Tesla turbine rotors appear lower than those of 
bladed turbines, which makes the selection of an electric 
generator easier. 


efficiency 


output 


Fig.1 Pictorial comparison of efficiency of small bladed and Tesla turbines 


Fig. 2 Tesla turbine [11] 
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The first bladeless turbine was designed and 
manufactured by Nikola Tesla (in 1913) [11]. The design 
makes use of the effects which occur in the boundary layer 
flow between the rotating disks placed very close to one 
another (Fig. 2). The gas flowing spirally from the outer 
to inner part transfers energy to the rotating disks. The 
supply usually takes place from several nozzles discretely 
located along the circumference. The medium flows out 
through the holes in the disks situated near the turbine 
shaft. 


Distances between disks are very small. The highest 
value of efficiency appears when they are approximately 
equal to the double boundary layer thickness [9]. 
The efficiency of the Tesla turbine depends on many 
parameters, namely on: pressure, temperature, inlet 
medium velocity, number, diameter, thickness and 
distance between the disks as well as on the state of the 
disk surface, rotational speed of the rotor, flow kinematics 
at the inlet to and outlet from the turbine, etc. In the 
subject-matter literature, examples of experimental 
research referring to the following models of Tesla micro- 
turbines can be found: 


* ~3.0kW output power, 15000 rpm, 32% efficiency, [5]; 
¢~1.0kW output power, 12000 rpm, 24% efficiency, [1]; 

* —].5kW output power, 12000 rpm, 23% efficiency, [8, 9]; 
*—50W output power, 1000 rpm, 21% efficiency, [7]; 

* —].5kW output power, 12000 rpm, 49% efficiency, [3]. 


Experimental works aimed first of all at establishing 
relationships between the turbine efficiency and 
parameters given below: 


* distance between the micro-turbine disks; 

* number and diameter of the micro-turbine disks; 
* number of inlet nozzles to the micro-turbine; 

* rotational speed of the rotor; 


* medium inlet pressure; 


* medium inlet temperature; 
* medium inlet velocity and inlet angle; 
* corrosion and erosion of micro-turbine elements; 


*constructional materials (composites, ceramic 
materials, bronzes, aluminium alloys); 


* kind of medium flowing through the micro-turbine 
(air, biogas, organic agents, exhaust gases, multi- 
phase media, etc). 


Main design features and operational parameters of 
Tesla turbines published in the literature referenced in this 
paper are presented in Tab. 1. Not only turbines but also 
compressors, pumps and gas turbine sets can be built on 
the basis of the same working principle. 


Preliminary design analysis 


CFD calculations of various models of Tesla disk 
turbines were carried out on the basis of the RANS 
model supplemented by the k-œ SST turbulence model 
available in the computer programme Fluent [12]. 


equations is performed using the finite volume method. 
The “segregated” solver with the sequential solving of 
the governing equations as well as the SIMPLE algorithm 
for correction of pressure and velocity are applied. 
Discretisation of convection fluxes is performed using 
an “upwind” scheme of the 2™4 order accuracy. The 
time-domain discretisation is made by an "implicit" 
scheme. The calculations are carried out until the 
stationary state 1s reached, lowering the residua of 
particular equations by 4+6 orders of magnitude. 


The calculation domain was prepared by means of 
the software Gambit [14]. The models of Tesla turbines 
contain the inlet nozzle, a half ofthe single interdisk space 
(under the assumption on flow symmetry in a central 
plane) as well as the outlet area. The models analyzed in 
this work represent a simplification of real geometry of 
Tesla turbines. The outlet area was simplified in a way to 
allow the medium outflow from the turbine flow passage 
along the entire circumference of the shaft. The presence 
ofthe tip clearance was neglected. The calculations were 
carried out in a fixed (motionless) reference frame, with 
turbine disk walls kept in rotational motion. In accordance 
with the assumed turbulence model, the calculation grid 


Numerical discretisation of the set of fundamental 
Tab. 1. Specification of Tesla turbine design solutions 
Pet Disk dnterdisks Pot Number | Rotational | Moment | Power | Medium | Mt Het Nozzle herd Nozzle | penciency 
isk thickness distance disk of disks speed pressure temperature diameter BEN angle 
No 
Di t b D5 Zz n T N PO t D S a n 
m] [m] [m [m] H1 [rpm] [Nm] [kW] [kPa] PC] [m] [mm] [] HH 
1 - - - 7 6 18000 2,983 Air 615,29 - - - 0,35 
2 - 8193 8,650 Air 593,00 - - 
3 0,17780 0,0024 0,0016 7 10000 2 Air 861,84 - - = 5 0,232 
4 | 0,17780 0,0024 0,001 11 9200 : Air 377,14 . : - 0 0,258 
5 0,20320 0,0005 0,0005 0,0335 24 17000 Air 615,29 - - : 20 0,35 
6 | 0,30480 0,0008 0,0008 | 0,07620 45 1100 0,447 Air 227,53 20,5 z 121 5 0,16 
7 | 0,30480 0,0008 0,0008 | 0,07620 45 6218 3,430 Gas 275,79 444 - - 5 0,123 
8 | 0,30480 0,0008 0,0008 | 0,07620 45 6284 3,206 Biomass | 275,79 391,6 - 5 0,11 
9 | 0,30480 0,0008 0,0008 | 0,07620 45 6500 9,247 | Saturated | 689,48 170 - 5 0,137 
vapour 
0 | 0,25240 0,0016 0,0016 9 6300 = Air s . - 206,67 5 0,41 
1 - 15000 2,983 - - - - 0,32 
2 - 12000 0,969 - - - 0,24 
3 : 12000 : 1,491 : B : 0,23 
4 z z 5 2 - 1000 - 0,045 s - 0,21 
5 - - a - : 120000 - 1,491 - 3 0,49 
6 | 0,30480 0,0008 0,0008 : 45 12000 " 4,4-50 Vapour/ | 689,48 537,78 0,0191 E : 0,25 
exhaust 
gas 
7 | 024765 0,0048 0,0016 . 7 : - 14,914 y - 7 7 E g 
8 | 0,40000 - - - 26 5400 - - 150 : : - 
9 | 0,25400 : S z - g - 74,569 - 2 T a : - 
20 | 0,45720 - : . - 3 z 149,138 g - : 3 : - 
21 | 1,52400 - - - E - - 503,341 - - T E 
22 | 0,24765 : - - 25 216 82,026 - " - z 
23 | 0,15240 0,0009 0,0012 - 7 1290 17,27 : - - - - 
24 | 0,15240 - - 3000 17,27 - - 517,11 - a : 
25 - - - 9240 2,908 Air 723,95 444 2 15 
26 - 10200 9,694 Air 723,95 3 - 15 
27 | 0,11430 24000 0,250 Vapour 10342 - - 
28 | 0,25400 z - z 29 18000 96,940 " " 2 
29 | 0,07620 f 3 15 10500 0,153 - 448,16 ` - 0,31 
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at the disk wall was refined so as to obtain the y, value 
equal to 1 + 2. This is a structural mesh divided into 
blocks, which contains 400 000 + 500 000 finite volumes, 
depending on a turbine model. The calculation grid was 
also refined in the inlet and outlet regions. Fig. 3 presents 
images of the calculation area for several turbine models 
with one, two and four inlet nozzles. Geometrical parameters 
of particular investigated models are given in Tab. 2. 


Thermodynamic parameters assumed for CFD 
calculations were found from preliminary 1D model 
calculations, making use of data from literature sources. 
It is assumed that the considered Tesla turbine models 
consist of 11 rotating disks (12 flow channels of the 
interdisk space). The outer diameter of the disk is equal 
to 10 cm. The nominal operating conditions are for 
the mass flow rate of 0.13 kg/s and pressure drop from 
14.8 bar to 1.9 bar. Solkatherm®SES36S was assumed 
to be a working medium. The perfect gas model was 
chosen for the calculations, assuming the individual 
gas constant and specific heat as average values from 
the given expansion range. Values of dynamic viscosity 
and heat conductivity coefficient were assumed in 
a similar way. Pressure boundary conditions relevant to 
compressible flow solution were assumed. However, it 
should be noted that the perfect gas model assumed for 
the calculations may be a poor approximation of working 
medium properties, especially for flow velocities close 
to the sonic velocity. The calculations were carried out 
for a range of operating conditions (by changing the 
available pressure drop) for two rotational speeds of the 
rotor, namely: 18 000 rpm and 9 000 rpm. 


Contours of static pressure and velocity for the 
investigated turbine models 1-3 are presented in Figs. 4, 5 


Fig. 3 The calculation domain for four models of the Tesla turbine 


Tab. 2 Geometrical parameters of Tesla turbine models 


Model 1 Model 2 | Model3 | Model 4 


| oo | 0.0005 0.00025 0.0005 
0.003 0.003 0.003 0.008 
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Contours of absolute Pressure (pascal) 
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Contours of velocity Magnitude (m/s) 
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Fig. 4. Contours of static pressure and velocity in flow through the four- 
nozzle model of the Tesla turbine; 6=0.25mm, n=18 000, pin= 14,8 bar, 
T,,=410K, G = 0,132 kg/s, P = 1177 W. 
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Fig. 5. Contours of static pressure and velocity in flow through the two- 
nozzle model of the Tesla turbine; ó—0.5mm, n=18 000, pin— 14,8 bar, 
T,,=410K, G = 0,124 kg/s, P = 1102 W. 
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and 6. The numbers below the figures show the interdisk 
gap, rotational speed, inlet pressure and temperature (at 
the outlet pressure of 1.9 bar) as well as the resultant 
values — rate of mass flow through 12 interdisk passages 
and output power generated on both sides of 11 rotating 
disks. It results from the pressure contours that a part of 
the available pressure drop is accomplished in the nozzle, 
the other part takes place in the interdisk space, i.e. in the 
rotor. In the case of higher loads, configurations of isolines 
characteristic for the occurrence of shock waves are 
observed some distance downstream of the inlet nozzle. 
When a shock wave occurs, an increase of pressure and 
a decrease of flow velocity takes place. Streamline 
patterns for the four-nozzle and two-nozzle models 
operating under part-load conditions (inlet pressure at 7.8 
bar) are presented in Fig. 7. It is clear that fluid elements 
make more than three or four rotations within the interdisk 
space before they reach the outlet section. 


A sample distribution of static pressure and velocity 
in flow through the model 4 of Tesla turbine (of a small 
shaft diameter) is presented in Fig.8. Such a system 
is characterized by a high reaction, which means that 
practically the whole pressure drop takes place in the 
interdisk space of the rotor and that the velocity at the 
nozzle exit is relatively low. This is due to the vortex 
motion induced in the interdisk space, which results 
from the conservation of angular momentum. Thus, 
flow transition to smaller diameters towards the outlet 
section causes a velocity increase and a pressure drop. 
This situation is unfavourable as far as the efficiency is 
considered — high outlet energy means high flow losses. 
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Fig. 6 Contours of static pressure and velocity as well as form of pathlines 
in flow through Tesla turbine model; d=Imm, n=18 000, pin= 10,8 bar, 
T,,=390K, G'— 0,024 kg/s, P’= 109,7 W 


Fig. 7. Streamlines in flow through the four-nozzle and two-nozzle models of 
the Tesla turbine; n=18 000, pin= 7,8 bar, T;,=390K. 


Therefore, the outlet section of a disk turbine should be 
located at an appropriately large diameter. 


The power output, inlet pressure required, flow 
efficiency and reaction of three investigated models of 
11-disk (12-passage) turbine are presented in Figs. 9, 10 
and 11 as a function of mass flow rate. The model with 
four nozzles gives the largest power and flow efficiency. 
Values of power and efficiency obtained from the two- 
nozzle model are slightly lower. The respective values 
for the single-nozzle model are considerably lower. The 
change in rotational speed from 9000 rpm to 18000 rpm 
leads to a significant power and efficiency increase. For 
the rotational speed of n=18000 rpm, at the nominal 
pressure drop from 14.8 to 1.9 bar, which refers to the 
nominal flow rate of about 0.13 kg/s, the output power 
of the investigated models of 11-disk turbine amounts to 
about 1020 W for a single-nozzle model, increasing to 
1100 W for the two-nozzle model, to 1180 W for the four- 
nozzle model. At n=18000 rpm, the calculated efficiency 
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Fig. 8. Contours of static pressure and velocity in flow through the model 4 
of the Tesla turbine with a small shaft small diameter; 6=0.5mm, n=18 000, 
pin= 14,8 bar, T,,=410K. 
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Fig. 9. The power output, inlet pressure, flow efficiency and reaction of an 
1I-disk turbine as a function of flow rate; the system supplied from four 
nozzles -ó—0.25mm. 
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Fig. 10. The power output, inlet pressure, flow efficiency and reaction of 
an 11-disk turbine as a function of flow rate; the system supplied from two 
nozzles -ó-0.5mm. 


ofthe models at the nominal operating conditions is equal 
to 25.5% for the single-nozzle model, increasing to 27% 
for the two-nozzle model, to 2896 for the four-nozzle 
model. For the rotational speed of n=9000 rpm, at the 
nominal flow rate of 0.13 kg/s (and nominal pressure drop 
from 14.8 to 1.9 bar) the output power of the investigated 
models of 11-disk turbine amounts to about 700 W for 
a single-nozzle model, increasing to 760 W for the two- 
nozzle model, to 780 W for the four-nozzle model. At 
n=9000 rpm, the calculated efficiency of the models at 
the nominal operating conditions is equal to 1796 for the 
single-nozzle model, increasing to 18.2% for the two- 
nozzle model, to 18.5% for the four-nozzle model. It is 
also important to note that for lower loads below nominal 
operating conditions the flow efficiency is increased. The 
reaction of the models at n=18000 rpm varies within the 
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Fig. 11. The power output, inlet pressure, flow efficiency and reaction of 
an 11-disk turbine as a function of flow rate; the system supplied from one 
nozzle -ó- Imm. 


range of 0.4 — 0.5, at n-9000 rpm changing between 
0.25 + 0.35. The characteristics of the model no. 4 are 
here omitted because of a low value of flow efficiency. 
The internal efficiency of this model is not greater than 
10% over the whole range of investigated loads. 


Closing remarks and direction of future works 
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Fig. 12. Contours of static pressure and velocity in flow through the large 
diameter two-nozzle model of the Tesla turbine; d=32mm, 6=0.5mm, 
n=9 000, pin= 7,8 bar, Tin=390K, G = 0,119 kg/s, P = 1312 W 


The flow efficiency of the investigated models 1-3 
of a Tesla turbine obtained by means of numerical 
calculations from the code Fluent fall within the range 
of efficiency values quoted in Introduction based on the 
literature search. The calculations also highlight some 
features of operation of the Tesla turbine from a point 
of view different than that presented in a majority of 
publications. First of all, the Tesla turbine operates as 
a reaction turbine with a large pressure drop within the 
interdisk space. The diameter of the outflow area from 
the intersdisk space should be still large enough to avoid 
large outflow velocities. It can also be expected that the 
lengthening of the friction path within the interdisk space 
would result in a more effective power transfer. 


The above conclusions lead to an idea of increased 
disk diameter of the Tesla turbine, also accompanied by 
a decreased cross-section area of the inlet passages. To 
verify this idea, preliminary calculations were performed 
for two variants of disk diameter of the two-nozzle model 
of the Tesla turbine scaled up to 16 cm and 32 cm. In 
both cases the distance between the disks was equal to 
0,5 mm. The supply of medium was accomplished by 
means of two nozzles located opposite to each other. 
Sample distributions of pressure and velocity fields in the 
interdisk space of the Tesla turbine model with the outer 
diameter equal to 32 cm are shown in Fig. 12, whereas 
power and efficiency characteristics are displayed in 
Fig. 13. Within the presented range of mass flow rate 
(between 0.1 and 0.4 kg/s) and pressure drop (between 
7.8 and 14.8 bar) the calculated flow efficiencies oscillate 
around 50%. 


It seems justified to conclude that precise optimization 
of the Tesla turbine design and geometrical parameters 
could ensure values of the output power and efficiency 
similar to those achieved by small bladed turbines. 
However, the rotational speed of the Tesla turbine rotor 
is several times lower as compared with small bladed 
turbines of rotational speed 100 000 rpm. This makes 
the Tesla turbine an attractive proposal for small heat and 
power stations, e.g. of heat capacity 20 kW. 
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Fig. 13 The power output and efficiency of an 11-disk turbine of a large 
shaft diameter (32 cm) as a function of flow rate; the system supplied from 
two nozzles -ó—0.5mm. 
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Abstract 


This paper presents results of a design analysis of turbines for co-generating micro-power plant working in accordance with 
organic Rankine 5 cycle and using biofuel. The heat power range from 25 kW to 100 kW with corresponding available electric 
power from 2kW to 12kW, was considered. Designs of axial-flow turbines (single-stage and multi-stage ones, also those partially 
fed), radial-flow and axial-radial -flow ones, were analyzed. Particular variants of the solutions were compared to each other. 


Keywords: micro power plant, microturbines, organic Rankine cycle, turbine design 


Introduction 


In the subject-matter literature steam micro-turbines 
are considered to be such devices which produce output 
power of the order of a few kW or even W. For last years 
greater and greater interest paid to such machines and 
their dynamic development has been observed. Hybrid 
co-generating systems, either three-generating or 
combined gas-steam ones can be most often met 
[1-3,10-13,16,18-21,24,25,27]. Steam turbines may find 
application to micro-power plants which co-generate 
electric power and heat. A micro-power plant based on 
water steam turbine, described in [14], can serve as an 
example. The system delivers electric power amounting 
from 0,5k W to 4,6 kW as well as heat power - from 
2 kW to 25 kW. In the literature can be met examples 
of steam micro-turbines of much lower values of output 
power, built with the use of the MEMS technology 
(Micro-Electro Mechanical System) [6+8]. Apart 
from the traditional medium, i.e.water vapour, also 
low-boiling media, as a rule organic ones, are taken into 
consideration. In this case the Organic Rankine Cycle 
(ORC) is dealt with. For many years such installations 
of electric power output of the order of several hundred 
kW or MW have been used in power plants based on 
geothermal sources [4,5,17,26]. On the market are 
already available the co-generating systems working 
on organic media, e.g. the power plant of 300+600 kW 
electric power and 1500+2800 kW heat power, or that 
of 200+1000 kW electric power and 1000-6000 kW 
heat power [9]. However only a few examples can be 
found of ORC installations of output power smaller than 
100 kW, which operate in co-generating systems (e.g. the 
ORC-CHP system intended for biomass combustion [23], 
or the systems generating: 72 kW [22] or 30 kW [15] 
of electric power). Under research are: a turbine set of 
5-12 kW electric power at 70+115°C inlet temperature 
of working medium (R134a, R245fa, R22, R6xx, 
R7xx) [15], and a power plant working on n-pentane 
and developing 1.5+3 kW power (at 90+100°C vapour 
temperature) [21 ]. 
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Problems faced by designers of such turbines are 
associated with very small volumetric flow rate of 
working medium. It leads to small values of height of 
flow-part blades (hence to increased loss of flow rate) and 
high values of rotor speed, namely from a few dozen to 
over a hundred thousand rpm, usually. 


Vapour 


Generator 


High-temperature 
heat exchanger 
i Low-temperature: 
heat exchanger 


Feed pump 


Fig. 1. Schematic diagram of micro-power plant cycle 

In all the considered variants the power plants operated 
on the Solkatherm®SES36 working medium. Calulations 
were performed in compliance with the power plant 
schematic diagram shown in Fig. 1, for four values of 
the cycle’s thermal power: 25 kW, 50 kW, 75 kW and 
100 kW. Particular variants, depending on a heat power 
value, differred from each other only by value of working 
medium flow rate, and its parameters at turbine inlet, as 
well as pressure behind the turbine were the same for all 
the considered cases: 


-the pressure at turbine inlet, Po = 1464 kPa, 
-the temperature at turbine inlet, tọ = 140 °C, 
-the pressure behind the turbine, p, = 189 kPa. 


Within the frame of performed design analyses the 
following types of microturbines were taken into account: 


e single-stage axial-flow turbine, 
e two-stage axial-flow turbine, 

* four-stage axial-flow turbine, 

e single-stage radial-flow turbine, 


e single-stage radial-axial-flow turbine. 


Fig. 2. Examples of the considered constructional types of micro-turbines: 


a) single-stage axial-flow turbine, b) two-stage axial-flow turbine, 


c) radial-axial-flow turbine, d) radial-flow turbine 


m 
V 


B» 


Fig. 3. Examples of flow part of micro-turbines: a) single-stage axial-flow 


turbine, b) two-stage axial-flow turbine 
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assumed values of the reactivityp, a) plants heat power — 100 kW, b) plants heat power — 75 kW 


POLISH MARITIME RESEARCH, Special issue 2009/51 


Fig. 4. Examples of design characteristics of single-stage axial -flow turbine: turbine output power in function of velocity ratio for three assumed values of 
the angle a, turbine output power in function of velocity ratio for three assumed values of the reactivity p, rotor speed in function of velocity ratio for three 
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Fig. 5. Examples of design characteristics of single-stage axial-flow turbine: turbine output power in function of velocity ratio for three assumed values of 
the reactivity p, rotor speed in function of velocity ratio for three assumed values of the reactivity p, 


a) plants heat power — 50 kW, 


Examples ofthe considered design solutions are presented 
in Fig. 2. For all the variants preliminary optimization of 
main design parameters was performed and appropriate 
blade profiles were selected. Examples ofthe flow part of 
single-stage and two-stage axial-flow turbines are given 
in Fig. 3. 


Results of design analyses of micro-turbines 


Calculations of flow-part of turbine were performed 
by assuming different sets of values of main design 
parameters such as: velocity ratio, reactivity and 
outlet angle of guide vanes grid. Examples of design 
characteristics are presented in Fig. 4 for plant's heat 
power values: 100 KW and 75 kW, and in Fig. 5 for plant's 
heat power values: 50kW and 25kW. For the greatest heat 
power value (100 kW) of the considered co-generating 
systems, the electric power value of ~11,5 kW can be 
obtained at the rotor speed of about 45000 rpm. For the 
smallest considered heat power value (25 kW) the above 
mentioned quantities amount to about 3.35 kW and 80000 
rpm, respectively. 


Possible elaboration of series of types of 
turbines which would be suitable for the power 
range of 25 kW + 100 kW and of the same design 
and similar gabarites, has been also considered. 

It was decided that all turbines, regardless of their 
power, would be built of the same rotor and the only 
difference would consist in adjusting the feeding arc of 
vane blades grid, respectively. It was assumed that the 
100 kW turbine would operate with full feeding arc, hence 
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b) plant 5 heat power — 25 kW 


turbines of smaller power values would operate with 
smaller feeding arcs, respectively. Examples of design 
characteristics of single-stage axial-flow turbines working 
with partial feeding are presented in Fig. 6. Application 
of feeding arc makes that power values of turbines 
working with partial feeding are smaller than those of 
turbines with full feeding (compare Fig.4 and Fig.5 with 
Fig.6), however in the variants of smaller power values 
rotor speed values were significantly lower, e.g. for 
the plant's heat power of 25 kW, when partial feeding 
has been applied, the turbine's electric power decreased 
from 3,35 kW to 2.65 kW, but with accompanying drop 
of rotor speed from about 80000 rpm to about 40000 rpm. 
Significantly lower production costs would be reached 
as a result of standardization, application of rotors of 
the same design to all variants, and possible selection of 
a generator intended for lower rotational speeds (definite 
influence on price of turbine set). 


Similar design calculations were carried out for the 
remaining constructional types of turbines in question; 
their results are presented in Tab. 1. For the presented 
solutions suitable electric generators were selected; their 
main particulars are attached also in Tab. 1. 


For the variant of 100 kW heat power of the power 
plant the largest electric power values (— 12,3 kW) were 
obtained in the case of application of the radial turbine 
or two-stage axial-flow turbine. In most considered cases 
the radial-flow turbine as well as two-stage axial-flow 
turbine turned out to be the most profitable as regards the 
obtained electric power value. From this point of view the 
radial-axial-flow turbine appeared the worst. 


a) 


8000 80000 
70000 
7500 
60000 
7000 T 
= & 50000 
& — p -0.05 E +> p -0.05 
3 6500 -2-p-0.10) & 40000 -2-p-0.10 
a =p -0.15 5 =0.15 
- 2 30000 =p 0.15 
6000 4 [2 
20000 
5500 
10000 
5000 0 1 
02 025 03 035 04 045 05 055 06 02 025 03 035 04 045 05 055 06 
Velocity ratio [-] Velocity ratio [-] 
b) 
5400 70000 
5200 60000 
5000 
.., 50000 
4800 & 
E "o 40000 
= 4600 3 > p =0.05 
a 
8 4400 ? 30000 penu 
La $ —p-0.15 
4200 [2 
20000 
4000 
5800 10000 
3600 0 
02 025 03 035 04 045 05 055 06 02 025 03 035 04 045 05 055 06 
Velocity ratio [-] Velocity ratio [-] 
c) 
2700 70000 
2600 
60000 
2500 
2400 = 50000 
& 
2300 = 
E -*-p -0.05 mg 40009 — p =0.05 
$ 2200 i 
3 =p -0.10 2 = p =0.10 
$ x 30000 
& 2100 p =0.15| 8 = p -0.15 
o 
[4 
2000 20000 
1900 
10000 
1800 
1700 0 


0,2 0,25 0,3 0,35 0,4 0,45 0,5 0,55 0,6 
Velocity ratio [-] 


0,2 0,25 0,3 0,35 0,4 0,45 0,5 0,55 0,6 
Velocity ratio [-] 


Fig. 6. Examples of design characteristics of single-stage axial-flow turbine with partial feeding: turbine output power in function of velocity ratio 
for three assumed values of the reactivity p, rotor speed in function of velocity ratio for three assumed values of the reactivity p 
plants heat power — 75 kW, feeding arc of ~ 0,75 
plants heat power — 50 kW, feeding arc of ~ 0,50 
plants heat power — 25 kW, feeding arc of ~ 0,25 


Summary 


Calculations of a few dozen of variants of turbines 
of various constructional types and different values of 
main design parameters, were performed. The following 
turbines were taken into consideration: single-stage 
axial-flow turbines, two-stage axial-flow turbines, 
radial-flow turbines and radial-axial-flow turbines. It 
was demonstrated that in the case of co-generating 
micro-power plant working with low-boiling medium (in 
accordance with the ORC principle) the following electric 
power values are possible to be obtained depending on 
a heat power value of the system: 


-10,4 kW + 12,8 kW (for 100 kW heat power), 
-7,6 kW =+ 9,4 kW (for 75 kW heat power), 
-5,0 kW =+ 6,1 kW (for 50 kW heat power), 
-2,4 kW = 3,3 kW (for 25 kW heat power). 


Rotor speed values of the considered micro-turbines 
were relatively high as they amounted to from 21000 
rpm to 134 000 rpm depending on working medium mass 
flow rate, constructional type of turbine and choice of its 
main design parameters. For many variants one managed 
to select suitable electric generators. However, because 
of high rotational speeds, novel innovative design, 
material and technological solutions should be applied 
in developing power plants of the kind. 
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Abstract 


This paper presents an evaluation of technical efficiency and time-dependent changes in productivity of Polish professional 
electric power and heat & electric power stations by using Data Envelopment Analysis non-parametric method. The research 
covered the enterprises whose total available electric power amounted to about 98% of that of all professional thermal electric 
power stations in this country. The analysis has concerned the years 2000-2004. Impact was considered of such factors as: 
scale of enterprise activity, form of ownership as well as type of carried out activity, on technical efficiency and time-dependent 
productivity changes. Changes in productivity of the power generation sector were determined by means of Malmquist index. 

Keywords: Data Envelopment Analysis (DEA), technical efficiency, Malmquist index, thermal power generation 


Introduction 


An important factor which determines enterprise 
competitiveness level is optimum use of its resources to 
generate higher profits. In other words, the obtaining of 
higher economic efficiency by an enterprise impacts its 
competitiveness level. One of the economic efficiency 
components is technical efficiency which is understood 
as the achieving of the highest output in presence of 
a determined number of production factors. Level of the 
technical efficiency is influenced by a.o. technical progress 
and economic scale. 


One of the methods for measuring the technical 
efficiency is the Data Envelopment Analysis (DEA) based 
on the linear programming where appropriate technical 
efficiency levels are obtained by determining optimum 
values of objective function at certain constraints imposed 
on it. 


The sector of heat and electric power generation 
industry, crucial for Polish economy, was selected for this 
analysis. More comprehensive research into technical 
efficiency of power production eneterprises has been 
lacking so far. In this paper an evaluation of relative 
technical efficiency (as well as its components) of domestic 
professional electric power and heat & electric power 
stations is performed by using the DEA method, and of 
time-dependent changes in their productivity in the years 
2000+2004 - by applying the Malmquist productivity index. 


For power industry enterprises effective use of their 
resources, i.e. increasing technical efficiency level, is 
of special importance. In the case of the industry sector 
technological changes and rationalization of internal 
organization are ofa great influence on technical efficiency 
level. 


Technical efficiency and productivity 


Technical efficiency determines enterprise capability 
of achieving maximum level of effects at given inputs or 
consuming minimum amount of inputs to achieve assumed 
effects, at a given technology, 1.e. structure of production 
inputs. 


Economists are unanimous that the technical efficiency 
is the most important efficiency component as its level is 
influenced mainly by technical progress and better use of 
applied technology. 


The reaching of maximum technical efficiency means 
that from a given combination of resources an enterprise 
manufactures the maximum possible amount of goods and 
services, in other words it cannot increase amount of its 
effects without increasing amount of its inputs or reduce its 
inputs without reducing its effects, which is unambiguous 
to the statement that resources are not wasted during 
production process. And, the technical efficiency measure 
orientated towards inputs shows to which degree to reduce 
consumption of inputs without reducing amount of effects, 
is possible. However the technical efficiency measure 
orientated towards effects shows by which amount the 
generated effects can be increased without engaging any 
greater amount of inputs. The technical efficiency can be 
assessed independing on prices. 


A way of measuring the eneterprise technical efficiency 
consists in relating an observed level of production at 
given inputs to their potential level. And, the potential 
level is determined by limit values, i.e. minimum amount 
of component inputs (or their combimation) necessary to 
produce a given amount of effects (or their combination). 
Distance from the analyzed enterprise production level to 
the curve which shows relevant limit values is a measure 
of technical inefficiency. If the point which symbolizes 
current production of a given enterprise lays on the limit 
curve then the enterprise is assumed fully efficient; and if the 
point lays below the curve then the enterprise is considered 
technically inefficient. The ratio of observed output and 
potential output at given inputs and definite technological 
conditions constitutes the efficiency level of the enterprise 
in question. 


The components of total technical efficiency are as 
follows: 


TE = PTE x SE (1) 
where: 
TE - technical efficiency; 
PTE -purely technical efficiency; 
SE - scale efficiency; 


TE, PTE and SE - take values from the interval <0+1>. 
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The purely technical efficiency is the part of technical 
efficiency, which cannot be considered as that caused by 
deviations from the most favourable scale of enterprise 
activity. 


The enterprise productivity is the ratio of quantity or 
value of one or several categories of the effects generated 
to a certain quantity or value of one or several caregories 
of the inputs by means of which the effects have been 
obtained. 


Description of Data Envelopment Analysis 
(DEA) method 
The DEA method is aimed at measuring efficiency 
of a certain number of enterprises. Entities belonging to 
a considered group are compared to each other. The DEA 
method finds the best entities and shows to which degree 
others are inefficient as compared with those fully efficient. 


The investigations presented in [2], [4] have formed 
a background to relative efficiency evaluation.Their 
authors have introduced the Data Envelopment 
Analysis method which consitutes operational research 
methodology based on linear programming technique and 
aimed at measuring relative efficiency of the decision 
making unit (DMU) when it is difficult to compare 
general efficiency measures in situation of many kinds 
of component inputs and effects. 


The DEA method belongs to non-parametric ones. 
Hence, in contrast to the parametric methods, it does not 
necessitates a functional relation between inputs and effects 
to be determined, neglects inflation influence on results, 
does not make taking into account prices of production 
factors in technical efficiency calculations, necessary. 

The Efficiency Measurement System (EMS), 
version 1.3, a special computer software (elaborated by 
Universität Dortmund, Germany, and kindly released for 
scientific use) was used to perform the technical efficiency 
analysis for the industrial sector in question. This paper 
presents the technical efficiency analysis based on DEA 
measures orientated towards inputs, which show how 
far it can be reduced without worsening the effects. 

The DEA model orientated towards inputs, written as 
a linear programming problem, is of the following 
form: 


min fox". y“ Jer} Q) 


In the model, T stands for the set of production 
capabilities, defined on the basis of available empirical 
data about n objects; X^,Y^ -vector of k-th entity 
inputs and effects, respectively. The parametr © 
determines technical efficiency level of a given object 
and simultaneously informs which perentage of inputs 
would suffice to be kept in the enterprise to obtain present 
effects in the case of using the technology implemented 
in the enterprises deemed efficient. Solving the above 
mention problem consists in finding minimum value of 
© . The value of © for the most efficient objects is equal 
to 1, for the remaining ones the value is adequate to the 
efficiency level achieved by them and is smaller than 1. 
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In calculating the efficiency a kind of economic scale 
can be assumed, namely it used to be assumed that either 
constant scale effects or variable ones occur. Depending 
on that, different results as to different kinds of efficiency, 
can be obtained. Let's attribute the following symbols 
to the different efficiency measures depending on an 
assumed kind of economic scale: 


*M CRS - stands for the efficiency measure in the 


case of constant scale effects; 


*M VRS - stands for the efficiency measure in the 
case of variable scale effects. 


If the measures differ from each other then one will 
be able to conclude as to occurrence of scale effects in 
a considered sector. It can be assessed by means of the 
following measure: 


E S VRS- M CRS/M VRS (3) 


If it yields a value smaller than 1, then the entity in 
question can be said inefficient with respect to scale. 


In this paper comparisons of time-dependent 
productivity were performed by using the Malmquist 
index. The index applicable to comparing changes in 
efficiency has the following form: 


t[ tl tH tH (t+ tH 1/2 
M; (x.y, 2. y")- oan | PIS yj ) (4) 


where: 


D; (z^ E ) -the efficiency in the case of making use 
of t-th year technology for data from 
(t+1)-th year, 


-the entity’s efficiency in the period t 
in the case of making use of the then 
available technology and data from 
the period t. 


D(x, y") 


D (x, y") -the entity's efficiency in the period 
(t1). 


D" (x' , y) -the efficiency in the case of making 
use of (t+1)-th year technology 
for data from t-th year. 


And, the notion of technology is deemed to represent 
a level of technological development ofa given enterprise 
in a given period. Such approach to the technology makes 
it possible a.o. to investigate productivity improvement 
resulting from technological development. 


When Me yx, y?) is equal to or smaller 
than 1, then productivity increases, remains 
constant, or decreases, respectively, in the period 
(t+1), as compared with that in the period t. 


Data and results of the investigations 


In this work 39 eneterprises were analyzed for the 
years 2000+2001 and 38 for the years 2002+2004. And, 
27 heat & electric power stations and their complexes 
were analyzed for the years 2000+2001 and 26 of them 
for the years 2002+2004. Also 9 electric power stations 
and their complexes were analyzed in this work. Three 
complexes combining electric power stations and 
heat & electric power stations were taken into consideration 
as well. 


The model’s variables were selected on the basis 
of the literature sources dealing with DEA method 
application to analyzing technical efficiency in the power 
industry sector. The model was so formed as to comply 
with specificity of Polish electric power industry. 


In this work for every investigated entity six categories 
of input data to DEA model were assigned, four of which 
concerned inputs and two — effects. The input categories 
taken into account were: 


* mean yearly number of employees [number of 
persons]; 


e available electric output power [MW]; 
e available heat output power [MW]; 


* amount of chemical energy loading contained 
in fuel, e.g. coal, consumed to produce heat and 
electric power [TJ]. 


Amount of sold heat energy, counted in [TJ], and 
amount of sold electric energy [TJ] appear in the model 
as the effects. 


In this work efficiency measures for the successive 
years: 2000, 2001, 2002, 2003 and 2004, were calculated. 
Asthese authors have been obliged to respect confidentiality 
ofthe input data achieved from the enterprises the obtained 
results were assigned only to symbols of particular 
enterprises. The example results of estimation of the 
efficiency measures are given in Tab. 1. 


The technical efficiency measures orientated towards 
inputs inform on which perentage of inputs would be 
sufficient to produce such amounts of heat and electric 
energy as to operate as efficiently as the electric power 
stations which obtained the highest marks. For instance 
to get 70% measure value in a given year means 
that to produce the same amount of energy by using 
70% of inputs for this aim will be possible if only the 
enterprise in question uses its resources as efficiently 
as the best enterprises among the analyzed ones. 
The productivity analysis is more suitable to present time- 
dependent changes in results of the sector's operation 
than the presenting of mean measures of relative 
efficiency as the latter quantity may increase whereas 
the productivity may decrease. Apparent improvement 
expressed in the form of mean relative efficiency tells 
mainly about degree of differentiation of enterprises with 
a view of their skill in using inputs to transform it into 
effects. Whereas changes in productivity shows whether 


relation of effects to inputs has really increased or not. 
Diag.1 shows the mean measure of Malmquist index, 
taken from the measures as well as their components 
for all the investigated entities in the years 2000-2004. 


Tab. 1. Technical efficiency measures orientated towards inputs obtained by 
means of DEA model for the years 2000-2004, at the scale effects assumed 
constant (M CRS) 


Symbol of 
entity 


2000 2001 2002 2003 2004 


100,00% 
95,39% 
73,8196 
70.7796 
00,0096 
72.9696 
00,0096 
00,0096 
92,5096 
87,1896 
00,0096 
00,0096 
00,0096 
00,0096 
00,0096 
00,0096 
90,6496 
86,6496 
90,3396 
69,1096 
94,8596 
100,00% 
63,82% 
100,00% 
100,00% 
97,95% 
88,17% 
83,30% 
93,09% 
100,00% 
100,00% 
81,1396 
70,1796 
100,00% 
80,1196 
68,5796 
100,00% 
74,04% 


00,00% 
87,60% 
84,66% 
81,11% 
00,00% 
77,50% 
00,00% 
00,00% 
90,21% 
88,39% 
00,00% 
95,64% 
00,00% 
80,82% 
100,00% 
100,00% 
92,2196 
87,87% 
91,34% 
67,36% 
93,24% 
100,00% 
65,34% 
100,00% 
97,15% 
81,62% 
97,93% 
100,00% 
73,70% 
100,00% 
91,42% 
65,80% 
81,55% 
59,47% 
100,00% 
79,50% 
74,3296 
71,1396 
83.7196 


96,7996 
92,7696 
74,2396 
71,7496 
100,00% 
73,8496 
100,00% 
100,00% 
90,8796 
81,2896 
100,00% 
97,01% 
100,00% 
100,00% 
100,00% 
100,00% 
88,5896 
89.2796 
90,5296 
78,1096 
96,2196 
99,8996 
60,9996 
00,0096 
00,0096 
93,2996 
86,5696 
85,2996 
82,4896 
00,0096 
00,0096 
73,7496 
69,7596 
00,0096 
782196 
57,9596 
00,0096 
73,4896 


FI 00,0096 
91,6296 
95,5096 
83,7796 
00,0096 
82,6696 
00,0096 
00,0096 
98,2596 
91,8096 
00,0096 
83,3296 
00,0096 
86,06% 
00,00% 
00,0096 
00,0096 
85,4996 
89,6096 
65,5896 
96,2896 
100,00% 
62,45% 
100,00% 
96,9596 
81,7496 
98,8596 
100,00% 
70,7896 
100,00% 
84,0296 
50,6996 
98,8696 
57,3096 
100,00% 
74,8796 
72,7996 
76.1396 
83,2896 


100,00% 
92,86% 
79,88% 
76,28% 
100,00% 
76,02% 
100,00% 
100,00% 
94,719% 
90,46% 
100,00% 
100,00% 
97,64% 
87,03% 
100,00% 
100,00% 
91,82% 
84,97% 
93,13% 
71,89% 
89,78% 
99,98% 
67,58% 
100,00% 
100,00% 
83,47% 
89,32% 
89,37% 
76,4496 
95,4796 
100,00% 
65,5296 
54,96% 
100,00% 
716,36% 
12,17% 
76,9896 
81,5396 


The performed analysis of Malmquist index shows 
that the high productivity level for the year 2001 has 
resulted from the large amount of heat energy sold in this 
period, that can be explained a.o. by different atmospheric 
conditions, namely low temperature values. And, for the 
years 2002-2003 the low level of the index (abt. 0,696) 
resulting from economy slowing-down, can be observed. 
In the year 2004 the productivity increase is greater than 
in the preceding years. And, changes in relative technical 
efficiency, as well as technological progress resulting from 
technological changes in the subsector of electric power 
generation, have been to a similar extent responsible for 
the increase. The technological progress, a component 
of the productivity index, results in that the limit line 


2001/2000 2002/2001 2003/2002 
= MOG 1,046 1,006 1,006 
mzm.TE 1,001 1,000 1,012 
mzm.TG 1,047 1,007. 1,001 


2004/2003 
1,027 ] 
i = 
3,01: 


Fig. 1. Mean measures of the Malmquist index (MOG) and its components: 
the relative efficiency (zm. TE) as well as the technological progress (zm. 
TG), observed in the professional heat & electric power and electric power 
stations for the years 2000-2004. 
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of the set of production capabilities moves upward in 
the diagram. It tells about changes in general tendencies 
for the investigated group of entities. The mean yearly 
increase of productivity of the investigated professional 
thermal electric power and heat & power stations installed 
in Poland in the years 2000+2004, was equal to 2,1%. 
For comparison, in China the mean productivity increase 
for the years 1995+2000, measured by TFP index, 
amounted also to 2,1%, whereas in Australia was equal 
on average to 1,6% for the years 1981+1991. 


The observed mean yearly increase, on 2,1% 
level, of productivity of the investigated power industry 
enterprises was mainly caused by technological progress. 
The productivity increase of the domestic power industry 
sector resulted from employment rationalization as well 
as modernization of permanent assets, and installation 
of modern devices in electric power and heat & electric 
power stations. As a result, higher efficiency of their 
equipment, manifested by lower specific fuel consumption 
for power generation, was obtained. 


In the investigated period the productivity increase 
was also influenced, apart from atmospheric conditions 
and economy trends, by other factors such as political 
decisions as to rebuilding and modernizing the power 
sources installed in this country, as well as multi-year 
continuation of the long-term contract system not including 
any requirements as to introducing improvements to the 
power industry enterprises, that could not be conductive 
to taking care of increasing energy generation efficiency 
in the enterprises. 


Apart from the general changes in productivity of 
the sector in question also impact of internal factors on 
efficiency and time-dependent changes in productivity 
were analyzed for the investigated group of enterprises. 
The factors selected for the analysis were as follows: 
scale of enterprise activity, form of ownership, as well 
as type of carried-out activity. Verification of impact of 
scale factor on activity results shows whether the making 
use of a greater scale of activity by power generation 
enterprise does allow to achieve inputs savings and this 
way higher efficiency measures. The influence analysis 
of ownership form verifies the opinion whether private 
power industry enterprises are more efficient and faster 
developing than the state ones. And, the activity type 
factor, i.e. differentiation of enterprises with a view of 
prevailing share of either heat power or electric power 
production in a given enterprise, reveals differences in 
specificity of activity of the enterprises. 


In order to analyze impact of size of the considered 
enterprises on their efficiency, and of changes in their 
productivity they have been split into three groups: those 
of the total available heat and electric power: below 500 
MW, between 500 and 1500 MW, and over 1500 MW. 
Mean values of particular efficiency categories were then 
calculated and time-dependent changes in productivity for 
the particular groups were evaluated (see Tab. 2 and Tab. 3). 
As results from the analyses the enterprises of the total 
available power below 500 MW have achieved the 
lowest scale efficiency measures in the particular years, 
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however the mean measure of their scale efficiency has 
increased in the years 2000+2004. For the entities the 
scale efficiency constitutes an important element of the 
total technical efficiency, and the efficiency limit line in 
the case of assuming changeable scale effects, is the most 
distant from that in the case of constant scale effects. 
The group forms the main source of scale inefficiency 


Tab. 2. Technical efficiency versus size of power industry enterprise, 
measured by its total available output power. 


Enterprises of total available heat & electric E S VRS 
output power of 
ow 500 MW 87,46% 


Year 


2000 


500 MW 
0 MW 


Z 


83,34% | 92,06% | 90,50% 
87,27% | 88339 
95,25% | 97,52% 
83,64% | 91,27% | 91,75% 
87,65% | 88,21% | 99,31% 
94,55% | 97,15% | 97,28% 
88,92% | 92,16% | 96,52% 
85,29% | 85,82% 
91,73% | 94,04% 
90,13% | 92,62% | 97,37% 
86,82% | 87,44% | 99,24% 
92,40% | 94,63% | 97,63% 


Tab. 3. Productivity versus size of power industry enterprise, measured by 
its total available output power. 
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Tab. 4. Technical efficiency versus form of ownership of power generation enterprise 


m 
Polish state enterprises 

Polish state enterprises 
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Polish state enterprises 
olish private enterprises 


Malmquist 
index MG 
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2001/2000 500 MW +1500 MW 
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over 0 MW 
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2002/2001 03,1096 
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over 1500 MW 
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Tab. 5. Productivity versus form of ownership of power generation 
enterprise 


almquist Relative Technological 
index efficiency rogress 
MG (zm. TË) (zm. T° ) 


Category of enterprise 


Y 
3 


ish private enterprises 


b 
° 


ish state enterprises 


oreign enterprises 


| 


Y 
3 


ish private enterprises 


| 
=] 


sh state enterprises 
Foreign enterprises 
ish private enterprises 


ish state enterprises 


01, 
8, 
3, 


Foreign enterprises 


Y 
fo) 


ish private enterprises 


for the years 2000+2002. And, almost in the whole 
investigation period, i.e. in the years 2001+2004, their 
technical efficiency, in the case of assuming changeable 
scale effects, is higher than that of the enterprises of 
the output between 500 and 1500 MW, and lower than 
that of the enterprises of over 1500 MW output. The 
enterprises of the total available power between 500 and 
1500 MW are characterised, almost during the whole 
investigation period, by the lowest technical efficiency 
measures in the case of assuming changeable scale 
effects, however they achieve the highest scale efficiency 
Tab. 6. Technical efficiency versus type of activity carried out by 


enterpriseenterprise 
Category of enterprise M_VRS E_S_VRS 
93,81% 


2000 
T 
2001 


"v 
° 


ish state enterprises 


oreign enterprises 


| 


Y 
fo) 


ish private enterprises 


v 
° 


ish state enterprises 


Foreign enterprises 


M_CRS 


Heat & electric power stations and 


D 
their complexes 86:36% 


91,79% 


Electric power stations and their 


95,16% 
complexes 


99,39% 


Complexes of electric power 
stations and heat & electric power 90,13% 
stations 


Heat & electric power stations and 


D 
their complexes 86A 


93,65% 


92,76% 
95,25% 


88,74% 96,35% 
91,98% 
95,02% 


90,19% 
91,57% 
91,68% 


87,82% 
92,26% 
92,67% 


88,79% 


Electric power stations and their 


93,02% 
complexes 


97,68% 
Complexes of electric power 
stations and heat & electric power 85,32% 
stations 


Heat & electric power stations and 


o, 
their complexes 87.10% 


94,71% 
Electric power stations and their 


93,06% 
complexes 


97,96% 
Complexes of electric power 
stations and heat & electric power 
stations 


84,76% 94,14% 


2002 
2003 
2004 


Heat & electric power stations and 


D 
their complexes 92:538 


97,8196 
Electric power stations and their 


90,5096 
complexes 


98,6896 
Complexes of electric power 
stations and heat & electric power 
stations 


80,1696 91,79% 


Heat & electric power stations and 


o 
their complexes 90,69% 


98,34% 
Electric power stations and their 


91,4996 
complexes 


98,7096 
Complexes of electric power 
stations and heat & electric power 
stations 


81,0896 91,7196 


measures. The enterprises of the total available power 
of over 1500 MW are the most efficient group as their 
technical efficiency at scale effects assumed constant, 
as well as its main component, i.e. technical efficiency 
at scale effects assumed changeable, have achieved 
the highest values in the whole investigation period. 
As results from the analysis of Malmquist indices, the 
improvement of productivity has resulted from activity 
rationalization by means of increasing the relative 
efficiency of electric power and heat & electric power 
stations of the total available power below 500 MW, as 
well as due to significant technological progress obtained 
by the entities of the power output greater than 500 MW. 
To compare efficiency and time-dependent productivity 
changes in function of ownership forms the enterprises 
were divided into relevant groups. To this end Polish 
private enterprises, state enterprises and foreign ones 
were distinguished. On the basis of the results presented 
in Tab. 4 it can be stated that in the years 2000+2004 
Polish private enetrprises as well as foreign ones 
achieved on average higher measures of technical and 
scale efficiency. Generally the state-owned entities 
obtained lower estimates in particular categories of 
efficiency. However disproportions between the state and 
foreign enterprises were decreasing in the investigated 
period. As results from the analysis of productivity 
indices (see Tab. 5) the productivity of Polish private 
enterprises was increasing a little faster, and that of 
the state-owned and foreign ones — a little slower. 
The higher technical efficiency measures achieved by 


Tab. 7. Productivity versus type of activity carried out by enterprise 


" Relative Technological 
: Malmquist $ 
Year Group of enterprises index M G efficiency rogress 
0 (zm. TP) (zm. T° ) 
Heat & electric power 
stations and their 105,01% 101,47% 03,64% 
complexes 
Electric power stations 3 n 
2001/2000 and their complexes 103,49% 06,02% 


Complexes of electric 
power stations and " d 
heat & electric power 10583% 10,34% 
stations 
Heat & electric power 
stations and their 99,72% 99,73% 
complexes 
Electric power stations " " 
2002/2001 and their complexes 02:59% 102,69% 
Complexes of electric 
power stationis ang 02,46% 99,25% 103,36% 
heat & electric power 
stations 
Heat & electric power 
stations and their 00,80% 98,31% 
complexes 
Electric power stations S 
2003/2002 and their complexes 00,5778 
Complexes of electric 
power stations and x o 
heat & electric power 3334A 24,316 
stations 


03,66% 
05,51% 
Heat & electric power 
stations and their 101,65% 01,38% 
complexes 
Electric power stations ð " P 
2004/2003 and their complexes iG 10103% 00,81% 


Complexes of electric 
power. stations and 102,91% 101,01% 01,86% 
heat & electric power 
stations 
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97,63% 


94,9396 


100,12% 


99,95% 


103,35% 


97,08% 


102,89% 


Polish private and foreign entities could result from 
the fact that privatization was mainly focused on the 
enterprises of good or average financial situation hence 
they had better production potential right from the 
beginning. The below presented conclusions were drawn 
as a result of the performed comparison of the efficiency 
and changes in productvity of the electric power stations 
and their complexes and those of the heat & electric power 
stations and their complexes, i.e. impact analysis of type 
of carried-out activity on effects (see Tab. 6 and Tab. 7): 
In the investigation period of 2000+2004 the mean 
measures of technical efficiency as well as scale 
efficiency for the group „electric power stations and their 
complexes” were higher than those for the group ,, heat 
& electric power stations and their complexes”. However 
the superiority was decreasing along with time in the 
period in question. The superiority of the electric power 
stations could result from an advantage they take from 
large-scale production. The decreasing disproportions 
could result from growing importance of co-generation. 


Summary 


The occurring economy transformations, energy 
marketisation as well as ownership transformations of 
power industry enterprises have resulted in that power 
producers more and more attention have focused on 
the category of efficiency and its measurement. The 
improving of efficiency of the entities is very important 
for their functioning in conditions of power market 
consolidation and liberalization. Technical efficiency, 
an important element of economic efficiency, is decisive 
to the future of domestic professional thermal electric 
power stations. The performed verification of impact 
of various factors on the effects obtained by power 
producers provides us with interesting conclusions as 
to influence of privatization, scale of activity, type of 
carried-out activity on chvanges occurring in the sector. 
This work verifies to some extent common opinions 
dealing with influence of selling-up power generation 
enterprises to domestic private and foreign entities. 
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Optimization of steam cycles with respect to 
supercritical parameters 


Marian Piwowarski. 


Abstract 


This paper contains an analysis of supercritical power stations have been presently built worldwide. The analysis concerns 
cycles with single and double interstage superheating and, additionally, with regenerative feed water preheating system 
equipped with six, seven, eight, nine and ten regenerative heat exchangers, respectively. Relevant calculations were performed 
for various values of fresh steam temperature and pressure, various pressure values of secondary superheating, as well as with 
value of condenser internal pressure maintained constant. The calculations show that to increase efficiency of steam cycle with 
double interstage superheating and extended regeneration even to a value greater than 51%, is possible. In the age of greater 
and greater demand for electric power and stronger and stronger limitations imposed on emission of noxious compounds 
to the atmosphere, the developing of power production technologies based on supercritical parameters seems inevitable. 


Keywords: analysis of supercritical power stations, optimization of steam cycles, supercritical parameters, concerns cycles. 


Description of supercritical power stations 


Investigations on the increasing of efficiency of 
steam turbine cycle have been carried out by many 
research centres of universities and turbine production 
concerns already for a few dozen of years. In the case 
of thermodynamic cycles with applied steam turbines 
the investigations have dealt with steam pressure and 
temperature values greater than: the 22.12MPa critical 
pressure and 647.28K critical temperature, respectively. 
Over the critical point the difference between liquid 
phase and gas phase disapears [20], [25], [43]. Already 
in the 1950s the first attempts to developing the turbine 
power units based on supercritical parameters have 
been made. In 1957 the first in the world turbine power 
unit of supercritical parameters was put in motion 
in the electric power station, Philo, USA. The unit of 
314bar/621?C steam parameters developed 125MW 
output power. Two years later 335MW power unit of even 
higher steam parameters (345bar/649°C/566°C/566°C) 
was built in the Eddystone I electric power station. For 
the high values of steam parameters special materials 
suitable to work in high temperatures were required. 
High investment cost, inadequate quality of materials, 
faults in assembling and operational problems forced 
the producers to resign from building the units for 
so high steam parameters. As late as in 1969 the unit 
designed for the supercritical parameters with double 
secondary superheating of 241bar/538*C/552?C/566?C, 
was set working. The application of such solution to 
the steam cycle has improved operation of the unit. 
Therefore in the 1970s, apart from the units with single 
superheating of 241bar/538/538°C steam parameters, 
were built the units with double superheating, which 
developed output power reaching from 350 do 1100 MW. 
Progress in material technology, started in the 1980s, 
has made it possible to apply higher steam parameters 
(310bar/593°C/593°C/593°C) to steam cycle. The 
application of better and better materials suitable for 
very high operational temperatures makes it possible 
to build the units for the supercritical parameters as 
well as ultrasupercritical ones (400bar/760°C). Modern 
3-D calculation programs used in designing fluid flow 


systems as well as highly efficient devices included into 
equipment of such units have guaranteed to achieve their 
high reliability and efficiency. In the 1960+1990s in USA 
159 units of the power range of 300+1400MW, pressure 
range of 230+260bar and temperature of 540+590°C, 
including 14 units of double secondary superheating, were 
set working. In the years 1990-1998 in Japan and China 
total output power installed in coal electric power stations 
has increased threefold. The growth was achieved 
mainly by applying the units of output power in the range 
of 400+700MW, working with the supercritical steam 
parameters of 255bar/570+590°C and the secondary 
superheating of 570+595°C temperature range. In the 
years 2000-2003 three 700MW units based on the 
246bar/593°C/593°C steam parameters, two 900MW 
units based on the 241 bar/593°C/593°C steam parameters, 
two 900MW units based on the 241bar/600*C/610*C 
steam parameters, three 1000MW units based on the 
245bar/600°C/600°C steam parameters and two 1050MW 
units based on the 250bar/600°C/610°C steam parameters, 
were put in operation. In Europe significant achievements 
in building the power units with supercritical parameters 
can be noted in Denmark and Germany. In 1984 in 
Denmark the first unit with supercritical parameters was 
set working in the Studstrupvaerket electric power station. 
At the beginning of the 21st century two 411MW units 
with the 290bar/582°C/580°C/580°C steam parameters, 
including one coal-fired and the other gas-fired, as well 
as 530MW unit with the 300bar/580°C/600°C steam 
parameters, working on combusted biomass, were put 
in operation in Avedore electric power station. Germans, 
basing on the achievements and operational experience 
of Danes in the area of building coal electric power 
plants with high steam parameters, have built mainly 
800+900MW power units fitted with brown coal-fired 
boilers. In 2002, 1012MW power unit based on the 
274bar/580°C/600°C steam parameters was built in the 
Niederaussem electric power station. A list of selected 
power units (presently installed or planned ones) 
intended for the operating with supercritical parameters 
is presented in Tab. 1. Patnów II electric power station is 
the first Polish electric power plant based on supercritical 
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Fig. 2 Cross-section of 28K460 turbine installed in Łagisza electric power station 


parameters. Its high energy conversion efficiency is 
associated with lower fuel consumption and limited 
emission of contaminations to the environment. And, the 
power unit built in Lagisza is fitted with the biggest in the 
world monotube boiler of circulation fluid-bed, operating 
with supercritical parameters. The modern ecological 
power unit would reach very high efficiency equal to 
about 45%. Its turbine-set parameters are as follows: 
electric output power of 460MW, fresh steam pressure 
of 275bar, fresh steam temperature of 560?C, secondary 
steam temperature of 580°C, secondary steam pressure 
of 54.6bar. The steam cycle schematic diagram as well 
as the turbine-set cross-section are presented in Fig. 1 
and 2, respectively. 


Fig. 1 Simplified steam cycle diagram of Lagisza electric power plant 


And, in the electric power station of Betchatow S.A., 
the largest conventional one in Poland and Europe, 
working on brown coal and being the biggest electric 
power producer in Poland, the building of a 833MW 
power unit based on supercritical parameters has been 
started. The unit with the 266bar/554/582°C steam 
parameters is intended for the fundamental mode of 
operation of 7500 h/year usage time of its rated power 
and the total working time of 8100 h yearly. By applying 
supercritical steam parameters it is already possible to 
obtain thermal efficiency of over 45% as compared with 
that of 41+42% achievable today by conventional thermal 
electric power plants. Contemporary pulverized-fuel 
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supercritical boilers deliver the steam of 25MPa pressure 
and 600°C temperature and abt. 610°C superheating 
temperature, to steam turbine cycle. It is expected that 
further development of the technology would be focused 
on the mastering of ultrasupercritical parameters as well 
as the increasing of efficiency of power units. This mainly 
depends on progress to be done in the area of material 
engineering [1], [2], [3], [5], [7], [9], [13], [15], [16], [17], 
[22], [28], [29], [30], [31], [32], [34], [39], [40], [48], [49], 
[50]. Investigations are also carried out in the frame of 
large international projects such as e.g. THERMIE 700 
Advanced Power Plant project financially supported by 
EU, whose simplified scheme is given in Fig. 3, aimed 
at the obtaining of fresh steam temperature of the order 
of 700°C and 37.5 MPa pressure. The setting in motion 
of the power unit is scheduled on 2015 [16]. 


> > 
235 "bar 


Fig.3 Simplified schematic diagram of THERMIE 700 steam cycle 


In the subject-matter literature can be found power 
plants fitted with turbines operating with supercritical 
parameters, in which fresh steam (of 26MPa and 570°C) 
is produced by cooling nuclear reactors. Efficiency of 
such systems exceeds 4596, as compared with that of the 
order of 33% in the case ofthe power plants co-working 
with light-water—cooled reactors [4], [38]. 


High steam temperatures in steam turbine cycles 
require the first and second stage of turbine placed behind 
the secondary superheating, to be cooled. The simplified 


Tab. 1. List of selected power units with supercritical parameters, installed in electric power stations worldwide 


[1], [2], [7], [9], [13], [21], [25], [29], [36], [37], [44], [45] 
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schematic diagram of the external cooling is shown in 
Fig. 4. The cooling steam taken before the first 
superheating is directed to the high-pressure (HP) part, 
to be mixed with the superheated steam. The cooling 
steam for the first stages of the intermediate-pressure 
(IP) turbine is taken from the third or fourth stage of HP 
turbine. 


- 


Application of the high parameters is associated with 
necessity of using special materials. Elements of HP and 
IP turbines, which are in direct contact with superheated 
steam, are made of a high-alloy steel (of high percentage 
content of such elements as Cr, Mo and V). The first 
stages of HP and IP turbine are made of Ni-alloy steel 
which is resistant to thermal loads resulting from the 
difference between temperature of superheated steam 
and that of steam after expansion. The low-pressure (LP) 
turbine module is made of NiCrMoV-high-alloy steel. 
The material is brittle-crack resistant due to lowered 
content of such elements as P, Sn, Mn and Si and 
increased Ni content. For the LP part it is very important 
to appropriately design the last stages exposed to high 
loads and erosion. 


> 


pean 


Fig. 4 Schematic diagram of external cooling 
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In present for steam temperatures up to 600°C ferritic 
steels, e.g. SAVEI2 steel, are used; for temperatures 
over 650°C — austenitic steels, e.g. SUPER304H steel 
and HR6W steel; and for temperatures in the range 
of 620+720°C — various alloys, e.g. Inconel 617, 
Haynes 230, Inconel 740 etc. Tab. 2 shows chemical 
element composition of selected kinds of materials used 
for the units of supercritical parameters. 


To produce the steam of suprecritical parameters, 
monotube boilers with elements made of special steels 
(e.g. high-temperature creep resisting, austenitic or 
ferritic-martensitic ones), are used. For its flow system 
a turbine-driven water-feed pump of high elevation 
head is necessary to overcome high flow drag resulting 
from higher density of medium and smaller internal 
diameters of pipes in the boiler. Another kinds of steam 
generators with supercritical parameters are atmospheric 
fluid-bed boilers as well as those of circulation fluid-bed. 


Boilers of pressurized fluid-beds make it possible to 
reach higher power concentration and higher efficiency 
values. Circulation fluid-bed technology in association 
with supercritical steam parameters constitutes a safe 
solution, more technologically and economically effective 
than the option of pulverized-fuel boiler. Advantages of 
the technology are: fuel flexibility, possible co-combustion of 
coal, sludge and biomass, compliance with WE80/2001/EU 
directive in the area of pollution emission, limitation of 
hazard to occurrence of high-temperature corrosion and 
erosion, increasing power unit’s cycle effciency due to 
heat recovery and lowering exhaust gas temperature, 
uniform distribution of heat flow in combustion chamber 
as well as improved dynamics of load changes [2], [3], 
[7], [8], [10], [11], [12], [14], [22], [23], [24], [27], [33], 
[41], [42], [46], [47]. 


Results of calculations of steam cycles 


The increasing of fresh steam pressure and temperature 
at steam pressure in condenser kept constant, is the first 
step to the achieving ofa higher cycle efficiency [6], [18], 
[19], [26], [32]. Successive figures present respectively: 
Fig. 5 — relation of ideal efficiency, Fig. 6 — real efficiency, 
and Fig. 7 — cycle efficiency -with wetness loss taken 
into account -of Clausius-Rankine cycle, all in function 
of fresh steam pressure (5+55MPa) and temperature 
(500+760°C) at condenser internal pressure kept constant 
(4kPa); and, in each of the figures the dryness degree 
limit line (xgr=0.85) is depicted. The figures illustrate 
the possibly obtainable increase of efficiency in function 
of values of fresh steam parameters at inlet to turbine, 
at condenser internal pressure kept constant. In the 
diagrams can be distinctly observed the limitation of 
upper value of steam initial pressure, resulting from the 
limit wetness. The C-R cycle ideal efficiency obtainable 
due to application of high fresh steam parameters exceeds 
52% (see Fig. 5), whereas the real efficiency of the cycle 
does not exceed 43% (see Fig. 6). It should be also noted 
that the real efficiency of the cycle in which steam wetness 
degree has been taken into account does not exceed 39%, 
and that C-R cycle optimum pressure values amount to 
about 100 bar (see Fig. 7). The applied here notion of the 
taking into account of wetness degree means that work 
is done only by steam and that water does not provide 
any work [26]. 


Application of the interstage superheating is the next 
way to increase efficiency of the cycle with steam turbine. 
Its efficiency depends, apart from fresh steam parameters 
and condenser internal pressure, on superheating pressure 


Tab. 2. List of chemical element composition of materials used for the power units of supercritical parameters [14] 
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Fig. 5 Relation of ideal efficiency of Clausius-Rankine cycle in function of 
fresh steam pressure and temperature at condenser internal pressure kept 
constant (4kPa), and the wetness degree limit curve xgr=0.85 


value too. There is one optimum pressure value for which 
the cycle efficiency reaches its maximum. In steam power 
units various configurations of the system with single 
interstage superheating applied in the boiler are possible. 
The solution in which HP part and IP part are placed in 
a common casing, is often applied. The integrated HP/IP 
module is connected with LP part which can be of single 
or double jet. Such system is characteristic of compact 
structure and lower investment and operational costs. 
The units which operate with supercritical parameters 
and have that kind of structure can reach output power 
of over 600 MW. The system in which HP and IP turbine 
parts are placed is seperate casings is only sporadically 
applied. In the units of a greater power the solution in 
which single-jet HP part and double-jet IP part are placed 
in separate casings and connected with double-jet LP part, 
1s applied. The LP module is consisted of one, two or three 
casings, depending on the output power the unit has to 
develop. In practice, finds also application the system 
in which the first casing contains HP part and single-jet 
IP part of countercurrent flow, and in the other casing 
double-jet IP part connected with LP module consisted 
of two casings, is placed. Its alternative is to place 
HP and IP parts in separate casings [2], [6], [9], [15], [18], 
[19], [26], [32], [33], [35], [36]. 


The profits resulting from the application of double 
superheating became obvious as early as in 1960s. In 
optimizing the cycle with double secondary superheating 
attention should be paid to appropriate choice of the 
pressure for the first and second superheating. 


In practice the Ist superheating pressure is selected 
depending on thermodynamic optimum, the 2nd 
superheating pressure is usually chosen depending on 
an assumed steam temperature at inlet to LP turbine. 
Maximum temperature of inlet steam to LP part is limited 
with respect to thermal strength of materials. Classical 
turbine-set consists of three separate modules designed 
for definite steam parameters, 1.e.: HP, IP, and LP parts. 
Electric generator is directly connected with the last part, 
i.e. LP. The arrangement of steam turbine for supercritical 
parameters depends first of all on choice of a kind of 
secondary superheating, unit’s operation range as well 
as special requirements as to regenerative preheating. To 
the arrangement with double superheating the solution in 
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Fig. 6 Relation of real efficiency of Clausius-Rankine cycle in function of 
fresh steam pressure and temperature at condenser internal pressure kept 
constant (4kPa), and the wetness degree limit curve xgr=0.85 
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Fig. 7 Relation of Clausius-Rankine cycle efficiency (with steam wetness 
degree taken into account according to Bauman formula) in function of 
fresh steam pressure and temperature at condenser internal pressure kept 
constant (4kPa), and the wetness degree limit curve xgr=0.85 
which HP part is placed in a separate casing and connected 
with casing of double-jet IP part after 2nd superheating, is 
often applied. LP part is consisted of one or two casings. 
In order to obtain higher power values the solution in 
which HP parts and IP1 part after 1st superheating are 
placed in common casing, is used. The second casing is 
intended for IP2 double-jet part after 2nd superheating, 
connected with LP module consisted of one, two or three 
double-jet parts [2], [6], [9], [15], [18], [19], [26], [32], 
[33], [35], [36]. 


Increase of cycle efficiency to a large extent depends on 
an appropriate choice of number and kind of preheaters. 
Regenerative preheating is a very important element of 
the power unit and it influences the unit's main elements, 
i.e. boiler, turbine, condenser and feed pump. The 
increased jet of steam produced in boiler, resulting from 
regeneration, is associated with the necessity of fitting the 
unit with feed pumps of greater capacity. High feed water 
temperature makes it difficultto maintain low exhaust gas 
temperature on which to a large extent depends boiler 
efficiency. It should be also remembered that along with 
increasing number of exchangers degree of complexity of 
the entire system also increases, that consequently leads to 
increasing investment cost. The first step in optimizing the 
system is to select an appropriate number of regenartion 
stages. To high power units 6-10 stages of superheating 
are usually applied. When considering profits due to 


POLISH MARITIME RESEARCH, Special issue 2009/51 49 


application of regenerative superheating attention should 
be first of all paid to possible achieving higher efficiency 
of turbine stages. Drop of steam jet in LP part makes blade 
system forming easier and, on the other hand, absorption 
of steam from interstage space of the turbine makes its 
design more complex and results in generating flow losses 
within the turbine [6], [18], [19], [26]. 


In this work the analysis is performed of the 
cycles with single and double interstage superheating 
and, additionally, with the regenerative feed water 
preheating system fitted with six, seven, eight, nine 
and ten regenerative heat exchangers. Respective 
calculations were performed for various values of fresh 
steam temperature (500°C, 600°C, 700°C) and pressure 
(5+65MPa), various pressure values of secondary 
superheating: p,=(0.24+0.36)*p,, and p,.=(0.06+0.12)*po, 
as well as with internal pressure value in condenser 
maintained constant (4kPa). The performed calculations 
indicate that the increase of efficiency of the steam 
cycle with double interstage superheating and extended 
regeneration, up to 51%, is possible. 


Comparison of cycle efficiency 
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Fig. 8 Comparison of efficiency of steam turbine cycles of various 
configurations, including values for Polish electric power stations operating 
with supercritical parameters 


Fig. 8 presents comparison of efficiency of the 
following cycles: 


* Clausius-Rankine cycle (marked C-R); 


* The steam cycle with single interstage superheating 
without regeneration (marked 1-prz); 


* The steam cycle with single interstage superheating 
and 8 regenerative heat exchangers (marked 
1-prz 8w); 

* The steam cycle with double interstage superheating 
without regeneration (marked 2-prz); 


* The steam cycle with double interstage superheating 
without regeneration and 8 regenerative heat 
exchangers (marked 2-prz 8w); 


* The steam cycle with supercritical parameters 
operating in Patnów electric power station (marked 
Patnow); 


* The steam cycle with supercritical parameters 
operating in Lagisza electric power station (marked 
Lagisza); 


* The steam cycle with supercritical parameters 
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operating in Belchatów electric power station 
(marked Betchatow); 


* The demonstrative cycle for supercritical 
parameters (marked Thermie 700). 


Fig. 8 was prepared for optimum values of interstage 
superheating pressure, fresh steam temperature (of 700°C) 
and constant value of condenser internal pressure (of 
4kPa) in function of fresh steam pressure (in the range 
of 20+60MPa). From the diagram it can be concluded 
that the application of double interstage superheating 
increases the cycle efficiency by about 1.5 % as compared 
with that ofthe cycle with single single superheating and 
by about 9 % as compared with the simple C-R cycle. 


Summary 


In Poland coal plays the most important role in the 
process of electric power production. With a view of 
it resources, gained experience and reliability of the 
coal-based technology of electric power production that 
fuel will be dominating for electric power generation 
in the years to come. Ecological and energy policy, 
both in Poland and EU, compels to apply low-emission 
technologies, e.g. clean coal-based technology. In Poland, 
because of its coal resources and its role for the state's 
energy balance, investments in clean coal technologies 
should be a natural phase of power industry development 
in this country. Other probable development directions 
of Polish power industry are a.o. the following: coal 
gasification integrated with high-temperature fuel cells, 
power systems with fuel cells combined with coal hydro- 
gasification, cola gasification and liquefaction, coal- 
nuclear synergy systems, pressurized coal combustion 
in fluid-bed boilers, combined production systems of 
electric power and hydrogen, or polygeneration. In the 
age of greater and greater electric power demand and 
stronger and stronger limitations imposed on emission 
of noxious compounds to the atmosphere, development 
of power production technologies based on supercritical 
and ultrasupercritical parameters, seems inevitable. 
The turbine power units fed with steam of supercritical 
parameters, which have been built so far, are characteristic 
of a higher efficiency and lower carbon dioxide emission, 
resulting from it. 
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A method for analysing ram pressure 
characteristics of impeller pump rotor 


R. Puzyrewski, 
K. Zochowski, 
P. Flaszyriski. 


Abstract 
This paper presents a method in which typical tests of centrifugal pump are used to obtain information on real value of discharge 
angle of flow leaving the rotor. The method can be applied to properly choose inlet angle to blade palisade of centrifugal guide 
vanes in the case when to perform measurements of velocity fields behind the rotor more precisely is not possible. 


Keywords: impeller pumps, characteristics, tests 


Introduction 


Characteristics ofimpeller machines such as efficiency 
or ram pressure of experimental stages, are determined 
on the basis of certain definitions and measurements 
of relevant quantities, performed on a tested machine. 
This work deals with the following characteristics: ram 
pressure of radial pump rotor, as well as liquid flow torque 
generated by pump rotor, determined on an experimental 
test stand. The characteristics were determined either at 
constant rotational speed or at constant volumetric flow 
rate and varying rotational speed. The problem has been 
formulated as follows: if and in which way the rotor 
channel characteristics in the form of triangle of velocity 
vectors at outlet from the rotor can be determined on 
the basis of experimentally obtained ram pressure and 
torque characteristics. An important parameter belonging 
to the rotor characteristics is the outlet angle of flow 
leaving the rotor, which is usually much different from 
the rotor blade outlet geometrical angle. For this reason 
the erroneous mating of rotor to blades of guide vanes 
channels often occurs. The correct mating is conditioned 
by information on the outlet angle from rotor and inlet 
angle to guide vanes channel behind the rotor. This 
paper presents a method to obtain such information by 
measuring the rotor characteristics on the tested pump. 
The problem is in opposition to the method commonly 
used for determining the pump rotor characteristics 
on the basis of the characteristics of rotor channels. 


Basic relations 


Fig. 1 presents schematic diagram of a test stand on 
which tests of a model pump were performed. The stand 
has made it possible to measure basic quantities necessary 
for determination of the characteristics and to perform 
their subsequent analysis. 


The basic relation intended to be used in analysing 
measured quantities is the energy conservation equation 


2 2 
SL gap E (1) 
2 p 2. p 


in the absolute coordinate system: 
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where: 
a -energy delivered to rotor, 
ep m -kinetic, pressure and internal 
p 2 energy at inlet (index ‘1’) and at 


outlet (index ‘2’), respectively. 


For single rotor the potential energy differences 
gz, =gz, are usually very small or can be taken 
into account by introducing the piezometric pressure 
ptpgz to Eq. (1). The applied notation complies with that 


used in classical technical literature on pumps [1], [2]. 


If to introduce the total pressure increase measured 
on the model machine (marked ‘m’): 


e; ci 
Ap, = ^ * i [» * 2 (2) 
then Eq. (1) can be expressed as follows 
a= a + Ae (3) 


In Eq. (3) energy dissipation (losses) is 
determined by the increase of internal energy: 


Ae =e, -e >0 
Eq. (3) states that the energy delivered by the 


rotor is transformed to the total pressure increase 
and energy losses inside the rotor. 


Fig. 1. Schematic diagram of model pump test stand (throttle valve not shown). 


Notation: pt — forcing pressure, ps — suction pressure, hzw — pressure 
difference on measuring orifice plate, F — force of action of driving motor 
pivotely fixed stator, n — rotational speed. 


According to Euler formula, for axial inflow 
direction to rotor the delivered energy is expressed 
by the relation: 


a-uc, (4) 


determined by the rotor discharge parameters 
(index ‘2’ assigned to the circumferential velocity 
u, and its projection to circumferential direction, c 
has been here omitted). 


u2? 


The energy conservation equation: 


u 


^ 
uc c3 t Ae (S) 


is expressed in the form in which appear two 
measured quantities: 


E nD, nm 
60 


and y 


AD 


where the rotational speed n, is measured. In 
Eq. (5) the quantities c, and Ae are unknown. 
For the torque M_ measured at the pivotely fixed 
stator of the electric driving motor of the pump, the 
following relation can be written: 


Mo l e,'Q, *M, (6) 


where: 
Q. -measured volumetric flow rate; 
M, -moment of friction of rotor disk and cover 
against liquid. 


In Eq. (6) the quantities c, and M, are unknown. 


By analysing the following characteristics: 


Ap =p(Q), where n =const (7) 
M,,= M(Q,,), where n,,=const (8) 


and 
Ap,, = Ap(n,,), where Q,,=const (9) 


M,-Mí(n,) where Q,-const (10) 


it is possible to determine parameters of the 
velocity triangle, under additional assumptions 
as to the velocity triangle at outlet from the rotor. 
Like in the above mentioned case, with exception 


N u; C 


ie 
Fig. 2 Velocity triangles at rotor inlet; for n,,=const and Q,, = var. 


of the characteristics (10), changes of M, values 
proportional to n, ? with an unknown proportionality 
coefficient, will be equivalent to changes of the 
rotational speed n,. For n, = constand varying Q, 
the velocity triangle shown in Fig. 2. 


The meridional velocity c, in the rotor outlet 
cross-section of the breadth h, at the diameter D,, is 
determined by the formula: 


c Qha 


= (11) 
mD,h, 


If to assume that change of Q. does not cause 
any significant change of the angle p then the 
differentiation of the relation: 


Q m 


tg B — ———c7—— —-const (12) 
i P mD,h, (u-c, ) 
with respect to Q.. yields: 
dc, _ C, -u (13) 
dQ. Qn 


For Q =const (c, —const) and varying u (n, —var) 
the velocity triangles are such as shown in Fig. 3. 


Fig. 3 Velocity triangles at rotor outlet; for Q,, 


—const and n m7 var 


Itis easy to observe that in this case the following 
yields from Eq. (12): 


(14) 


Eqs. (13) and (14) will be further used for analysis 
of the characteristics. 


Analysis of the ram pressure characteristics 
Ap, -Ap(Q,). 


In Fig. 4 is presented the pump rotor characteristics 
Ap=f(Q) where: Ap [Pa] -ram pressure rise in 
function the volumetric flow rate Q [m?/h] at the 
rotational speed n-const, obtained by means of the 
measurements carried out on the model test stand. 
In Fig. 4 the indices ‘m’ appearing in the notation 
of coordinates, are omitted. 
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Fig. 4 The measured pump rotor characteristics Ap = f(Q) 


The measurement points were interpolated by 
using a 3rd order curve which yielded a satisfactory 
approximation within the range of flow rate of 40+70 
[m?/h]. 

Eq. (5) in the form: 
p:u-c,=Aptp -Ae (15) 
after differentiation with respect to Q, yields: 

" dc, dAp $p dAe (16) 
dQ dQ dQ 
By making use of Eq. (13) the following equation 

is obtained: m 

uc =p u+ p 
puc,-pu *Q do 
which, after dividing by puc,, makes it possible 


to estimate the following parameter crucial for the 
velocity triangle: 


C 1 


U a 


p 


dAe 
+ (17) 
pQ do 


18 
uU oq Q dAp pQ d^e e» 
puc, dQ puc, dQ 
By introducing the approximations: 
puc, =Ap (19) 
and: d 
gae i (20) 
dQ 
the following formula is achieved: 
C VL 1 
u , Q dAp 
Ap dQ 
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Eq. (21) makes use of information taken from the 
characteristics, inclusive of derivative in a selected 
point of it. 

After determination of the quantities: 

ie mD,n Q 


and c„= 
mD,h, 


(22) 


for known values of: the flow rate Q, rotational 
speed n, rotor outlet diameter D,, and breadth h,, the 
following angles appearing in the velocity triangle, 
can be determined: 
a=AreTan *- and B-ArcTan $ (23) 
C u-c 


u u 


Analysis of the rotor torque characteristics 
M, -M(Q) 


As the pump motor stator is pivotely fixed there 
is possible to measure the rotor - generated torque. 
Example results of the torque measurements are 
depicted in Fig. 5 where the experimental points has 
been interpolated by using a 2nd order curve. For 
their analysis the relation (6) will serve as the basis. 
Let's observe that at the constant rotational speed 
n and only the flow rate Q variable the following 
approximation seems to be good: 


dM, _ 
dQ 
Hence by differentiating the relation (6) the 
following expression is obtained: 


Q4) 


dM D, dc 
Lg hui (25) 
do ^2 (e, *Q ao? 
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Fig. 5 The measured pump rotor characteristics M = f(Q) 


From Eq. (25), by making use of Eq. (13), the 
following formula can be obtained: 


Mop D eea) 29) 
dQ 2 
from which it yields: 
u 2 pD,u dQ 


The right hand side of the relation is determined 
on the basis of the experimental characteristics M(Q) 
and measured values of the rotational speed n, . The 
angles are determined from the expressions (23). 


Analysis of the ram pressure characteristics 
Ap,,-Ap(u). 


A similar analysis can be performed for the 
pump rotor characteristics Ap=f(u) at Q=const. 
The characteristics are shown in Fig. 6 where the 
experimental points were interpolated by using a 2nd 
order curve which provides a good approximation. 


dc, ldAp,dAe 
u = ar 


(28) 
" du 


p du du 


Differentiation of the relation (5) with respect to 
u gives: 


u u dAp , pu dAe 


1+ (29) 
c, peu du pce,u du 
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Fig. 6 The measured pump rotor characteristics Ap — f(u) 


By making use of Eq. (14) the above given 
relation can be transformed into the following: 


c 1 


u D u dAp —] 
Ap du 
Finally the approximation (19) and the assumption 


on a small change of losses resulting from u (see 
Eq.(20) yields the formula: 


(30) 


whose right-hand-side values can be determined 
from the characteristics such as given in Fig. 6. 
Next from the expressions (23) the angles can be 
determined. 


Discussion of results of analysis of experimental 
characteristics 


The calculation results of the outlet angles a 
and p for particular characteristics, are presented 
in Figs. 8 through 10. The values obtained from the 
calculations according to Eqs. (21), (27), (30) are 
burdened both by measurement systematic errors 
and simplifying assumptions. 


The characteristics Ap=Ap(Q) are ,, contaminated” 
by backflow drag behind the rotor and guide vanes 
drag in the centripetal channel behind the rotor. In 
the conditions of the performed experiment it was 
not possible to obtain seperate characteristics for the 
rotor only. Drag rise behind the rotor at greater values 
of Q leads to lower values of measured ram pressures, 
that - in consequence - gives greater values of the 
angle a. If behind-the -rotor drags were reduced and 
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Fig. 7 Comparison of the experimental and hypothetical pump rotor 
characteristics 
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a greater ram pressure rise were this way obtained 
as it indicate the hypothetical characteristics shown 
in Fig. 7, then good comformity between the angles 
a calculated from Eqs. (21) and those obtained 
from torque measurements, acc. Eq. (27), would be 
achieved. The comparison of the values calculated 
according to the formula (21) for the experimental 
and hypothetical characteristics demonstrates 
influence of the error residing in the characteristics 
on the determined values characterizing the velocity 
triangle. In Fig. 8 the point determined on the basis 
of the calculations performed by using a 3D code, is 
presented. This is the value of the angle a, calculated 
as the surface average. The value relatively well fits 
those determined from the characteristics. 


The determined angle amakesitpossibleto calculate 
value of the rotor outlet angle B which constitutes 
an important feature of rotor blades palisade. 
Fig. 9 presents the calculation results. The change of 
the angle p, observed in the diagram, demonstrates 
an approximate character of the assumption on its 
invariance. Only variability ranges for Q — values, 
within which the changes would be sufficiently 
small, can be indicated. 


Analysis results of the characteristics Ap-Ap(u) 
are given in Fig. 10 and11. The experimental point 
shown in Fig. 10 is situated relatively close to the 
point obtained from the analysis of experimental 
characteristics. The small changes of the angle p, 
seen in Fig. 11, seem to satisfy the assumption 


on its invariance, better. However, by comparing 
30 
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Fig. 8 Comparison of calculation results of the discharge angle a, obtained 
from experiment with those obtained from 3D-calculations 
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Fig. 9 The discharge angle p, calculated on the basis of experiment 
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Fig. 10 Comparison of calculation results of the discharge angle a acc. Eq.(30) 
with those obtained from numerical experiment 


25 
[Us 
20 
acc. Eq. (30) e 
15 FS 5 e. e . 
10 
5 


3000 
n [rpm 


2500 2600 2700 2800 2900 


Fig. 11 Values of the discharge angle p, calculated acc. Eq. (30) 


Fig. 9 with Fig. 11, the discrepancy in the angle p 
values of the range from 14° to 20°, was revealed. 
The value ofthe geometrical angle ofthe rotor outlet 
edge amounted to about 30?. 


The angle a determined from the analysis makes 
it possible to more properly design the inlet to guide 
vanes channel behind the rotor for a given nominal 
value of Q. 


Final remarks 


The presented method for determining the 
velocity triangle parameters on the basis of the pump 
rotor characteristics can be used in the case if other 
possibilities to measure flow parameters are lacking. 
Information on the angle of flow discharge from the 
rotor is especially important. It much differs from 
the rotor blade geometrical outlet angle. This is an 
important parameter for correct forming the inlet 
to guide vanes ring behind the rotor. In this paper 
attention has been also drawn to the necessity of 
exact determination of discharge pressure from the 
rotor. In the described experiment the pressure was 
contaminated" by flow drag occurring in the space 
between the rotor outlet and the forcing pressure 
measurement point. For this reason the experiment 
was not sufficiently ,,pure” for the presented method. 
Nevertheless its relatively good conformity with the 
experiment was achieved at some scatter of values 
of the determined angles. 


The analysis demonstrates in which sense 
the rotor characteristics constitute ,,carriers” of 
information on flow parameters. 
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Nomenclature 
a  — work determined by Euler equation; 
c  —absolute velocity; 
e internal energy; 
h -—rotor breadth; 


n -rotational speed; 
p - pressure; 


— circumferential velocity; 


c 


z  —location height; 

D -rotor diameter; 

M - torque; 

Q -volumetric flow rate ( volumetric flow capacity); 


a  —angle between directions of the absolute 
velocity (c) and the circumferential velocity (u); 


P ^ - blade angle; 
A -rise of a quantity; 


p liquid density; 


Indices 


m 


2  — parameters at outlet edge of rotor blade; 
m - parameters of model machine; 


u - projection towards direction of circumferential 
velocity; 


t — quantities associated with friction. 
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— (stands for) parameters at intlet edge of rotor blade; 
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Remarks on aerodynamic forces in seals of 
turbine stages 


Robert Stępień, 
Krzysztof Kosowski 


Abstract 
This paper presents results of numerical examination of flow through over-shroud seals of turbine stages. Various 
labyrinth seals of different configurations and number of sealing teeth were considered. It was demonstrated 
that results of investigations of isolated seals cannot be directly used for analyzing turbine stage operation. Such 
approach may lead to relatively large errors in determining value of aerodynamic force and direction of its action. 


Keywords: rotor dynamics, self-excited vibrations, aerodynamic forces, blade seals, turbine shrouds 


Introduction 


Turbine rotor often moves along a complex trajectory 
during operation, usually it is displaced by the eccentricity 
e relative to turbine cylinder centre, it changes its 
logitudinal position a, and its axis can be situated aslant 
relative to cylinder axis, forming with it the angle ọ 
(Fig.1). This results in circumferential and axial change of 
clearance in turbine stage seals (over-shroud ones as well 
as external glands), that in turn causes that distribution of 
flow both in seals and grids of blades becomes no longer 
circularly symmetrical. Unbalanced aerodynamic forces 
appear (called blade-ring forces if occur in passages of 
grid of blades, and pressure forces if in seals). The forces 
can lead to self-excited vibrations of turbine rotors and 
-in some cases -prevent turbine form developing its rated 
power. 


For more than fifty years the phenomenon has been 
investigated in many research centres worldwide. Many 
results of experimental, theoretical and numerical 
investigations have been so far published. In the 
beginning, analyses based on the so called ,bulk-flow 
theory" have been developed ([11]-[13], [15]). In the 
method, simplifications are implemented in solving 
three-dimensional, non-stationary equations of flow 
through the clearance, consisting in linear association 
of shear stresses on the walls with average flow velocity 
(that is only possible at rather small rotor displacements 
against casing). This way flow behaviour equations 


Fig. 1 Changes in position of rotor against casing, a — change of axial 
clearance, e - eccentricity, 9 - angle of aslant displacement of rotor axis 
against casing axis 
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were simplified to 2*4 order ones as well as it was made 
possible to consider the flow to be axially symmetrical 
and stationary and that in which only small disturbances 
described by linear relations occur. As a result of solving 
such set of equations distributions of velocity and pressure 
fields in the clearance are achieved; they are then used 
for determining forces and moments generated in the 
clearance. The end of the 1980s and the 1990s brought 
development of computer calculation techniques used 
also in the modelling of flows through the clearance 
between rotor and casing, as well as in the determining 
of dynamic coefficients for such systems. The simulations 
are based on three-dimensional Navier-Stockes equations 
with Reynolds averaging (RANS), closed by means ofthe 
k-e turbulence model. For calculations algorithms based 
on the methods of finite elements [7], finite differences 
[9], [14], [17], [18] and finite volumes [1]-6], [8], [10], 
are used. Analysis of forces generated in isolated seals 
is only of a limited usefulness for design practice. Flows 
through seal are also influenced by phenomena which 
occur both in passages of vane blades and rotor blades 
grids. Therefore results obtained from investigations of 
isolated seals cannot be directly applied to analysis of 
operation of a complete turbine stage. 


Vane blades grid 


Effective diameter 374 mm 
Blade height 19 mm 
Chord of airfoil section 74 mm 
Pitch of grid 41,9 mm 
Number of blades 28 
Blade-angle setting 39°42’ 
Outlet angle 9°42” 
Rotor blades grid 

Effective diameter 374 mm 
Blade height 22mm 
Chord of airfoil section 37,2 mm 
Pitch of grid 30,9 mm 
Number of blades 34 
Blade-angle setting 74?30* 
Outlet angle 14°50’ 


Fig. 2 Vane-blade and rotor-blade profiles as well as scheme of turbine 
stage flow part 


Object of the investigations — a turbine stage 


The schematic diagram of the flow part 
of the investigated turbine stage is shown in Fig. 2. 
Three-dimensional flow analysis of the invstigated 
turbine stage was performed by means of the CFD 
Fluent code in which the method of finite volumes is 
used for solving the set of Navier-Stokes equations. The 
calculations were performed for compressive, viscous 
(with constant viscosity coefficient) and turbulent flow 
without any heat exchange with environment. Geometry 
of the turbine flow part (Fig.2) was modeled by means of 
a computational mesh split into four blocks: inlet to the 
stage together with guide vanes, rotor, over-shroud seal, 
and outlet diffuser, which are so connected to each other 
as to form one computational entity (the applied kind of 
computational mesh and interfaces between particular 
blocks are described in detail in [16]). 


All the flow calculations were performed for the 
nominal clearance L = 0.5 mm and the eccentricity 
ew = 60% (vertical). 


Influence of type of seal on pressure forces 


To investigate influence of geometry of an isolated 
seal on forces generated in it, a simplified computational 
model was elaborated; in the model guide vanes were 
replaced with a smooth pipe and rotor blades grid was 
omitted. Fig. 5 shows cross-section through the upper part 
of the computational system geometry. The calculations 
were carried out for the constant value of working 
medium overpressure of at inlet to the pipe, equal to 5kPa. 
This value is equivalent (roughly) to the pressure at inlet 
to the seal (after expansion in guide vane), which occurs 
in the case of analysis of flow through the real turbine 
stage operating with the rotor angular frequency of 566 
rad/s and the inlet overpressure of 20kPa. 


Fig. 5 Cross-section through geometrical model of ‘isolated’ seal 


The run of pressure in the investigated seals along shroud 
breadth is shown in Fig. 6, and the pressure fields in the lower 
cross-section of seal, (D), and the upper one (G) —in Fig. 7. 
All the presented characteristics are in compliance with 
the Lomakin effect [5] according to which at inlet to 
unsymmetrical clearance the maximum pressure value 
occurs in the place where the lowest clearance value 
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Fig. 4 Schemes and meshes of the investigated seals as follows: 
a) with 4 teeth (marked U4), b) with single tooth at inlet (marked U1), 
c) with 2 teeth (marked U2), d) with single tooth in the middle of shroud 
(marked UIs), e) smooth one — without teeth (marked U0) 


In Fig. 3 an example of the applied computational 
mesh is presented, and in Fig. 4 are shown schemes, 
notations and computational meshes of the considered 
over-shroud seals, namely: 


* the seal with 4 teeth (marked U4), 
* the seal with single tooth at inlet (marked U1), 
* the seal with 2 teeth (marked U2), 


* the seal with single tooth in the middle of shroud 
(marked U1s), 


* the smooth seal, i.e. without teeth (marked U0) 
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Fig. 6 Pressure distribution along shroud breadth in the following seals: 
a) U4, b) U2, c) U1, d) UIs, e) UO, 
(f — arrangement of the cross-sections along seal circumference) 
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appears (cross-section D), and the minimum pressure 
value — in the place where the maximum clearance 
value appears (cross-section G). The resulting pressure 
difference causes the force acting in opposite direction 
to eccentricity, called also righting force, to appear. In 
the plane perpendicular to that direction the pressure 
distributions on both sides of the examined seal 
(cross-sections: L and P in Fig. 6f) are symmetrical, 
therefore no horizontal component of the resultant force 
appears. 


Also, velocity fields in the analyzed cross-sections 
of the seal appear to be different respectively. 
Velocity of working medium flowing through the seal 
reaches its maximum value in the place where the 
maximum clearance occurs (cross-section G), and its 
minimum value where the minimum clearance occurs 
(cross-section D). Direction of the flow is axial. The 
pathlines of the working medium flow through the seal 
UO are presented in Fig. 8. Red colour of the presented 
lines, which appears in the upper part of the seal, stands 
for the maximum flow velocity of working medium. For 
the sake of clarity, the flow through inlet pipe as well as 


Fig. 7 Pressure field in the upper cross-section G (of maximum clearance) 
and lower cross-section D (of minimum clearance) of the following seals: 
a) U4, b) U2, c) UI, d) Uls, e) UO 
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Fig. 8 Pathlines of the working medium flowing through the seal UO 
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pathlines in outlet diffuser have been omitted in the figure. 


Analyzing the obtained results one can concluded that 
type of applied seal decisively influences distribution of 
pressure field occurring in its clearance, consequently 
also pressure forces generated in it. In the same boundary 
conditions, values of the forces in particular seals differ 
several times or even a dozen or so times, to each other 


(Fig. 9). 


a) Pressure force Fc | —»- ML) 


b) 100% 
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Pressure force Fe/Femax [96] 


U1 
Type of seal 


Uts uo 


Fig. 9 Vectors, (a), and values, (b), of the pressure force Fc. occurring in 

the investigated seals (the force values are related to that which occurs in 
the seal U0) (red colour stands for the seal U4, green colour- the seal U2, 
black colour- the seal UI, violet colour- the seal Uls, blue colour- the seal UO) 


The Lomakin effect is observed in the cases in which 
shroud is motionless or working medium non-viscous. For 
viscous working medium and rotating shroud viscosity 
forces make pressure field in the seal clearance and 
location ofits extremes, changing. Pressure distributions 
in the plane perpendicular to eccentric displacement, on 
both sides of the seal (cross-sections: L and P), become 
no longer symmetrical, and, as a result, a force vector 
deviation from direction of eccentric displacement 
of rotor occurs. The pressure field distribution in the 
considered seals, for the case of viscous working medium 
and shroud rotating with the frequency of 566 rad/s, is 
presented in Fig. 10, and in Fig. 11 — the vectors and 
values of the calculated pressure forces. 


The greatest deviation of vector of generated forces, 
resulting from shroud wall rotation, occurs in the seals 
characterized by greater number of teeth (U2 and U4 
type). Moreover, in the seals friction effects due to motion 
cause a significant increase of generated pressure forces. 
For instance, in the seal of four teeth, U4, the force is 
increased almost ten times. In the remaining seals, smooth 
one, UO, and that of single tooth, U1, values of the forces 
do not undergo so significant changes as those observed 
in the motionless seal system. Even a small decrease of the 
force FC was observed in the seals of Uls and UO type. 
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Fig. 10 Pressure distribution in the seals along shroud breadth for viscous 
working medium and the shroud angular frequency equal to 566 rad/s: 
a) U4, b)U2, eI, d)Uls, e) U0 (f— arrangement of the cross- 
sections along seal circumference) 

In Fig. 12 flow through the four-teeth seal U4 is 
compared for the variants: when the seal shroud is 
motionless and when it rotates. In the latter case the flow 
through over-shroud clearance is no longer only axially 
directed because circumferential velocity components 
due to viscosity forces, appear. 


From additional numerical analyses of non-viscous 
flow through the examined seals no deviation of the 
direction of action of pressure forces from the eccentricity 
direction, resulting from shroud rotatary motion, was 
stated. Character of pressure field distribution in the seal 
as well as values of generated forces were in compliance 
with results of the analyses in the case of motionless 
shroud wall. 


Forces in the turbine stage with seals of various 
types 

Numerical investigations of operation of the turbine 
of 7.5kW output power at the angular frequency of 566 
rad/s (at 20kPa inlet working medium overpressure) 
were performed in order to make comparative analysis 
and determine resultant forces acting on the rotor in the 
real turbine stage fitted with labyrinth over-shroud seals 
of various types. Fig. 13 presents pressure distributions 
on the shroud of the examined seals, and Fig. 14 -vectors 
of the pressure forces FC, blade-ring forces FW and 
resultant force FP, acting on the rotor. As values of the 
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Fig. 11 Vectors, (a), and values, (b), of the pressure force FC occurring 
in the investigated seals for viscous working medium and the shroud wall 
angular frequency equal to 566 rad/s 
(the force values are related to that which occurs in the U0 seal of 
motionless shroud) 


b) 


Fig. 12 Pathlines of viscous working medium flowing through the seal U4: 
(a) in the case of motionless shroud, (b) in the case of rotating shroud with 
the angular frequency equal to 566 rad/s 


blade-ring forces are small they are shown in five times 
greater scale than the remaining — pressure and resultant 
ones. 


In Fig.15 are presented the values of the forces 
calculated for the investigated seals, related to the value 
of the resultant force FP occurring in the two-teeth 
seal U2. Blade-ring forces are marked blue, pressure 
forces — brown, and resultant force acting on turbine 
rotor — yellow. 
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Fig. 13 Pressure distribution in the seals: a) U4, b) U2, c) U1, d) Uls, 
e) U0 -for 20kPa working medium overpressure at inlet to turbine stage, 
and the rotor angular frequency equal to 566 rad/s 
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Fig. 14 Vectors of: the pressure force FC (on the left), blade-ring force FW 

(in the middle) and resultant force FP (on the right), acting on the rotor, for 

the model IV at 20kPa inlet overpressure and the rotor angular frequency 

of 566 rad/s (red colour stands for the seal U4, green — the seal U2, black- 
the seal Ul, violet- the seal UIs, blue- the seal U0) 
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Fig. 15 Values of the blade-ring forces FW (blue), pressure forces FC 
(brown), resultant force FP (yellow), at 20kPa inlet overpressure and 566 
rad/ s rotor angular frequency 


Analyzing the obtained results one can state that 
influence of seal design on blade-ring forces generated 


62 POLISH MARITIME RESEARCH, Special issue 2009/S1 


in rotor passages is small (comp. Fig.14 and Fig.15). For 
all the examined systems generated blade-ring forces are 
of comparably small values regardless of type of applied 
seal. 


Type of over-shroud seal decisively influences 
pressure forces generated in it. As results from analysis 
of the diagrams shown in Fig. 13, working medium 
pressure at inlet to and outlet from the seal is almost of 
the same value for all the compared systems, and it has 
been determined by values of the turbine stage operational 
parameters and eccentricity value. However type of 
over-shroud seal decisively influences pressure field 
distribution in its clearance (between inlet and outlet), 
and consequently pressure forces generated in it. 


As it turned out, the smallest pressure forces (Fig. 15) 
occurred in the seal of single tooth placed at inlet to 
the seal U1 and were nine times smaller than those in 
the seal U2 of two teeth, and over seven times smaller 
than those in the seal U4 of four teeth. This results 
from a favourable pressure distribution occurring in the 
clearance, determined by the conditions prevailing behind 
the rotor blade ring, and which is levelled out to a constant 
value over whole breadth of the shroud. Therefore no 
large circumferential differences in pressure appear 
(in contrast to the seals of a few teeth), and the generated 
forces are relatively small. 


In should be also observed that the pressure forces 
obtained in the real turbine stage differ from those 
generated in the isolated seals, i.e. without taking into 
account other elements of turbine stage. In the case of 
isolated seals the largest forces occurred in the smooth 
seal, UO, as well as that of single tooth placed in the 
middle, Uls. In the real turbine stage the largest forces 
are generated in the seal of two teeth, U2, the seal of 
smooth shroud, UO, as well as the seal of four teeth, U4. 
In all the performed analyses the seal of single tooth at 
inlet to the clearance, U1, appeared to be characterized 
by the smallest pressure forces. 


However it should be taken into consideration 
that the lowering of number of teeth in over-shroud 
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Fig. 16 Mass flux of leakages flowing through the seal U4, U2 and U1 at 
the rotor angular frequency of 566 rad/s and 20kPa overpressure of working 


medium at inlet (the values are related to the mass of leakage from the 
single-tooth seal U1) 


seal, favourable due to smaller values of generated 
aerodynamic forces, leads to the increasing of leak 
flux of medium flowing through seal clearance and 
in consequence to the lowering of total turbine stage 
efficiency, in the clearance. The working medium mass 
fluxes in the seal related to the largest leakage occurring 
in the single-tooth seal U1, are shown in Fig. 16. From 
the efficiency point of view the most favourable appears 
the four-teeth seal U4 which simultaneously - as results 
from the performed investigations - is, apart from that 
of two teeth, U2, the least favourable from the point of 
view of generation of pressure forces in it, which are able 
to trigger self-excited vibrations of rotor. Direction of 
action of generated pressure forces depends - as results 
from the previously performed investigations - a.o. on 
type of seal and number of teeth applied in it. For single- 
tooth seals (of U1 and Uls type) directions of action 
of pressure forces amount to about 38°, counting from 
eccentricity direction (Fig. 14). However for the two-teeth 
seal U2 such deviation from eccentrity direction amounts 
to 63°, and for the four-teeth seal U4 — as much as 80°. 
The smaller number of teeth in the seal the less deviated 
direction of action of pressure force from eccentricity, 
however the conclusion does not concern the smooth- 
shroud seal U0 because of another mode of expansion of 
working medium in a long and narrow clearance. 


Summary 


From analysis of the obtained results the following 
conclusions can be drawn: 


* The making use of force coefficients for isolated 
seals without taking into account the remaining 
elements of turbine stage can lead to very large 
errors in determining aerodynamic forces occurring 
in a real device. Therefore the coefficients 
determined for isolated seals should not be used in 
analyzing operation of a turbine stage. 


* For all the considered seals, direction of action of 
pressure force in the systems in which influence 
of working medium expansion in rotor blades grid 
or guide vanes has been not taken into account, is 
significantly different (even by 180?) from direction 
of action of pressure force in a real turbine stage. 


* For all the analyzed systems, pressure distributions 
in the planes perpendicular to turbine axis, at inlet 
to and outlet from seal clearance, are similar and 
determined a.o. by assumed boundary conditions, 
nominal clearance and eccentricity of the system. 
To a much lower extent they depend on design 
type of applied seal which decisively influences 
distribution of pressure field in its interior. 


The results presented in this paper are 
obtained from numerical calculations, however 
they have been confirmed by results of the 
investigations carried out on experimental stand. 
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Analysis of application of feed-water injector 
heaters to steam power plants 
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Roman Kwidzinski, 
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Abstract 


Steam-water injectors are devices in which exchange of mass, momentum and energy between two fluids being in direct contact, occurs. 
They can operate as pumps, mixers or direct contact heat exchangers. In the last aspect their use as feed-water heaters in Rankine 
thermal cycle of steam power plants both in land and sea applications (to merchant and naval ships) is very interesting. This 
paper presents selected results of heat-and-flow investigations of a supercritical steam-water injector, obtained in Institute of Fluid 
Flow Machinery, Polish Academy of Sciences (IMP PAN). On their basis value of average heat transfer coefficient for mixing 
chamber was determined; the obtained values were even a few dozens greater than those for classical shell-and-tube heaters. 
In the theoretical part of this work is presented an original injector model based on balances of mass, momentum and energy, 
written for control volumes containing separately particular elements of injector. On the basis of the model flow parameters in 
characteristic cross-sections of injector were determined. The calulations were performed for two different injectors tested in 
IMP PAN (Gdansk) and SIET (Piacenza, Italy), and their good compliance with experimental data was achieved. 


Keywords: power plants, feed-water heaters, supercritical steam-water injector, injector model based on balances of mass. 


Introduction 


Steam-water injectors are devices in which exchange 
of mass, momentum and energy occurs due to direct 
contact between two fluids. In 1901 an injector was used 
by Charles Parsons to air extraction from condensers, 
and in 1910 Maurice Leblanc applied an injector to 
a cooling device. Generally, injectors can operate 
efficiently as pumps, mixers or direct contact heat 
exchagers. In engineering a special role is given to double 
phase steam-water injectors. In recent years interest to 
such injectors has increased because of new proposals of 
their application to Rankine cycle of thermal power plants 
both land-based and marine ones (on merchant and naval 
ships). The most important advantages of such injectors 
consist in their possible operation without electric supply, 
rather small gabarites and lack of movable parts, resulting 
in their high operational reliability. 


Two-phase steam-water injectors can be generally 
split into two kinds: (a) -subcritical ones (Fig. 1a) and 
(b) -supercritical ones (Fig. 1b). In both the cases water 
steam is their driving medium. In subcritical injector 
steam expands in the tapered steam nozzle SN and its 
discharge velocity wy, is smaller than or at most equal 
to the sonic velocity a, (wy, <a), hence Mach number 
reaches value: M x 1. In the second case (Fig. 1b) 
water steam expands down to the pressure smaller than 
critical one, as a result the discharge velocity from the 
nozzle is greater than the sonic one, hence Mach number 
reaches value: M > 1. Discharge steam from the nozzle 
flows into the mixing chamber MC and contacts there 
with cold water. Due to difference of temperatures and 
velocities between both the phases heat exchange (steam 
condensation) and momentum exchange takes place. 
The processes terminate in vicinity of mixing chamber 
throat, that is manifested by a sudden rise of pressure, 
characteristic for shock wave. Behind the wave, only 
the liquid flows through the diffuser. From the cognitive 
point of view the steam-water injector exemplifies 
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a device which, though of simple design, is characteristic 
of a variety of flow phenomena. In the two-phase flow 
zone (covering the mixing chamber and condensation 
shock wave) drastic change in flow structure takes place 
— from stratified (annular) flow to homogeneous (bubbly 
one). Also flow velocity undergoes large changes — from 
a high subcritical value (smaller than about 450 m/s for the 
version a) or supercritical one (of the order of 500-1000 
m/s for the version b) at inlet to the mixing chamber, 
down to the value of the order of 20 m/s in the diffuser. 


In this paper the results are presented of the authors’ 
experimental investigations which covered measurements 
of pressure and temperature distribution along injector 
at different flow parameters, between inlet to and outlet 
from the injector. The investigations were carried out 
for the supercritical injector, (Fig. 1b), [1], i.e. the case 
considered more general. The critical injector (M = 1) 
is finally obtained by shortening the divergent part of 
the steam nozzle at maintaining the same inlet steam 
parameters. If steam pressure at inlet to the nozzle SN 
is now decreased then subcritical flows in the nozzle are 
obtained, and at the inlet: M « 1. 


On the basis of the investigations were determined 
the maximum values of the injector discharge pressure 
(Pamax) at which the device can operate stably. Also, 
value ofthe average heat transfer coefficient for mixing 
chamber as well as injector efficiency, were determined. 
The injector's theoretical model based on balances 
of mass, momentum and energy, written for control 
volumes containing flow through particular elements 
of the injector, i.e. steam nozzle, water nozzle, mixing 
chamber, condensation shock wave zone and diffuser, 
shown in Fig. 1, has been also presented. By solving the 
model’s equations for given parameters of driving steam 
and sucked-in water, flow parameters were determined 
in the injector characteristic cross-sections located on 
the boundaries of the above mentioned control volumes. 
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Fig. 1. Schematic diagram of steam — water injector: a) subcritical injector, 
b) supercritical injector. Main parts of injector : SN — steam nozzle, WN — 
water nozzle, MC — mixing chamber, DF — diffuser, t — throat 


Experimental investigations of steam-water 
injector 


In the Institute of Fluid Flow Machinery of 
Polish Academy of Sciences, Gdansk, experimental 
investigations of the supercritical two-phase steam — 
water injector driven by water steam (Fig. 1), were 
performed with the use of the laboratory model [1]. 
They were carried out by using the so called ,,short” 
injector of the mixing chamber pressure p „c =20kPa, with 
the aim to determine two quantities characteristic for 
injectors, namely: the maximum pressure p,,,, and 
value of the average heat transfer coefficient a for 
the mixing chamber MC. The tests were a continuation 
of the work carried out previously with the use of the 
so called ,,long” injector (of the pressure p,,. =8kPa) 
[2, 3] with the aim to investigate their applicability to 
feeding nuclear reactor in failure states. Inlet parameters 
of the tested laboratory injector could be controlled within 
the following ranges: steam flow rate between 85-130 
kg/h, inlet steam pressure down from 0.5 MPa, steam 
superheating temperature between 0+40 °C, water flow rate 
between 1+6 m?/h, water temperature between 15+20 °C. 
It was also possible to change conditions of two-phase 
flow through mixing chamber by controlling the water 
nozzle gap size within the range of 0.5+1.5 mm. 


During the steam-water injector investigations, 
apart from the above mentioned inlet parameters, also 
temperature and pressure distributions along mixing 
chamber and diffuser, were recorded. Moreover, 
value of counterpressure at outlet from the injector, 
possible to be set by means of control valve, was 
measured. In Fig. 2 typical distributions of pressure 
and temperature along mixing chamber and diffuser, 
measured at different values of discharge pressure 
(Fig. 2a) and water flow rate (Fig. 2b), are presented. 
During investigations of injector performance 
characteristics the maximum discharge pressure 
at different values of the injection coefficient 
U-m,m,, for three selected values of water nozzle 
gap, was determined. Results of the measurements 
are presented in Fig. 3. On the basis of the obtained 
characteristics it can be concluded that an increase of 


injector discharge pressure can be achieved both by 
increasing the coefficient U and mass flow of steam and 
or water at inlet to mixing chamber. 


a) 


z [m] 
b) 


TTTS 


P, = 200 kPa 


0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 
z [m] 


Fig. 2. Example profiles of pressure, (a), and temperature, (b), measured 

in mixing chamber and diffuser of a steam-water injector operating with 

selected values of the outlet pressure p4 and water flow rate G. The inlet 
steam mass flux. 


a) 


35 45 55 65 75 
U[-] 
b) 


Fig. 3. Values of the dimensionless maximum outlet pressure p g,,,,/P yo in 
function of the injection coefficient U, measured for three selected values of 
the water nozzle gap ô at the inlet steam flux equal to: 

a) 85 kg/h, b) 125 kg/h. 
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Heat exchange in mixing chamber 


The average heat transfer coefficient for the injector 
mixing chamber was determined on the basis of 
experimental data by using the method proposed in 
[3, 4]. By assuming that the flow of both phases through 
injector mixing chamber is one-dimensional and balanced 
it is sufficient, in order to determine thermodynamical 
parameters of water, to measure temperature only as 
steam parameters are depending only on pressure. Taking 
the above into account one can write the following 
energy balance equation for the mixing chamber: 


my hy, +m, hj, = (mj, +m, 12 (7, » (mj, =M, Jy (p;) (1) 


from which it results that the condensed steam mass flow 
is equal to: 


ii Hy, [hy -h (p, JH Mr [hu -h, (7, Jl Q) 


h, (,)- hp.) 


Now the average heat transfer coefficient for mixing 
chamber of the investigated steam—water injector can be 
calculated as follows: 


mh 


AucAT, 


MC MC 


(3) 


The so determined Q values were generalized and 
presented in dimensionless form with the use of similarity 
numbers resulting from dimensional analysis of the 
equations of mass, momentum and energy conservation, 
written for flow through injector [1,6]. The dimensional 
analysis resulted in a correlation function expressed in 
the form of Nusselt number dependent on dimensionless 
similarity parameters, for the average heat transfer 
coefficient. 


With the use of the authors’ experimental data for the 
two injectors which have been so far examined in IMP 
PAN [1,2,3] (comprising altogether 2051 measurement 
points from testing the ,,short” injector as well as the 
long” one) the following correlation function was 
obtained for the average Nusselt number: 


Nu = 0.00284 Re’? S 9.186 U91? E 0.605 (4) 
where: 


a D 
Nu = ——— 
Àz 
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The correlation (4) is valid for values ofthe dimensionless 
parameters contained within the following intervals, respectively: 


88240 < Re” < 234950, 24 < S' 211,11 xU € 76, 50<€<111 


together with the correlation coefficient R = 0.91. Physical 
parameters of water and steam as well as flow velocities in 
both the phases are calculated at the inlet cross-section of 
mixing chamber. 


In order to make use of the correlation (4) to know 
geometry of steam and water nozzles is necessary 
because pressure and velocities at inlet to mixing 
chamber are associated with it. In Fig. 4 are presented 
example calculation results of heat transfer coefficient in 
compliance with the correlation (4), for two values of 
steam mass flow rate and one value of water flow rate, 
which correspond with the conditions of the experiments 
carried out in IMP PAN. The points representing two 
values of the mixing chamber pressure Puc, equal to 
8 and 20 kPa, are depicted on the diagram. The steam 
flow velocities wy, as well as the steam nozzle outlet 
diameters D,, corresponding with the pressures were 
determined on the basis of thermodynamical calculations 
of expansion in steam nozzle (from the initial pressure 
Pyo to the final pressure Puc). As the water nozzle 
outlet area was known, the water flow velocity w,, 
and the slip ratio S could be determined. The term 
expressing the physical properties € was determined for 
the given pressure Puc from the water vapour tables. 
From the calculations two important observations result: 


* In mixing chamber the heat transfer coefficient a 
obtains very large values. 


e Values of the coefficient œ are increasing along 
with pressure increasing in mixing chamber. For the 
pressure Puc = 8 kPa its average value amounted to 
350 kW/m?K, and for the pressure p. = 20 kPa 
- as much as 725 kW/m?K. 


Such influence of the pressure on intensity of steam 
condensation on water film was confirmed by results of 
the injector investigations performed in Japan [12]. 


1.00E+6 + 
9.00E+5 + 
8.00E+5 + 
7.00E+5 + o 
6.00E+5 = ae —À 
5.00E+5 + a xut 
4.00E+5 + — 
3.00E+5 € 4277 
2.00E+5 Fe po pago 
68 10 12 14 16 18 20 22 
P uc [kPa] 


Fig. 4. Run of mean values of mean heat transfer coefficient, à, for mixing 
chamber, calculated acc. (4) in function of the pressure P gç for two 
selected values of the steam mass flux my, at a set value of the outlet nozzle 
diamater De, , and the water mass flux m,= 3000 kg/h.. 


LUI mV1 = 130 kg/h o" 
Hbc mV1 = 85 kg/h : 


E 


a [W/m K] 


Calculation of maximum discharge pressure 


Modeling the flow through steam-water injector faces 
great difficulties especially in the zone of two-phase 


flow in mixing chamber and condensation wave zone. 
The difficulties result first of all from the complexity of 
flow structures occurring successively in various parts 
of mixing chamber (annular, drop-like and bubble-like), 
presence of phase transformations (condensation on water 
film and in shock wave) and large changes in flow velocity 
(interphase slip at inlet to mixing chamber, transition from 
supercritical to subcritical shock wave). However it is 
possible to determine the maximum injector discharge 
pressure by using a simple model based on balances of 
mass, momentum and energy, written for control volumes 
containing mixing chamber and condensation shock wave 
zone, separately. 


In the proposed model equations which describe flow 
through steam nozzle, and water nozzle as well, can 
be solved for parameters given at inlet to injector, and 
consequently yield the flow parameters at inlet to mixing 
chamber (marked ‘1’). Next, the solution of the balance 
equations of mass, momentum and energy for two-phase 
mixture in mixing chamber provides data to solve the set 
of equations for the successive control volume (between 
the cross-sections ‘2’ and ‘3’), i.e. the shock wave zone. 
Now, on the basis of Bernoulli equation for liquid (water) 
flow through diffuser the pressure and temperature at 
outlet from injector, (marked ‘4’), can be determined. 


The equations which describe driving steam expansion 
in steam nozzle were obtained under the assumption 
that the flow is adiabatic and irreversible. The steam 
parameters at outlet from the nozzle (cros-section ‘1’) 
result from the following mass and energy balances: 


where s -index stands for isentropic expansion value 
(i.e. when s,, = Sy). In real flow irreversibility losses 
contribute to lowering the discharge velocity w,,,, that can 
be taken into account by applying the velocity coefficient 


Ayo PyoWyo = APM © 
2 2 2 
Ww Ww Ww 
hy + = ha, —u- = hp + © 
2 2 2 


c, defined as follows: 


Ww 
yl 
c, =—— U 
Wris 


The coefficient takes into account not only losses 
resulting from friction but also due to condensation 
(when expanding steam parameters exceed Wilson’s 
line) and formation of shock wave in the vicinity of the 
steam nozzle outlet (which takes place when at inlet to 
the mixing chamber a difference between steam and water 
pressure appears). As results from the investigations 
presented in [2], the above mentioned irreversible 
phenomena cause that value of the coefficient c, for 
injector steam nozzle is relatively low, i.e. equal to about 
0.9. 


Under the relation (7) the double equation (6) can be 
transformed to the following form: 


2 2 
Wris 2 Wy 
ha = hys +( = c) 3 = chy, + ( = Zu + 3 | din 


In the typical steam-water injector the driving steam 
is superheated at inlet to the nozzle but during expansion 
it enters wet-steam zone, hence the balance equations 
make the below given constitutive relationships to be 
complete: 


Pro = P (Pros hys) > Syo =8(Pyaslyo) (8a) 
hy, = Xi My sat (pri )+ ( a Xi ui (pri) (8b) 


x us Syo — SL sat (pri) 
T Sy sa (Pri) Sirsa (Pri) 


l- x 


X 
Py, = + 
" PY sat (pr) P L sat (pr) 


qm hy, E (pri) ee 
hy sat (py, ) hy sat (py, ) 


where x stands for steam quality and the index sat 
concerns parameters in state of saturation. 


(8c) 


(8d) 


Water nozzle discharge parameters can be 
determined under the assumption that pressure of water 
is the same as that of steam: py, — p,, — p,, and its 
temperature is the same as that at inlet to injector. Then 
PLi = P(Pvi Tio) and the mass conservation equation: 


Mro = Ay PuWi 9) 


is sufficient to determine the water nozzle discharge 
velocity w; ,. 


Flow parameters on the mixing chamber boundaries 
satisfy the balances of mass, momentum and energy, 
counted between the inlet cross-section 1 and the throat 
2 (see Fig. 1), which are expressed by the following 
equations: 


4,p,W, = Ay py Wy, + App; Ww, Or m, =m m,00 


m,w, + A, p, t (A, — Aj) Puc = my Wy, t mj wy + Ap, an 


2 2 2 
mh x = ZD 23 DD 2: i 


It was assumed here that steam phase is always in state 
of saturation and in the cross-section 2 both the phases 
(steam and water) move with the same velocity w,. 
Therefore to determine the enthalpy h, the relationship 
valid for homogeneous mixture can be used: 


h, = Al pT) -xJhy((p,t-(ü-x)h.((pT,) a» 


Moreover, it was assumed that the average pressure 
which acts on the mixing chamber conical surface and 
appears in the momentum balance equation (11), is equal 


to pyc = (P,+p,)/2. 


When the mixing chamber inlet parameters (cross- 
section 1) are known to solve the set of balance equations 
(10)=+(12) is still necessary in order to know the liquid 
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phase temperature in the throat, T,,, which is lower than 
the steam temperature Ty. The temperature T}, can be 
determine on the basis of heat flow which rises water 
temperature during condensation of driving steam in 
mixing chamber. The heat flow is connected with the 
average heat transfer coefficient Œ by the relation (4). 


The balance equations for the shock wave zone 
(between the cross-sections 2 and 3) are similar to those 
for mixing chamber and take the following form: 


A,p,W; =A,p,w, or m,=m, (14) 


mW; + Ap, + (4, — Aj) psy = mw, + Ay Py a5) 


2 2 


W. W. 
m, RES =m, as (16) 


In the shock wave zone complete steam phase 
condensation takes place, hence in the cross-section 3 
only liquid (water) appears, whose enthalpy amounts to 


h, = h(p3,p3). 


The average pressure which acts, within the 
shock wave zone, on the diffuser conical surface, 
was determined from the approximate relation: 
Psw = (p5*p3)/2. Shock wave width and location of the 
cross-section 3, resulting from it, was determined on 
the basis of pressure profile measurements. 


On the basis of the balance equations (14)+(16) 
e.g. water pressure in the cross-section 3 can be 
determined. By taking into account the pressure 
increase determined form the Bernoulli equation: 


+ Path = pu (- D 3 


—.w a7) 
A, 


Pat 


where p, =P, W, = 


value of the injector discharge pressure p, can be 
achieved. Value of the pressure loss coeficient ¢ in (17) 
can be determined acc. [5]. As in the equations (14)-(16) 
it was assumed that shock wave is located in direct 
vicinity of mixing chamber throat the determined 
pressure p, obtains its maximum value to keep injector's 
operation stable. 


The presented model was used to determine the 
maximum discharge pressure for the flow conditions 
complying with the measurements performed in IMP 
PAN, [1], and presented in Fig. 3b, as well as with the 
investigations carried out in SIET (Piacenza, Italy), [7,8]. 
Relatively good calculation conformity was obtained 
between the above mentioned experiments (Fig. 5 and 6), 
namely with the error of about 15%. 


Injector efficiency 


In the subject-matter literature different definitions 
of injector efficiency can be met [9,10]. It is usually the 
ratio of compression work done within injector and steam 
energy delivered to the injector. As steam-liquid injector 
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Fig. 5. Results of calculations of the maximum outlet Pressure Pamay, for 
the flow conditions given in Fig. 3b: a) calculated dimensionless maximum 
outlet pressure, b) relative difference between the calculated value, Pyyqy + 

and measured one, p exp. 
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Fig. 6. Measured and calculated values of the pressure p and water 
temperature T on wall of mixing chamber and diffuser : the left column (a) 
concerns IMP- PAN conditions, the right one (b) concerns SIET conditions 

(Piacenza, Italy) [7,8]. 


can operate as a heat exchanger or jet pump therefore 
acceptance of two definitions of injector efficiency is 
justified. 


In the case of injector operation as a heat exchanger it 
is proposed to use the notion of the exergetic efficiency n, 
defined as the ratio of exergy of working fluid discharged 
from injector and exergy delivered to it: 
(myo + mro ) b4 -(1 We 
Tig byo + mro bro 


(18) 
= 


Calculation results of the efficiency for two versions 
(,,long” and ,short") of the injectors examined in IMP 
PAN are presented in Fig. 7a. 


In the case of injector operation as a jet pump, water 
exergy after compression is not so much important as its 
pressure part [11]. Hence the jet pump efficiency np can 
be defined as the ratio of compression mechanical work 
and delivered exergy: 


_(mo+tmro) |— Pa- Puc 
PL My byo + M9 bro 


_1+U P4-PwC ag 
Pz bwo tU bro 


p 


Calculation results of the jet pump efficiency n, are 
presented in Fig. 7b. From comparison of the efficiency 
N, with the np it can be observed that the first is much 
greater. 
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Fig. 7. Efficiency calculated for two versions of the injector experimentally 
tested at IMP PAN: a) exergetic efficiency, b) injector pump efficiency. 


Final remarks 


In this paper the comprehensive results have been 
presented of model tests of two-phase steam-water 
injector, aimed at its application to steam power plants 
operating in compliance with Rankine cycle. In such 
cycle the so called regeneration heating of water feeding 
the boiler takes place within shell-and-tube heater. It is 
planned to use steam-water injectors instead of shell-and- 
tube heaters, that could significantly lower investment 
cost of power plant. 


The model tests of two-phase steam-water injector 
performed in IMP PAN were aimed at determination of 
the coefficient of heat transfer between steam and water 
in mixing chamber as well as the maximum pressure of 
water discharged from the injector. 


Inthe theoretical part ofthe work the injector's model 
based on the balance equations of mass, momentum 
and energy, written for particular injector elements, 
was elaborated. In the model correlation of heat 
exchange within mixing chamber was used to determine 
temperature difference between steam and water in the 
mixing chamber's throat. The model makes it possible to 
determine, apart from flow parameters in characteristic 
control cross-sections, also the maximum pressure p,,.,. 
at outlet from the injector. Correctness of the model, in 
the aspect of the pressure, was examined with taking 
into account the authors' experiments as well as those 
performed in SIET (Piacenza, Italy). Good comformity 
between results of the calculations and experiments was 
achieved as the error amounted to about 15% within 
a wide range of flow paramters at inlet to injector. 

The experimental test results of heat exchange in 


injector mixing chamber confirmed that to achieve 
very high values of the heat transfer coefficient ©, 
reaching even 700 kW/m? K, is possible. They are over 
100 times greater than the heat transfer coefficient k in 
shell-and-tube heaters (which reaches the value of about 
3 kW/m? K). The so large values of @ could make it 
possible to decrease, by about a hundred times, surface 
area of steam-water contact in mixing chamber (at the 
same value of transferred heat flow rate), that could 
potentially result in the decreasing of gabarites of an 
alternative heat exchanger in the form of injector, by about 
10-times. Hence, such application of the injector could 
significantly lower investment cost of steam power plant. 


The lowering of gabarites of steam power plant devices 
is especially important in case of merchant and naval 
ships where shortage of space takes place. By applying 
injector water heaters their number could be increased, 
that could result in increasing Rankine cycle efficiency 
and -in consequence -lowering operational costs. Another 
factor which leads to improving operational parameters 
of the cycle fitted with feed-water injectors is lack of 
thermal degeneration. In injectors, in contrast to shell- 
and-tube heaters, heat transfer occurs by direct contact of 
vapour phase with liquid one, hence in such devices heat 
transfer effectiveness is not worsened due to accumulation 
of contaminations on heat transfer surface during their 
long-lasting operation. Preliminary estimations indicate 
that the increase of Rankine cycle thermal efficiency by 
about 0.2 96, resulting from the above, can be expected. 
Another gain can be achieved due to injector pump effect 
which makes it possible to decrease feed-water pump 
power by a few dozen percent. 


This scientific research project (the ordered project 
No. PBZ-MEiN-4/2/2006) has been financed from 
science budget resources for the years 2007-2000. 
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Nomenclature 
A — cross-section area, m? B 
a — sonic velocity, a = (dp / dp dus rs : 
Awc - surface area of mixing chamber wall, m? 
b — specific exergy, kJ/kg; 
D,d - diameters, m; 
G — volumetric water flow rate, m?/h; 
m — mass flow rate, kg/s; 
M — Mach number, M = w/a; 
Nu  -Nusselt number, Nu 2 ad, / À y; 
h — specific enthalpy, kJ/kg; 
hfg — heat of condensation, kJ/kg; 
p — pressure, Pa; 


Re — Reynolds number, Re=wpd,/u; 
S — thermal capacitance, kJ/kg K; 


S — slip coefficient, S = wy, / wr; 
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T — temperature, K; 

AT,.  —logarithmic mean temperature in mixing 
chamber, K; 

U — injection coefficient, U = m, o/Myp; ; 

w — flow velocity, m/s; 

x — steam quality, x = my /( my +m; ); 

a — mean heat transfer coefficient for mixing 


chamber, W/m2 K; 
— water nozzle gap, mm; 
— thermal conductivity coefficient, W/(m*K); 
— injector efficiency; 


— dynamic viscosity, kg/(m*s); 


ow E 3s m» o0 


— parameter of physical properties, 
0.5 0.1 
$-(Pn/pyi) "(Mi/up) 
p — density, m3/kg. 
lower indices: 
— at injector inlet; 
— at inlet to mixing chamber; 


— at mixing chamber throat; 


0 
1 
2 
3 — behind shock wave; 
4 — at outlet from injector; 
b 


— (stands for) exergy; 


e — equivalent; 

max -— maximum; 

L — (stands for) liquid (water); 

p — (stands for) pump; 

SN  -(stands for) steam nozzle; 

V — (stands for) vapour (of water). 
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On the modelling of aerodynamic force 
coefficients for over-shroud seals of turbine 


Wojciech Włodarski, 
Krzysztof Kosowski. 


stages 


Abstract 


This paper presents experimental investigations which made it possible to determine dynamic coefficients of labyrinth 
over-shroud seal of a model air turbine. The coefficients associate pressure forces with turbine rotor displacement, velocity 
and acceleration respective to turbine casing (linear model) and play important role in analyzing turbine-set dynamics. 

The obtained results indicated that involving serious errors can be expected in the case of application of the 
simplification consisting in neglecting inertia coefficients, proposed in the literature. It was simultaneously demonstrated 
that seals can be also met of weak damping qualities, for which to neglect damping coeffcients is allowable. 


Keywords: rotor dynamics, self-excited vibrations, aerodynamic forces, blade seals, turbine shrouds 


Introduction 


In engineering many cases are known of occurrence 
of aerodynamically generated vibrations of turbine 
rotors. Such phenomena can lead to dangerous vibrations 
which make correct operation of turbine set impossible. 
Self-excited aerodynamic forces occur as a result of 
phenomena associated with non-stationary flow of 
working medium through a turbine stage in the case 
of eccentric position of turbine's rotor respective to its 
casing. During operation of a rotary machine whose 
rotor axis is deviated from casing axis, local flows through 
particular passages between rotor blades are changed as 
well as working medium leakages along circumference 
of radial gap of over-shroud seal are also changed. As 
a result, a circularly unsymmetrical distribution of forces 
acting on rotor blades and shroud surface occurs. Hence in 
publications devoted to the problem two main components 
of aerodynamic forces resulting from different character 
of their occurrence, are distinguished: 


e resultant force due to non-uniform distribution of 
circumferential forces in rotor blade ring, called 
blade-ring force, 


e resultant force due to non-uniform distribution 
of pressure in radial gap over rotor blades, called 
pressure force. 


In most cases of design solutions of turbine stages the 
pressure forces are of decisive influence on occurrence of 
self-excited rotor vibrations of aerodynamic origin. For 
this reason to know dynamic features of over-shroud seals 
in designing process of a new turbine set is very important. 
The features are tried to be determined in a theoretical 
way, however results of experimental investigations 
appear extremely valuable because of a highly complex 
character of flow through turbine stage. 


Pressure forces — linear model 


Mutual relationship between pressure force and 
displacement of rotor against casing is usually described 
by means of the so called ,,seal dynamic coeffcients" 
[1-8, 12-14] which can be incorporated into whole model 


of turbine set during its designing process. Three kinds of 
the dynamic coefficients can be distinguished: stiffness, 
damping nad inertia ones. They are associated with 
displacement, velocity and acceleration of rotor relative 
to casing, respectively. 


In such system, forces and moments generated in radial 
gap over turbine rotor blades are described by using the linear 
model which takes form of the following vectorial equation: 


x x X 
y y y 

Fk] Hel; MEZ | œ 
Ey E, Ey 


where: 


F — vector of aerodynamic forces and 
moments generated in seal, 


x,y — rotor linear displacements in the 
assumed frame of coordinates, 
ex, £y — angular displacements of rotor axis, 
K —matrix of stiffness coefficients, 
C — matrix of damping coeffcients, 
M - matrix of inertia coefficients, 


- og g  — the first time-related derivatives of 
X, y,£ ,,€, : : 

linear and angular displacements, 
X.y.g .g, — the second time-related derivatives 
X, y 20 xy» y . , 

of linear and angular displacements. 
In the subject-matter literature a simplification which 
consists in neglecting the inertia coeffcients in the model, 
can be found [5,6], [8+11]. The vectorial equation of the 
linear model takes then the following form: 


x X 
y y 
Fl-k]_ +c]; (2) 
£c E 
£, £, 
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Research stand 


The above mentioned investigations were performed 
by using a research stand of model turbine installed 
in the Laboratory of the Ship Automation of Turbine 
Propulsion Department, Faculty of Ocean Engineering 
and Ship Technology, Gdansk University of Technology. 
Schematic diagram of the stand is presented in Fig. 1. 
The turbine is fed with compressed air delivered from the 
compressor through the piping system. The compressor 
is driven by the electric motor through the reduction gear. 
Th electric brake absorbs the turbine’s output power. 
The electric energy generated in the electric brake is 
transformed into heat in the set of resistors. 


The examined turbine is a single-stage impulse turbine 
with labyrinth over-shroud seal having two teeth in the 
shroud and casing (Fig. 2). Construction of the turbine 
makes it possible to change mutual position of casing and 
rotor during operation of the stand. 


The measuring system of the stand makes it possible 
to measure non-stationary distribution of pressure 
along circumference of the over-shroud seal. The 
turbine’s casing is adjusted to fastening the ENTRAN 
EPE miniature pressure sensors (Fig. 2). Altogether 88 
pressure measurement points are provided for: 4 in each 
row, circumferentially distributed every 150, except for 
the parting plane of the casing. 


During the carried-out investigations a displacement 
of the turbine’s casing (precisely, of the vane 
blade-ring with seal ring) against the turbine’ rotor, was 
forced (during operation of the turbine). The casing 


Fig. 1 Schematic diagram of the model air turbine stand: 1 — turbine, 
2 — brake, 3 — compressor, 4 — reduction gear, 5 — electric motor, 
6 — set of resistors 
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Fig. 2 Schematic diagram of the measurement system of pressure in over- 
shroud gap of the model turbine: 1 — axial cross-section of the turbine stage, 
2 — pressure sensor, 3 — measurement amplifier, 4 — recording system of 
measurement data, I, II, III, IV — measurement planes perpendicular 
to rotor axis 
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displacement was measured by using an induction sensor 
of linear displacements. The sensor was fastened to the 
motionless part of the stand and its measuring tip was 
leaned against outer surface of the movable casing of 
the turbine. 


Experimental investigations 


The research was aimed at verifying validity of the 
simplification of linear model, consisting in neglecting the 
inertia coefficients (Eq. 2), as proposed in the literature. To 
this end the dynamic coefficients ofthe over-shroud seal of 
the model air turbine were experimentally determined in 
compliance with the following vectorial equation: where: 


2 


di ds 
Fy(0) 2 Kyy Kyy | x(t) 4 Cix Cy | dt}, My Myy | d£ (2) 
FQ} |Kyx Kyy [yO] [Cyx Cyy] | W| | Myx Myy | dy 

dt dt 

Fy, Fy — components of pressure 
force vector, 

— time, 

X,y — rotor displacements in 
the assumed frame of 
coordinates, 

Kyy, Kyy, Ky, Kyy — stiffness coefficients, 

Myx, Myy» Myy, My, X — inertia coefficients, 


Cx: Cy. Cyy, Cyx 


— damping coefficients. 


During the investigations changes of pressure in 
the sealing, which result from mutual transverse 
displacement of rotor and casing, were recorded; next 


on the basis of the experimental data pressure forces 
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Fig. 3 Displacement of rotor against casing as well as change of pressure 
in the over-shroud seal, resulting from it and measured in 4 points 
distributed every 900 along circumference of the turbine stage; a. — angle 
which determines location of measurement point along circumference of the seal 
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Fig.4 Increase of the resultant pressure force at the step change of mutual 
position of rotor and casing by 0,74 mm and different turbine s loads 


and dynamic coefficients of the sealing were calculated. 
The investigations were carried out at the constant rotor 
speed equal to 3600 rpm, and the turbine power values 
of 2 kW, 4 kW, 6 kW, 8 kW, 10 kW, 12 kW, 14 kW, 
16 kW, as well as the values of eccentricity of rotor 
against casing, changed by 0,16 mm, 0,35 mm, 0,61 mm, 
and 0,74 mm. The nominal radial gap in the seal was 
equal to 1 mm (rotor co-axially aligned in casing). 


In Fig. 3 is presented the example run of position of 
rotor against casing as well as that of change of pressure 
in the over-shroud seal, resulting from it and measured 
in 4 points located every 90° along circumference of the 
turbine stage, at the turbine’s power of 12 kW and its 
rotor speed of 3600 rpm. 


On the basis of the measurements the forces acting on 
the turbine’s shroud were determined; the example of the 
obtained results is presented in Fig. 4. 


As pressure force changes with time and rotor 
displacements against casing are known, it was possible to 
determine the stiffness, damping and inertia coefficients 
appearing in the linear model equations. The coefficients 
were determined by searching for solution of Eq. 3 with 
the use of the least squares method. The example results 
are shown in Fig. 5, 6, and 7. As observed on the basis of 
the obtained results, the relationship between the dynamic 
coefficients of the considered over-shroud seal and the 
turbine’s output power, is linear. 
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Fig. 5 Stiffness coefficients of the seal in function of output power of the 
turbine 
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Fig. 6 Damping coefficients of the seal in function of output power of the 
turbine 
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Fig. 7 Inertia coefficients of the seal in function of output power of the 
turbine 
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Fig. 8 Pressure force calculated by means of the linear model and that 
calculated by means of the model neglecting inertia coefficients. 


By making use ofthe determined dynamic coeffcients 
of the seal it was possible to determine pressure force 
for a given value of displacement of turbine's rotor 
against its casing. The example run of the pressure 
force calculated by means of the linear model with 
the use of the determined coefficients is presented in 
Fig. 8. Inthe same figure is also drawn the similar run but 
determined with the use of the model neglecting inertia 
coefficients. 


The diagram in Fig. 8 shows that the neglection 
of inertia coeffcients in the linear model equations can 
lead to distinct differences in runs ofthe calculated force. 
In the presented example the differences reached as far 
as 15 96. 


On the other hand, rather small influence ofthe inertia 
coeffcients on the determined forces has been revealed. In 
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Fig. 9 Pressure force calculated by means of the linear model and that 
calculated by means of the model which neglects damping coefficients. 


Fig. 9 are presented the example runs ofthe pressure force 
calculated by means ofthe complete linear model, i.e that 
in which all three kinds ofthe coefficients have been taken 
into account, and that calculated by using the simplified 
model neglecting the term associated with damping 
coefficients -they differ only a little to each other. 
Such situation repeated within all the range of values of 
turbine's load and change of rotor displacement against 
casing, for which measurements have been performed. 
This suggests that the considered type of seal (i.e. the 
labyrinth over-shroud seal with two teeth on the shroud 
and two on the casing) is of very low damping qualities. 


Therefore in the case of this kind of seal it was 
proposed to simplify the model which describes pressure 
forces, by omitting its part associated with velocity 
of rotor displacement against casing. In such case the 
equation which describes generation of pressure forces 
would take the following form: 


Conclusions 


(4) 


*in the case of labyrinth over-shroud seal, values 
of dynamic coefficients of the seal change 
proportionally along with turbine's power changing; 


ein the case of seals of low damping qualities, 
omittance of inertia coefficients can lead to noticeable 
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errors in determining the aerodynamic forces; 


*the labyrinth over-shroud seal with two teeth on 
the shroud and two on the casing has as much 
low damping qualities that in the case of the seal 
in question the linear model can be simplified 
by omitting the term associated with damping 
coefficients. 
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